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INTRODUCTION 


I IS WELL KNOWN in the theory of elasticity that the 
problem of uniform twist of a cylindrical column 
possesses a simple solution. This solution, commonly 
known as the Saint-Venant solution, gives a system of 
strains and stresses uniform along the axis. Moreover, 
the axial stresses vanish throughout the column and the 
shearing stresses arising from the uniform twist produce 
a warping of the cross section uniform along the axis. 
This rather simple situation is absent if the warping is 
prevented or limited in some way, or the twist is not 
uniform. ‘Then the strains and stresses will vary from 
section to section, and there are axial stresses in addi- 
tion to the shear stresses. Thus the problem becomes 
extremely complicated. 

An approximation is possible if the variation of twist 
along the axis of the cylinder is small, i.e., the twist is 
almost uniform along the axis. Then as a zero order 
approximation, the warping of the section can be cal- 
culated from the Saint-Venant theory corresponding 
to the local value of twist at the section considered. 
The axial rate of change of the warping so calculated 
gives an induced axial stress. ‘This induced axial stress 
will then serve as the starting point of a first order ap- 
proximation of the so-called ‘“‘secondary shear stresses”’ 
as a correction to the Saint-Venant solution., In mathe- 
matical language, this method of approximation is an 
iteration process. S. Timoshenko,? H. Wagner,* R. 
Kappus,‘ J. N. Goodier,’ and others have applied the 
method for determining the torsional stability of open 
thin-walled channels under axial compression. Their 
results agree well with the experiments and the treat- 
ment is indeed satisfactory as far as this particular type 
of structure is concerned. 
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However, if the problem to be solved is that of torsion 
of a closed thin-walled cylindrical column with one end 
built in, such as the twisting of an airplane wing, the 
convergence of the iteration method is doubtful. In 
fact, as shown in a previous analysis by the senior 
author and Christensen,® the secondary stresses com- 
puted in the above-mentioned manner are considerably 
greater than the primary shear stresses. Therefore, 
it is necessary to re-examine this problem to obtain a 
satisfactory solution. 

From the viewpoint of the theory of elastic bodies, 
such a thin-walled cylindrical column should be con- 
sidered in general as a thin shell, or if the wall thickness 
is extremely small in comparison with the section di- 
mensions, as a membrane. A thin shell possesses local 
stiffness with regard to bending and twisting of the mid- 
dle surface as well as to tension, compression, and shear 
in the middle surface, whereas the bending and twisting 
stiffness of a membrane can be neglected. Hence, the 
stress system ascribed to a membrane is not complete 
and can only be in equilibrium under certain specific 
conditions. It is easy to see that in the absence of ex- 
ternal forces normal to the surface of a cylindrical mem- 
brane, no circumferential stresses can exist, since they 
cannot be in equilibrium because of the curvature of 
the membrane. Now, in the case of uniform twist 
without axial restraint, no such stresses occur, and the 
structure can be considered as in equilibrium under a 
membrane stress system without further complications. 
However, this is only an exceptional case. In general, 
the problem of nonuniform twist of a thin-walled cylin- 
drical column should be treated by the general theory of 
cylindrical thin shells with the complete system of 
stresses. 

Unfortunately, the exact theory of thin shells in- 
volves, in general, mathematical complications that 
restrict its application to geometrically simple cases. 
However, in many practical applications—for example, 
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in many airplane wing constructions—most of the force 
normal to the shell is resisted by internal reinforce- 
ments, such as ribs or webs. For such cases one can 
develop an approximate theory of torsion using the as- 
sumption that a web system is continuously distributed 
and has perfect rigidity against distortion of the section 
-nits plane. However, axial warping, or axial displace- 
ment of the section is not prevented. In other words, 
the displacement of the section due to the twist of the 
structure can be considered as a rigid body rotation 
around the shear center in its plane together with a cer- 
tain amount of axial warping. The displacement can 
be expressed in two components, the rotation and the 
warping. The bending moments of the shell need not 
be considered, as the normal forces necessary for equilib- 
rium with the membrane stresses are supplied by the 
webs. In this case, only the condition of equilibrium 
in the axial direction enters into the problem; the other 
two equilibrium conditions are satisfied automatically 
by taking into account the reaction of the web system 
acting in the plane of the section. It is believed that 
this assumption yields a fairly good approximation for 
a thin-walled cylindrical structure having a consider- 
able number of webs and ribs. Such a thin-walled 
structure stiffened with an internal web system under 
nonuniform twist will be studied in this paper. 
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Fic. 1. Coordinate system. 
NOTATIONS 


The following notation is used (Fig. 1): 


2,3 = the linear coordinates parallel to the axis and along 
the centerline of the curved wall, respectively 
r = the normal distance between the center of rotation 


of the cross section and a tangent to the center- 
line of the wall at an arbitrary point . 


t = the wall thickness which is assumed to be constant 

Gr, 03,7 = the stress components acting on the faces of a wall 
element 

€z, €s, Y = the corresponding components of strain 

u = the axial displacement of an arbitrary point of the 


wall centerline, or warping function 
= the rate of twist 
7" = Young’s modulus and Poisson’s ratio, respectively 
= shearing rigidity 
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I = moment of inertia about the center of symmetry of 
the cross section 
M = total twisting moment 


FUNDAMENTAL EQUATIONS 


Consider a thin-walled cylinder of constant cross 
section twisted by an external moment J/ under various 
possible end conditions, It is required to find the rate 
of twist and the warping of various sections. Through- 
out the following treatment, two assumptions will be 
made: 

(a) The bending stiffness of the thin wall is ne- 
glected. 

(b) The deformation of an arbitrary section con- 
sists of a rigid rotation in the plane of the section and a 
warping, i.e., displacement parallel to the axis of the 
cylinder. 

According to these assumptions, the stresses—the 
stress components—parallel to the axis are in equilib- 
rium, whereas the strain in the circumferential direction 
vanishes. The zero value of the strain in the circum- 
ferential direction is enforced by the web system, which 
is supposed to be infinitely rigid and continuously dis- 
tributed along the beam. Hence, the stress compo- 
nents acting in the circumferential direction on the faces 
of an arbitrary element of the thin wall are balanced by 
forces from the webs acting on the inner surface of the 
wall, 

These assumptions yield two conditions for the de- 
termination of the stress distribution: the equilibrium 
condition in the x direction and the condition of vanish- 
ing strain in the circumferential direction. The 
equilibrium condition is expressed as 

(Oc,/0x) + (07/O0s) = 0 (1) 
The other equilibrium conditions are satisfied by taking 
into account the forces from the webs. In fact, these 
additional equilibrium conditions can be used to com- 
pute such web forces, if necessary. The condition of 
vanishing strain in the circumferential direction can be 
written as 


e, = (1/E)(o, — woz) = 0 (2) 


The two other strain components e, and y are given by 
the stress-strain relations 





ated (1/E) (oz = os), a. 1/G (3) 
From Eqs. (2) and (3), it follows that 
E E 
Cee eee Seek ee (a) 
Li-, : 2(1 + p) ’ 


Substituting the values of ¢, and r in Eq. (1), one ob- 
tains : 

2 Oe . OY 

——-+-—— = 0 5 
1 — 2O, a Os (5) 





Now ¢, = 0u/Ox and according to Eq. (2) 
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TORSION 


y = (0u/ds) + 18 (6) 


By introducing these values in Eq. (5), the following 
fundamental equation is obtained: 
2 & 0? d 

ro a « if -. (7) 

1 — p Ox? ~~ Os? ds 
On the right side of this partial differential equation, 0 
is a function of x only, whereas dr/ds is a given function 
of s. 

The relation between the stress and strain in the 
cylinder and the external moment M is applied by inte- 
grating the moments rr ds about the shear center due to 
each element ds. The forces between the web system and 
the wall do not contribute to the resultant mo- 
ment. 





t$rrds = M (8) 
Or by substituting + from Eqs. (4) and (6) 
t$(0u/ds)r ds + I0 = M/G (9) 


Thus where J is the moment of inertia of the section 
about the center of shear, i.e., 


I=tfr'ds (10) 


Eqs. (7) and (9) with appropriate boundary conditions 
determine the solution of the problem. For a closed 
section with continuously varying r(s), Eq. (9) can be 
transformed by integration by parts into the following 
form: 


—t§u(dr/ds) ds + I8@ = M/G (11) 
Eliminating 6(x) from Eqs. (7) and (11), one obtains 
the equation for u, 


10% 0% . tdr M dr 

pastas tg, EGS 
This is a linear integrodifferential equation for u(r, s), 
where for simplicity, k? stands for (1 — y)/2. 


=0 (12) 


POLYGONAL SECTIONS 


Consider first the case of a box section consisting of 
a number of flat sheets. Then, for each side, dr/ds = 
0, while at the corner r changes discontinuously. 
Hence, for every point except for the corner points, u 
satisfies the equation 


(1/k?)(O%/Ox2) + (0%u/ds%) = 0 (13) 


This equation can be reduced to the Laplace equation 
by introducing § = kx as a new variable. One obtains 
by integration of Eq. (7) around the corner point 


A(du/ds) = 


where the symbols A(0u/0s) and Ar denote the dis- 
continuous jumps in the values 0u/0s and r. 

For this case, Eq. (11) yields for the rate of twist @ 
the expression 


—6Ar (14) . 
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6 = (M/GI) + (t//I) © uA (15) 
t-1 
where u, .... u, and Ayr .... A,r are the values of u 


and Ar at the m edge points. Combining Eqs. (14) and 
(15), one obtains the edge conditions for u: 


a($*) = -= Ar — Ar- >> u,Ay (16) 

It is seen that this equation gives linear transition 
conditions between the jumps in the values of the func- 
tion u and its derivative 0u/0s at the edge points. 

For equilateral polygons, where Ay = 0, the right- 
hand side of Eq. (14) vanishes, i.e., not only the warp- 
ing function u(x, s), but also its first derivative, are 
continuous at all the corner points, and consequently 
u(x, s) is continuous and regular everywhere in the 
region. Hence, u(x, s), satisfying Eq. (13), is a har- 
monic function periodic in one of the variables, s. If 
the cylinder is restrained at one end and free at other, 
ie., u(0, s) = 0, (Ou/Ox),., = 0, then u(x, s) can be 
shown to vanish everywhere as follows: The harmonic 
function u(x, s), periodic in s and vanishing at x = 0, 
must be of the forms 

A, sin = sinh — aT B, cos — = sinh = k 

l l l l 
where » is any integer, and / is the circumferential length 
of the section. By the free end condition at x = L, 
i.e., (Ou/Ox)2., = 0, A, and B, vanish. Thus, u(x, s) 
= 0. This shows that for an equilateral polygon the 
warping function vanishes everywhere and the uniform 
twist rate is given by Eq. (15), 


6 = M/GI (17) 
This is well known from Saint-Venant’s theory. 


TORSION OF A THIN-WALLED COLUMN OF A 
RECTANGULAR Box SECTION 


The application of the method described above will 
be evident by the following simple example. 

Consider a rectangular box section with sides 2a and 
2t, where a < 3, as indicated in Fig. 2. Let the cir- 
cumferential parameter s be taken as zero at the point 
y =0,2=b. Then for0 < s <a,r = 5, and fora < 
s < 2b+ a, = ra. Because of symmetry considera- 
tions, the warping u vanishes ats = 0,s =a+b,s = 
2(a + 6), ands = 3(a+ 5). Therefore, it is sufficient 
to consider only one corner. At the corner s = @, or 
y = —a,2 = b, Ar = a — b. Furthermore J = 
4(a + b)abt. Hence, Eq. (16) reads 


u\ | Ma@a—b)  (@-—bd) 
45;)- 


4G(a + b)abt (a + d)ab 
where U is the value of u at s = a. 
As a special example, assume that at the end section 
= 0, warping is completely suppressed, or u = 0, 


U (18) 
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Fic. 2. Rectangular section. 


while at x = ©, the warping is not restrained, or ¢,; = 
0. 

Hence, by Eq. (4), ¢- = Ou/Ox = Oatx = ~. Con- 
sequently, according to Eq. (13), 0°u/0s? = 0, i.e., u 
is a linear function of s. Then if U; is the value of u 
ats = aforx = o, the condition that u vanishes at the 
centers of four edges requires that 











Ou Oni a+ob_ 
ae: Fag git egg 
From Eq. (18) it follows that 
atb _ M(a — b) _. o~ oF — 
saith as 4Gabt(a + b) ab(a + 5) Ue Ge) 
or 
U, = U(@) = (M/16Gt)/[(a — b)/ab] (20b) 


Hence, U; at s = a is negative fora < 3. 
Then the warping function at the far end x = 
may be written in terms of Uj as 


: 
m(s) = U(,s) = s — for —ags Ka 
(20c) 


u,(s) = u(~,s) = —s* for —bgs* <¥b 


where s* = s — (a+ 5). And, also, similar expressions 
for the other two edges. Finally, from Eq. (14), 


(©) = M/16Gt [(a + b)/a*b?] (21) 


Eqs. (20b) and (21) are well known from the theory of 
uniform twist with no axial thrust. The solution as 
specified satisfies the free end condition at infinity. In 
order to meet the condition at the restrained end, one 
can write “4 = u, + we, where wu is a function of x and 


s to be determined by substituting u into the differential 








Eq. (13). Thus, 
2. 2 
a ot = 0, -a Ss Ka,8 30 
O7u2* O7u2* (22) 
dE? Os*2 -. —b g a” < b, é 20 


such that uw, u2* are odd in s and s*, respectively, and 
satisfying the following end conditions and transition 


conditions 





oe = or = ate o (23a) 

u2(0, s) = = u,(—-a<s<a)  (28b) 

uit(0, s*) =< Uy (BS s* <b) (Be) 

ue(t,a) = u2*(t, —b), & 2 0 (24a) 

(Sr) id (=). = —Im(t,a),§>0 (24b) 


Here U; is given as in Eq. (20b) and I is the abbrevia- 
tion 
T = (6 — a)*/ab(a + 5) (25) 


The solution of restrained torsion of a rectangular box 
column is then given by 


u(é, s) = -U; + um(t,s), -a Qs Ka 


* (26) 
u*(é, s*) = — put u2*(E, s*), —b S s* Fb 
Consider solutions of the type 
u2(&, s) = Ae~™ sin As (27) 
uo*(t, s*) = A*e~™ sin ds* o 


where \ > 0, the conditions in Eq. (23a) are satisfied. 
They are odd functions in s and s*, respectively. The 
homogeneous transition conditions (Eqs. 24a, b) give 
two homogeneous equations for the determination of 
the two unknowns A and A*. The condition that a 
solution different from zero exists is the relation 


Asin A(a + b) — [sin Xa sin Nb = O (28) 


Hence, one obtains the following set of characteristic 
values 


An = [x/(a + 5) ] yn, m = 1, 2, 3,... (29) 
where y, are given by 
FYn Si Ty, — I’ sin ry,8 sin ty,(1 — B) = 0 (30) 
The abbreviations are 
I’ = (a + 6)r = (1 — 28)*/8(1 — 8) = (31) 
B = a/(a + Dd) 432) 
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It should be noted that +, is a function of Bonly. The 
corresponding characteristic functions are 


= Csc d,@ sin d,5, o,* = —Cscr,b sin \,s* (33) 


It can be shown from the differential equations and the 
transition conditions that (¢,, ¢,*) form an orthogonal 
set. The orthogonality conditions are 


Sa" OnPm IS + So” On*Gm_* ds* = 
The normalized factors for m = 1, 2, ... are 
On = So on? ds + f} on*? ds* = 

a Csc* h,a + b Csc? Xb — 


0, form ~m (34) 


<3 (85) 
This orthogonal set is complete in the following sense: 
for any given pair of functions f(s) and f*(s*), being odd 
in—a ££ s Ka, —b < s* < |, respectively, one can al- 
ways expand them formally into the following pair of 
series 


f(s) = LAsex(6), — sga 
(36) 
f*(s*) = TAnea(s*), —bgs* gb 
n=1 
where the coefficients A, are determined by 


i, 5, Sifisels ds + S2f*(s*)\en*(s*) ds*] (37) 


An elaborate system of conditions and theorems can be 
worked out for this series in parallel with that of other 
familiar orthogonal series, for example, Fourier series, 
However, for the present simple purpose this discussion 
is unnecessary and will not be pursued. 

Now the Saint-Venant solution described can be ex- 
panded into this series at once. The coefficients of the 
series are given by 


A, a 8U,/(a + b)An*Qn (38a) 
or 
AY ts 8U, 
ah » (a + b)An?On * 
Csc \,@ sin d,s, —a fs fa 
‘ ‘ ant +. (38b) 
s 1 
a ae ee Del es 
‘ b 2 (a + b)dn?Qn * 





Csc d,b sin ys*, —b < s* f bj 


Therefore, the solution for the restrained warping func- 
tion for the present case is 
8, < 1 
u(é,s) = — U, — —————- x 
) (a + b) n=T Qn Xn? ‘" 
Csc Aya sin d,se~™™, (—a < s < a) 


Ss «. “1 
7+ eae FON 
Csc “yi sin \,s*e~*™, (—b < s* < b) |} 


(39) 





xX 


It is recalled that Uj, \,, Q, are given, respectively, by 
Eqs. (20b), (29) and (35). From Eg. (29), it can be 
shown that for large values of n, \,, approaches asymptot- 
ically to rn/(a + 6). Hence, the infinite series in Eq. 
(39) converges as n~* for — = 0, and, in fact, u(0, s), 
u*(0, s*) vanish by Eq. (38). 

From Eq. (14) one finds at once the rate of twist. 

86 - a)? 1a, 

we) = (2) 1 — PEEPS em] co) 
where 6() is the rate of uniform twist in the Saint- 
Venant theory, as given by Eq. (21). By Eq. (38b), one 
has ats =a 


> 8/(a + B20, = 1 (41) 


Hence, at the end — = 


(6 — a)? 4ab B M 
eral eer a 
(42) 
It should be noted that this result is different from 
that given in the previous paper,® where the over- 
simplified approximations give vanishing @(x) at the 
restrained end x = 0. Present work agrees with the 
physical fact that the twist is the necessary consequence 
of application of moment. As a matter of fact, this 
result can be shown to be true in general by Eq. (9). 
At x = 0, u(x) = 0, then Eq. (9) immediately gives 


6(0) = M/GI (43) 


(0) = @(@) [1 - 


Eq. (43) holds for any cross section. 

For the variation of @(x) against x where kx = &, the 
characteristic values \, must be computed from Eq. 
(30). This is a transcendental equation. A series of 
roots have been computed numerically for 8 = 0.1, 
0.2, 0.3, 0.4. They are tabulated in Table 1. These 
characteristic values approach integers as m becomes 
large. The values of 6(x) are computed for various 
values of 8. They are plotted in Fig. 3. It should be 
noted that 6(x) approaches a constant as 8 approaches 
1/,. For B = '/s, one has a = 8, i.e., a square section; 
therefore, as shown in the previous section by Eq. 
(17), @(x) is constant throughout. 

The circumferential distribution of axial stress at the 
restrained end section o,(0, s) has also been calculated. 
These stresses are plotted in Fig. 4 in dimensionless 








TABLE 1 
The Characteristic Values 7, for Various Values of 8 
B 0.1 0.2 0.3 0.4 
I’ 7.111 2.250 0.762 0. 1667 
‘1 0.836 0.891 0.944 0.984 
2 1.778 1.886 1.965 1.997 
Ys 2.786 2.933 2.981 3.0 
v% 3.818 3.981 3.993 4.0 
Ys 4.859 5.000 5.0 5.0 
6 5.901 5.986 6.0 6.0 
v1 6.940 7.0 7.0 7.0 
Ys 8.0 8.0 8.0 8.0 
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Fic. 3. The rate of twist at various distances from the restrained 
end. 
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Fic. 4. Circumferential distribution of axial stress at the 
fixed end of a rectangular thin-walled section with 2a X 2b cross 
section, 8 = a/a + b. 


form. It is found that the axial stress is maximum at 
the four corners. This maximum stress is also plotted 
in Fig. 5 against various sizes of the rectangular thin- 
walled section. 

The foregoing treatment is based on the assumption 
that a continuous distribution of web system prevents 
any circumferential strain. It is evident, however, 
that a plane sheet itself is capable of transferring any 
plane stress system without the help of any internal 
_ web system; thus, the assumption of the web system 
mentioned above is not necessary for this special case. 
For such a cylindrical box column without the help of 
an internal web system, local reinforcement is required 
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at the junction of the plane sheets by stringers. These 
stringers must be capable of resisting the normal force 
resulting from the stresses in the plane sheets meeting 
at the edges. The theory of such a rectangular box sec- 
tion structure reinforced with edge stringers was worked 
out by D. Williams.’ Unfortunately, Williams’ work 
is extremely complicated, and no numerical results have 
been explicitly given for the case of a box section with 
uniform thickness. Therefore, no comparison can be 
made between the present case of a thin-walled box 
section with internal web system and the Williams’ 
case of the same section reinforced with edge string- 
ers. 


TORSION OF A THIN-WALLED COLUMN WITH 
CONTINUOUS ?(s) 


Consider a thin-walled column with continuous r(s), 
As a special example assume that at the end section at 


x = 0, warping is completely suppressed (i.e., a fixed , 


end), while atx = ~, o, = 0 (i.e., a free end). 

A particular solution of the Eq. (12) is well known as 
Saint-Venant’s solution for unrestrained torsion. Con- 
sider the case where u, = 1;(s), 0; = constant, and 7; = 
constant. Integration of Eq. (6) gives 


ns = G[m(s) + AS’r ds] (44) 


where the point s = 0 is so located that 1(0) is equal to 
zero. From the periodic property of u(s), i.e., w(/) = 
u,(0) = 0, where / is the total circumferential length of 
the section, the value of @, is determined from Eq. (44) 
by putting s = 1. This gives 


0; = 71/2GA (45) 


Here A is the area bounded by the centerline of the 
cross-sectional ring. 








A= /,$rds (46) 
Finally 7; is determined from Eq. (8) 
71 = M/2tA (47) 
4.0 | — 
dw 5; e 
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Fic. 5. Maximum axial stress at corner for various rectangular 
thin-walled sections. 
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Hence, one has at once from Eqs. (44), (45), and (46) 
by eliminating 7, 61, 


Mi 
ai 4GtA? 


It is easily verified that this is a particular solution of 
Eq. (12), well known as Saint-Venant’s solution of uni- 
form torsion. This solution evidently satisfies only the 
unrestrained end conditions at ends of the column. 

In order to meet the condition at the restrained end 
x = 0, asolution of a homogeneous equation correspond- 
ing to Eq. (13) must be found such that it vanishes at 
infinity and equals the negative of the Saint-Venant 
solution atx = 0. The solution of restrained torsion is 
then given by the sum of this solution and the §aint- 
Venant solution. In other words, the problem is now 
reduced to findings the solution of the homogeneous 
equation 


u(s) = Ser ds (48) 


2. 2. 
a we ~ is us ods =0 (49) 
such that 
We = —u atx = 0 (50) 
Ou2,/0x = Oatx = @ (51) 
Hence, the solution of restrained torsion is given by 
u(x, Ss) = m(s) + u(x, s) (52) 


The solutions of Eq. (49) satisfying the condition at x = 
© may be written as 


uz = Ae~™*U(s) (53) 
where U(s) satisfies the equation 
@U td 
MU + > +55 = 0 (54) 
Ids 


and the condition 
£,U(s)(dr/ds) ds = 1 (55) 


It should be noted that \ is the characteristic value of 
the problem to be determined. The periodic solution 
of Eq. (54) can be obtained at once: 
t cos A(s — 5 + 4/2) : 
U(s) = —- — ’(s1) 3 - 
(s) r’(s1) ain */gd ds, 


8 


r'(s:) sin Ms — s;) ds; (56) 





I 
where r’(s) denotes the first derivative of r(s) with re- 
spect to s. This solution must also satisfy the condi- 


tion in Eq. (55). Hence, one obtains the following 
equation for the determination of X: 


- — mr ot 1/21 £ Gf. ’(s)r’(s1) cos A(s — sy) dsyds— 


‘¢ Pi (s}) sim X(s — s,) dsyds = 1 (57) 


These characteristic values and the corresponding char- 
acteristic functions will be denoted by \, and U,(s). 
The orthogonality of these functions can be easily es- 
tablished in general by means of Eq. (54) and the condi- 
tion in Eq. (55). Furthermore, let NV, be the normaliz- 
ing factor defined by 


N, = $:U,(s)U<s) ds (58) 


Then the Saint-Venant solution u,(s) can be ex- 
panded in terms of this set of orthogonal functions. In 
other words, 


m(s) = SA UiAs) 69) 
where 
A,= (1/N;) Sf, uy(s)U,(s) ds (60) 


And consequently, the solution for the problem of re- 
strained torsion is, by the condition in Eq. (52), 


u(x, s) = m(s) — y Ae U5) (61) 


This method can be easily extended to finite columns 
with various end conditions. In such cases, the funda- 
mental characteristic functions will be more compli- 
cated than for the particular case of the infinite column 
given. 


TORSION OF A THIN-WALLED SECTION WITH DUMBBELL- 
SHAPED SECTION 


A simple example of the above method will be dis- 
cussed. Consider the case of a cross-sectional ring of 
the form 


r = % + 7 cos 4x(s/I) (62) 


which is a dumbbell section with major axis 7; + ro and 
minor axis 7 — 7; (Fig. 6). The Saint-Venant solution 
can be obtained from Eq. (48) as 


sin 4a — (63) 





u(s) = — 


From Eq. (58) one obtains 


4r ¢ sin 4x(s/l) 
_ if 
1 ‘I (4x/l)? — 2 








Uy) = - (64) 





lof 














Fic. 6. Dumbbell shape section. 





510 JOURNAL OF THE AERONAUTICAL SCIENCES—OCTOBER, 1946 


where \ is determined by the condition in Eq. (18). 
After simplication, this condition yields 


4r\? tlr;? 4r\? 7,2 Lo 
2 = _- Oo = _ 1 oe 
F)0-F)-Gete) © 
where 
I = It(ro? + 1/27”) (66) 
Consequently, by the condition in Eq. (50), one obtains 


Mn ’ 4r —kaz 
i 1 - 67 
SrGir,? sin (= s)¢ aoe ee 





u(x,s) = — 


This is the warping function of a semi-infinite column 
of cross section given by Eq. (62) with one end com- 
pletely prevented from warping. With u(x, s) known, 
the rate of twist @(x) and the longitudinal tension ¢, 
can be easily calculated. 


aa M ne om | rn 
a(x) = x3 haga 3 #0 (68) 








CO; = 


p Mp, . {4r — 
aes a Dar? hk sin (7 s)e (69) 


From the above analysis it is obvious that the end effect 
of the restrained condition extends to a depth directly 
proportional to the length of the circumferential arc. 
In fact, this conclusion can be generalized for thin- 
walled columns of any cross section. Let s = lo, r = 
7n, where 7 is the maximum value of r over the circum- 
ference, and 0 < o < 1,|n| < 1. Hence, 


l=tfrds = FS do (70) 


With these substitutions, Eq. (57) can be reduced to a 
dimensionless form. The solution of \ will be, there- 
fore, of the form 


hat = Ny, m =1,2,... (71) 


where N,, are discrete characteristic values. Hence, 
the warping function is the sum of the Saint-Venant 
solution and all the characteristic functions of the type 
given by Eq. (53), which vanish exponentially as x in- 
creases. The component that has the deepest pene- 
tration is the one corresponding to the lowest charac- 
teristic value, ie., Nj. This component vanishes ex- 


ponentially as exp. (—WN,kx/l). Hence, the above 


conclusion follows at once. 


CONCLUSIONS 


The simple theory presented in this paper gives a 
reliable method for computing the stress distribution in 
a thin-walled column of uniform thickness which is sub- 
ject to variable twist. The differential equation for 
the warping function and its solution have been dis- 
cussed in detail. In the case of a box section con- 
sisting of a number of flat sheets, the differential equa- 
tion of the warping function is reduced to the Laplace 
equation together with certain transition conditions at 
the edge points. The case of a rectangular box column 
has been worked out in detail as an example. 

The results of this theory deviate considerably from 
that given in the previous paper® where, as a result of 
the simple method of approximation used, the rate of 
twist is equal to zero at the restrained section. The 
present theory, however, yields in general the value 
M/GI for the rate of twist at the restrained section, 
which checks the physical picture of restrained tor- 
sion. 
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Nondimensional Buckling Curves—Their 
Development and Application 


F. P. COZZONE* anp M. A. MELCON? 
Lockheed Aircraft Corporation 


SUMMARY 


Basic expressions for column buckling, initial buckling in com- 
pression, initial buckling in shear, interrivet buckling and local 
crippling are reduced to a common parameter—namely, the 
ratio of tangent modulus to buckling stress (E;/F). 

A nondimensional buckling parameter B is plotted as a func- 
tion of the parameters F/ Fo. and m, where Fy.z is a secant yield 
stress and m is determined by the shape of the knee of the stress- 
strain curve. 

Constants for evaluating the buckling parameters for the 
various types of buckling phenomena, values of the secant 
yield stress Fo.;, and the shape parameter m are tabulated for the 
materials commonly used in aircraft. 

It is shown that these nondimensional buckling curves,.when 
entered with the proper buckling parameter, may be used to 
evaluate the buckling values required for design. 


INTRODUCTION 


Tl? PREVALENT USE of semimonocoque aircraft 
structures has brought with it the problem of 
determining the buckling characteristics of a large 
variety of stiffeners and of stiffened sheet combinations 
for various materials. The type of data found neces- 
sary to adequately determine these buckling character- 
istics are as follows: 

(1) Compressive stress-strain data based on mini- 
mum guaranteed values. 

(2) Tangent modulus versus compressive stress 
data. 

(3) Column curves based on various yield stresses. 

(4) Initial buckling in compression for flat and 
curved elements. 

(5) Initial buckling in shear for flat and curved ele- 
ments. 

(6) Interrivet buckling curves. 

(7) Local crippling curves for flat and curved ele- 
ments. 

Such data are primarily dependent on typical com- 
pressive stress-strain data. However, up to afew years 
ago, no compressive stress-strain data were available for 
thin sheets. In such cases, compressive stress-strain 
properties were assumed equal to the tensile stress- 
strain properties. 

With the improvements in the technique of obtaining 
compressive stress-strain data on thin sheets, much of 
the buckling data became available for the relatively 


Presented at the Structures Session, Fourteenth Annual 
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few materials then in use. Then came the use of 5 
per cent Alclad, postaging, and the new high-strength 
aluminum alloys, together with an upward revision 
in the minimum guaranteed values of most of the 
former aluminum alloys and the adoption by the 
Army and the Navy of design values higher than mini- 
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FIGURE § 
BUCKLING RELATIONSHIPS 
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NONDIMENSIONAL BUCKLING CURVES 


mum guaranteed. Where formerly the aircraft struc- 
tures engineer and designer had available but a few sets 
of physical properties, he now has available a total of 
63 sets for the aluminum alloys alone. With 16 sets 
or the stainless steels and five for the chrome- 
molybdenum steels, the total is 84 different sets of 
physical properties. 

This imposed a tremendous task on the research de- 
partments of the various manufacturers since they had 
to furnish the data necessary for design and analysis. 

Fortunately, at this time Osgood proposed his equa- 
tion for the stress-strain curve. Ramberg expressed 
the two arbitrary parameters in this equation in terms 
of two points of the stress-strain curve. The equations 
developed are of such a nature that the tangent modulus 
relationships and the buckling relationships of all air- 
craft materials can be expressed and plotted in a non- 
dimensional form. [The balance of this paper will deal 
with the application of the Ramberg-Osgood expression 
to the development of the ‘‘buckling allowables’’ for 
aircraft materials. ] 


GENERAL 


Basic Nondimensional Relationship 


The basic Ramberg-Osgood relationship may be ex- 


pressed as 
e-“+ieG) | 
aw (ase DAL) ve wedi 
 tanlien “lee ie,” 


(plotted as Fig. 8 of reference 2), where 


E, = the tangent modulus of elasticity 

E Young’s modulus of elasticity 

n = the stress-strain shape parameter expressed 

as 1 + {0.3851/[logio (Fo.7/Fo.ss)]} (plotted as 

Fig. 10 of reference 2) 

F = the allowable stress 

Fo.7 = the secant yield stress corresponding to the 
stress at the intersection between the 
stress-strain curve and a secant line through 
the origin having a slope of 0.7E 

Fo.35 = the secant yield stress corresponding to the 
stress at the intersection between the 
stress-strain curve and a secant line through 
the origin having a slope of 0.85£ 


In the case of Alclad materials, the initial or Young’s 
modulus is not constant throughout the so-called elastic 
range because of the yielding of the Alclad coating at a 
telatively low stress. This gives a primary modulus 
and a secondary modulus. Since the use of Alclad in 
aircraft structures is widespread, some means was 
necessary to utilize the basic Ramberg-Osgood rela- 
tionship. The following approximation is recom- 
mended for the initial modulus EZ in Alclad materials: 

10,000 | Fy — oe) 


pee / 
. Few ( E, : E, 
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where 


F.y = compressive yield stress as determined by the 
stress corresponding to the intersection be- 
tween the stress-strain curve and a line 
parallel to the initial modulus at an offset 
strain of 0.002 
E, = primary modulus of elasticity with a value of 
10.5 X 10° 

secondary modulus of elasticity with a value of 
9.5 X 10-* for 10 per cent Alclad; and 
10.0 X 10* for 5 per cent Alclad 


E; 


This approximation has been found adequate, since 
it does not have an appreciable effect on the buckling 
allowables. 


Buckling Relationships 


The basic buckling phenomena are usually expressed 
in the following form: 


Column, F, = 2? E,/(L’/p)?. 

Initial compression buckling, Fc, = KE, (t/b)?. 

Initial shear buckling, Fs,, = K,E,, (t/b)*, where 
Es, = shear tangent modulus. 

Interrivet buckling, F;, = 4? E,/(s/0.29t)? where 
s/0.29t is comparable to L/p of an element with L = s, 
the rivet spacing. 


All of the above can readily be put into the form 





AY 2 
E,_(L'/p)? _ (b/t)? _ (b/t)? _ (au) 
fF & ae 
|. “Interrivet buckling 
Initial shear buckling 


L Initial compression buckling 
Column 


Therefore, the problem resolves itself into obtaining 
an expression for E,/F from the basic nondimensional 
relationship. 

To do this, multiply both sides by Fo.7/F and equate 


to B?, 
wah 
FRA2/) .8..3 (F?*¥ 
tie) 
Fo.z 


For 7 


— 








Fig. 9 is a plot of this expression and shows F/F.7 
versus B for various values of n. 

The shape of the knee of the stress-strain curve is 
given by the shape parameter m, and the abscissa B 
incorporates the particular properties of the material 
Fo.7/ E. 

Insert in turn each of the expressions, E,/F, for the 
various buckling phenomena in the basic expression for 
B? in order to determine the buckling constants for 
evaluating the buckling parameter B. 





Boasrt 





bst 


Column Buckling Constant 
B?, = (Fo.1/E) (1/x*) (L’/p)? 


WFO)-+€) 
B, ws 9 —E re = k, > 


where the column buckling constant, k, = (1/7) X 
V Fox/E. 


Initial Compression Buckling Constant 


er = (E)(E)G) 
Boe = V8 (Yee) ~ VE (i), ~ * (i), 


where the initial compression buckling constant, 


keep = WFox/E and (b/t), = (b/t)/WKe 


Initial Shear Buckling Constant 


Nadai® states that, according to tests by American 
engineers, it appears that the yield point for soft metals 
in pure shear possesses a value of 0.6F,,. He presents 
an analytical approach that shows this to be 0.577 Fy. 

Since for many materials used in aircraft Fy, is ap- 
proximately equal to 0.61F,,, it was felt that the shear 
stress-strain curve could be represented by .a curve 
having the same shape factor as the compressive 
stress-strain curve but having values of F,, and (Fo.7), 
equal to 0.61F,, and 0.61(Fo.7)-, respectively. 
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Therefore, the corresponding value of B for initial 
shear buckling may be obtained from the initial com- 
pression buckling by merely substituting 0.61F .7 for 
Fo.7. 

This gives 


Bs a 0.61 Fo.7 (7) ae ks (*) 
cr ‘ E t . cr t a 


where the initial shear buckling constant, ks,, = 


V0.61Fo.7/E and (b/t), = (b/t)/WK,. 





Interrivet Buckling Constant 
(By)? = (Fo.x/E) (1/4x*) (s/0.29t)? 
By, = (1/22V0.29) V Fox(s/t)/E = kys/t 


where the interrivet buckling constant, ky, = 


(1/2"V 0.29) V Fo2/E = 0.551V Fo.x/E. 


Local Crippling or General Buckling Constant 


In the case of buckling phenomena where empirical 
buckling curves are derived from testdataon a particular 
material, conversion to any other material may be 
made as follows: 

Consider the case of crippling curves such as are 
used at Lockheed. These curves are empirical curves 
obtained from test data on 24SRT-10 per cent Alclad. 
Such curves are shown in Fig. 6. 

Consider one of these curves as shown in Fig. 11 
where the ordinate scale has been converted from F,, to 
F.¢/ Fo.1. 

For any value of b/t a corresponding value of F,,/Fo.1 
may be obtained, and for this value of F,,/ Fo.7 a corres- 
ponding value of B may be obtained as shown in Fig. 9. 
This value of B may be designated as Basspr. 

Such values of Boyspz have been found in this manner 
from the original empirical crippling curves and plot- 
ted in Fig. 12. 

Now to convert from 24SRT material to some other 
material, use can be made of the basic buckling relation- 
ship (E,/F..)1 = (E,/Fec)oasrt and B? = (E,/F)(Fo.7/E). 
Hence, 


[B?(E/Fo.7) hi, = [B? (E/Fo.7) lessrt 
and 


(Fo.7/E): 
(Fo.7/E) asset 


Bee = Ree Bossrt 


Bi = Bossrt 


where the local crippling constant, k,, = V (Fox/E): + 
V Foa/E)asne- 


Physical Properties 


For the various materials used in aircraft structure, 
a new set of physical properties has been determined 
ranging from shape parameter m to various buckling 
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constants k. These values have been obtained: from 
the average of available compressive stress-strain data 
for the materials in question. Since sheet data can 
only be considered typical, it was necessary to reduce to 
minimum guaranteed values. This was accomplished 
by means of an affine transformation. The mechanics 
of this transformation can best be described with the 
help of Fig. 14. 

For any typical curve, construct a line parallel to the 
initial modulus line through the 0.002 strain point. 
The intersection of this line with the typical stress-strain 
curve is the yield point F,, of the typical curve. The 
point A’ on this line is the corresponding minimum 
guaranteed yield point F,,. A line drawn through the 
origin 0 and A’ will intersect the typical stress-strain 
curve at point A. For any other points on the typical 
curve such as B and C, the corresponding points on the 
minimum guaranteed curve may be found as follows. 
Draw radial lines OB and OC. On line OB lay off OB’ 
= OB (OA’/OA) and on line OC lay off OC’ = OC(OA’/ 
OA). The curve through the points C’, B’, and A’ is 
the desired minimum guaranteed curve. 

These physical properties, together with the buckling 
constants for all the commonly used materials, have 
been summarized in Table 1. 


Typical Applications 


(1) Typical stress-strain curve is obtained from 
manufacturer of material or research laboratory. 

(2) Affine transformation gives stress-strain curve 
based on minimum guaranteed properties. 

(3) Compute ” from equation 


n=1i+ 0.3851/logio (Fo.7/ Fo.ss) 


(4) To obtain column curve of F, vs. L’/p 

(a) For each value of L’/p compute B = k,(L’/p). 

(b) Enter nondimensional column curve with 
these values of B to obtain for the proper u-value the 
corresponding values of F/F.7. 

(c) Multiply these values of F/Fo.7 by Fo.7 to 
obtain F, corresponding to originally assumed values 
of L’/p. 

(5) To obtain tangent modulus curves use Fig. 8 
of reference 2, which gives a plot of F,/E vs. F./Fo.1. 

(6) To obtain critical (initial) compression buckling 
curves for Fc, vs. (b/t),. 

(a) For assumed values of (b/t), compute B = 
Reer(b/t)e. 

(b) Enter nondimensional column curve with 
these values of B to obtain for proper m-curve the 
corresponding values of F/ Fo.7. 

(c) Multiply these values of F/Fo.. by Fo.7 to 
obtain F:,, corresponding to originally assumed values 
of (b/t),. 

(7) To obtain critical shear buckling curves of 
Fs,, vs. (b/t),. 
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(a) For assumed values of (5/t), compute B = 
Rser (b/ t)e. 

(b) Enter nondimensional column curve with 
these values of B to obtain for the proper m-curve the 
corresponding values of F/Fo.7. 

(c) Multiply these values of F/Fo.7 by 0.61F.7 to 
obtain F:,, corresponding to originally assumed values 
of (b/t),. 

(8) To obtain interrivet buckling curves of F;, vs. 


s/t. 


(a) For assumed values of s/t, the ratio of rivet 
spacing to skin thickness, compute B = k,,(s/t). 

(b) Enter nondimensional column curve with 
these values of B to obtain for the proper n-curve the 
corresponding values of F/ Fo.7. 

(c) Multiply these values of F/Fo.7 by Fo.7 to ob- 
tain F,, corresponding to originally assumed values of 


s/t. 


(9) To obtain crippling curves of F,, versus b/t and 


R/t. 


(a) Enter Fig. 12 with assumed values of b/t 
and R/t to obtain from proper 24SRT crippling 
curves the corresponding Bospr values. 

(b) Compute corresponding B values for each of 
Bossrt Values where B = k,, Bogspt. 

(c) Enter nondimensional column curve with 
these values of B to obtain for the proper m-curve the 
corresponding values of F/ Fo.7. 

(d) Multiply these values of F/Fo.. by Fo.7 to 
obtain F,, corresponding to originally assumed }/t 
and R/t values. 

(10) To reduce compression test results to an al- 


lowable based on minimum guaranteed properties. 
(This application taken from reference 6.) 


(a) Take a plot of stress-strain curve based on 
minimum guaranteed properties. 

(b) Draw the 0.002 yield offset line parallel to 
initial modulus and extend this line to the yield 
stress of actual material tested. 

(c) Draw a secant line through origin and actual 
yield stress of test material. 

(d) Draw an 0.7E secant line and by affine trans- 
formation determine F).; value for actual test mate- 
rial. 
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TABLE 1 
_ 24ST—Sheet, Plate, 24ST Sheet, Plate 24SRT Sheet, Plate 
Strip—Heat Treated by User As Received As Received 
Non- | Non- 
Alclad Alclad Alclad Non-Alclad Alclad Alclad 
0.012 to 0.064 All 0.012 to/0.064 to| 0.189 | 0.012to | 0.189 to | 0.012 to 0.064 All 
0.063 up gages 0.063 0.188 up 0.188 1.00 0.063 up gages 
n 15 15 15 10 10 10 10 10 10 10 10 
Fo.t 36,000 37,000 39,000 38,000 | 39,000 | 37,000 41,000 39,000 46,000 48;000 50,000 
E x 10-6 9.76 10.1 10.5 9.75 10.1 10.1 10.5 10.5 9.70 10.1 10.5 
Fey 7 ,000 38,000 40,000 39,000 40,000 38,000 42,000 40,000 46,000 48,000 50,000 
Column, ke — 0.0193 0.0193 0.0194 0.0199 | 0.0198 | 0.0193 0.0199 0.0194 0.0219 0.0219 0.0220 
Initial Buckling (Comp.), ke,, .060' 0.0605 0.0609 0.0624 | 0.0621 | 0.0605 0.0621 0.0609 0.0689 0.0689 0.0690 
Initial Buckling (Shear), ks, 0.0474 0.0473 0.0476 0.0488 | 0.0485 | 0.0473 0.0488 0.0476 0.0538 0.0538 0.0539 
Interrivet Buckling, kir 0.0335 0.0333 0.0336 0.0344 | 0.0342 | 0.0333 0.0344 0.0336 0.0379 0.0380 0.0380 
Crippling, kee 0.882 0.879 885 0.907 0.902 0.879 0.907 885 1.000 1.001 1.002 
24ST Extrusion 24STS80 Sheet, Plate, Strip 
Heat Treated 
by user As received Alclad Non-Alclad 
All Up to 0.250 to | 0.750 to 1.50 0.012 to | 0.032 to 0.064 0.012 to 0.125 
gages 0.250 0.749 1.499 up 0.031 0.063 up 0.124 up 
n 10 10 10 10 1 Sagoed 25 20 Jae 
For 37,000 37,000 40,000 44,000 50,000 cael 47,000 49,000 cies ae 
Ex 10-7 10.5 10.5 10.5 10.5 10.5 | 9.73 10.1 | 10.5 
cy ,000 38,000 41,000 44,000 : 47,000 49,000 50,000 
Column, ke 0.0189 0.0189 0.0196 0.0206 0.0220 Not 0.0221 0.0222 Not gave 
Initial Buckling (Comp.), ke,, 0.0594 0.0594 0.0617 0.0647 0.0690 to be 0.0695 0.0697 to be 
Initial Buckling (Shear), ks,, 0.0464 0.0464 0.0482 0.0506 0.0539 Aged 0.0543 0.0544 Aged 
Interrivet Buckling, kir 0.0327 0.0327 0.0340 0.0357 0.0380 | 0.0383 0.0384 | et 
Crippling, Ree 0.862 0.862 0.896 0.940 1.002 1.009 1.011 os 
24ST81—Sheet, Plate 24ST86—Sheet, Plate 
Alclad Non-Alclad Alclad Non-Alclad 
0.012 to | 0.032 to | 0.064 to 0.012 to 0.125 to 0.012 to | 0.032 to 0.064 0.012 to 0.125 
.03 0.0 0.188 0.124 0.188 0.031 0.063 up 0.124 up 
” ‘ 25 20 gids nese 25 20 sabe 
Fo.1 | 55,000 58,000 | Be es | 63,000 67,000 | plnaee 
E % 10-8 onaae RR onat an eats Pe 
r 57,000 ,000 ,000 ,000 ’ 
Column, ke _ st. | 0.0240 | 0.6241 Fag a ph Net. | 0.6257 | 0.0259 aa: ake 
Initial Buckling (Comp.), ke,,| Aged 0.0754 0.0758 Aged Aged 0.0808 0.0814 Aged 
Initial Buckling (Shear), ke, f 0.0589 0.0592 | | 0.0631 0.0636 | 
Interrivet Buckling, kir 0.0415 0.0418 0.0445 0.0449 
Crippling, kee 1.095 1.100 1.173 1.183 
24ST8X Extrusion 14ST Extrusion | 75ST—Sheet, Strip 75ST Extrusions 
Heat Treated, 
aged from Aged as 
from 24SO received As received Clad 
All Up to 0.125 to 0.500 to 0.750 To 0.040 to All 
gages 0.749 0.49 0.749 up 0.039 0.249 gages 
n 20 20 35 35 35 10 10 20 
Foz 50,000 52,000 50,000 55,000 58,000 62,000 63,000 71,000 
E x 107% 10.5 10.5 10.5 5 10.5 9.67 10.1 10.5 
Fey 50,000 52,000 ,000 55,000 58,000 60,000 62,000 70,000 
Column, ke — 0.0220 0.0224 0.0220 0.0230 0.0237 0.0255 0.0251 0.0262 
Initial Buckling (Comp.), ke,, 0.0690 0.0704 0.0690 0.0724 0.0743 0.0801 0.0790 0.0822 
Initial Buckling (Shear), ke,, 0.0539 0.0550 0.0539 0.0565 0.0580 0.0625 0.0617 0.0642 
Interrivet Buckling, hir 0.0380 0.0388 0.0380 0.0399 0.0410 0.0441 0.0435 0.0453 
Crippling, kee 1.002 1.022 1.002 1.051 1.079 1.163 1.147 1.194 
CrNi (Stainless) Steel 
With Grain 
1/,H 1/2 H 8/,H FH 
1/,4H 1/,H 3/,H FH B.T. ET. =F. a. 
n 3 3 3 3 6 
Fo.t 49,000* 58,000* 66,000* 80,000* 61,000* 101,000* 119,000 142,000* 
E xX 10-¢ 26 26 26 26 27 2 2 
cy 65,000* 71,000* 83,000* 91,000* 66,000* 104,000* 120,000 139,000* 
Column, ke 0.0138 0.0150 0.0160 0.0177 0.0151 0.0195 0.0211 0.0231 
Initial Buckling (Comp.), ke,, 0.0434 0.0472 0.0504 0.0555 0.0475 0.0612 0.0664 0.0725 
Initial Buckling (Shear), ks,, 0.0339 0.0369 0.0394 0.0433 0.0371 0.0478 0.0519 0.0566 
Interrivet Buckling, kir 0.0239 0.0260 0.0278 0.0306 0.0262 0.0337 0.0366 0.0400 
Crippling, kee 0.630 0.686 0.732 0.806 0.690 0.888 0.964 1.053 
CrNi (Stainless) Steel 4130 and 8630 Steel 
Cross Grain 
1/4H 1/2H 3/4H FH Nor- {100,000 | 125,000 | 150,000 | 180,000 
1/,H 1/27H 3/,H FH H.T. HT. H.T. H.T. mal SoA H.T. s H.T. 
n 4 5 6 8 6 7 = 20 25 35 40 e 50 
Fo.1 93,000 | 118,000 | 142,000 | 177,000 | 65,000* | 124,000 | 165,000 | 183,000* | 74,000 | 94,000 | 114,000 140,000 | 165,000 
E x 10+ 27 27 27 27 28 28 28 28 29 29 29 29 29 
F 100,000 | 120,000 | 140,000 | 170,000 | 75,000* | 125,000 | 160,000 | 177,000* | 75,000 | 80,000 | 100,000 135,000 | 165,000 
Column, ke 0.0187 | 0.0210 | 0.0231 | 0.0258|0.0153 | 0.0212] 0.0244] 0.0257 |0.0161/0.0181 | 0.0200; 0.0121 0.0240 
Initial Buckling (Comp.), ke,, | 0.0587 | 0.0661 0.0725 | 0.0810 | 0.0482 0.0665 | 0.0768| 0.0808 |0.0505|0.0569 | 0.0627 | 0.0695 | 0.0754 
Initial Buckling (Shear), ke, | 0.0458 | 0.0516 | 0.0566 0.0632 | 0.0376 0.0520 | 0.0600 |} 0.0631 |0.0395|0.0445| 0.0490 | 06.0543 0.0589 
Interrivet Buckling, kir 0.0323 | 0.0364} 0.0400] 0.0446 | 0.0265 0.0367 | 0.0423} 0.0445 |0.0278|0.0314| 0.0345 | 0.0383 0.0416 
Crippling, kee 0.852 0.960 1.053 1.176 |0.700 0.966 1.115 1.174 0.734 |0.827 0.910 1.009 1.095 





























* Values obtained from stress-strain curves which did not come up to minimum guaranteed values as specified in ANC-5. 
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TABLE 2 
Procedure and the Source for each Step in the Determination of the Various Buckling Allowable Stresses 
F. Fe F,.. Fir Fre 
Initial buckling— Initial buckling Interrivet eae 
compression 
Column ~ (flat plates) ———— shear buckling-—_—~ —Crippling———~ 

Item Source Item Source Item Source Item Source Item Source 

ke Table 1 Reer Table 1 Recr Table 1 Rir Table 1 Ree Table 1 

n Table 1 n Table 1 n Table 1 n Table 1 n Table 1 
Fo.1 Table 1 Fo.r Table 1 Fo.1 Table 1 Fo.7 Table 1 For Table 1 
A OPER (7) b/t (7) bit s _ spacing PS aes 

p t}e VK. t/e VK, t thickness ss 
pruas Pi tlie ektaty. b sig rw tert ots een Ore ee hans 0 we Basrt Fig. 12 
B ke =) B ke?) B b(?) B &(*) B RecBassrt 

F F F . F . F ‘ 

= ig. rig. a .9 = Fig. 9 os Fig. 9 

Fo.r Fig..0 Fo.7 Fig. 9 Fo.1 Fig Fo.7 ed Fo.2 % 

F - 7) ( ) (+) 

F. Peal 5) Feer Fo. z.) Freer 0.61.Fox( Foz F,, Fo.7 i; Fre Fo.7 1 





(e) On plot of F, vs. L’/p based on minimum 
guaranteed properties draw new column curve based 
on Fo.7 of actual test material, the shape factor of 
minimum guaranteed material, and k, = (1/r) X 
V Fo.1/E (actual test material). 

(f) Fig. 3 indicates such a plot for test specimens 
having various yields. 

(g) Enter this figure with actual test compression 
stress value at rupture and obtain an effective L’/p 
value from actual column curve. 

(h) For this (L’/p) value, read corresponding 
allowable compressive stress value from column curve 
based on minimum guaranteed properties. 





Table 2 summarizes the procedure involved in ap- 
plying the nondimensional buckling curves and the 
physical properties given in Table 1 to the problem 
of finding the various buckling allowables commonly 
needed by the aircraft structures engineer and designer. 
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SUMMARY 


The increased use of heavy wing skins, resulting from various 
aero-structural considerations at high speeds, has necessitated 
more exact plate instability calculations. Since heavy skins con- 
stitute a large percentage of the wing bending material, accurate 
determination of critical stress and effective width are necessary 
from a weight standpoint. 

Theory and tests show that the critical compressive stress of a 
stringer panel is intermediate between that for simple and 
clamped edge support; some function of the elastic restraint pro- 
vided by the stringer element along the unloaded edges of the skin 
panel. Theoretical methods for determining the degree of re- 
straint have not been well developed. Consequently, critical 
compressive stresses were determined on a large number of typical 
stringer panels. Analysis of the results indicated that as the 
critical stress increased, a transition occurred from clamped to 
simple edge support. 

Since the effective width problem initiates at the critical stress 
of the panel, the effective width equations were modified to ac- 
count for the empirically determined critical stress equations. 
This method yields results that are consistently in good agreement 
with experimental data. 


SYMBOLS 


plate length, in. 

plate width, in. 

plate thickness, in. 

plate instability coefficient 

w*k/(12(1 — v?)] 

critical compression stress, Ibs. per sq.in. 

strain, in. per in. 

Poisson’s ratio, vy = 0.3 for aluminum alloys 

modulus of elasticity, Ibs. per sq.in. 

secant modulus, Ibs. per sq.in. 

restraint coefficient, Yy = 45)b/D (reference 3) 
stiffness per unit length of stringer at unloaded edge 
flexural rigidity per unit length of plate, Z##/[12(1 — v?)] 
sheet edge stress, Ibs. per sq.in. 

sheet edge strain, in. per in. 

effective width at each side of rivet line, in. 

effective area at each side of rivet line, in.? 


arr os 


SUhtpm 


A, 


INTRODUCTION 


FFICIENT DISTRIBUTION OF BENDING MATERIAL, flut- 

ter considerations, and the importance of delaying 
buckling to maintain smooth aerodynamic contours at 
high speeds have necessitated the use of heavy skins on 
wing surfaces. The weight penalty that may be in- 
curred by inaccurate instability analysis of such a 
wing structure has required a reconsideration of the 
instability action of the various component structural 
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elements under compression loads. Determination of 
critical compression stress and effective width of 
panels that are stiffened by stringer elements fastened 
along the unloaded edges is discussed in the present 
paper. 

In panels that are simply supported at the loaded 
edges and stiffened by stringers along the unloaded 
edges, the torsional rigidity of the stringer acts as an 
elastic restraining medium against rotation. Conse- 
quently, it is conservative to assume that the panels are 
simply supported at the unloaded edges. This fact 
was recognized early in the period of theoretical develop- 
ment of plate instability since it is associated with 
boundary conditions for the differential equation of a 
plate under compressive stress. 

Timoshenko! cites references that ascertain the in- 
fluence of the stiffening element torsional rigidity upon 
critical compression stress. This work was expanded 
by Dunn? and Lundquist and Stowell:*:* these refer- 
ences present curves of the panel instability coefficient, 
k, as a function of the torsional] restraint provided by the 
stiffening element. Furthermore, Lundquist and Sto- 
well‘ indicate several deficiencies in the theoretical 
determination of the stringer torsional rigidity and 
present a method for ascertaining this property. 

Unfortunately, the assumptions made in reference 4 
in deriving the theory of stringer torsional rigidity are 
questionable when applied to the present analysis. 
Typical wing construction in current use consists of 
panels of relatively high bending rigidity stiffened by 
formed Z- or hat-section stringers. Consequently, the 
stringers cannot be assumed to be “sturdy” as required 
by theory. In addition, corrections to the assumption 
of perfect elasticity which are presented as part of the 
theory are subject to debate. 


CRITICAL COMPRESSION STRESSES OF TEST PANELS 


As a‘result of the limitations of the theory, it was 
hoped to evaluate the torsional rigidity of stiffening 
elements indirectly by determining buckling stresses of 
typical flat sheet-stringer panels. 

Fischel> presented buckling data for Alclad 24S-T 
test panels and noted a transition from a theoretical 
fixed-edge support at low critical compression stresses 
to a simply supported edge condition at higher stresses. 
Experimental data on a large series of long sheet- 
stringer panels tested at the N.A.C.A. and a smaller 
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TABLE 1 


Test Panel Construction 
Long Panels a/b > 5 








————Stringer Element -——— — Skin — Materials—-———— 
Group Source Type Thickness Range b/t Thickness Range Stringer Skin 
1 rm -A.C.A. Hes 
-R.R. L4F : 
01 0.064 to 0.187 in. 35to220 0.040 to 0.187 in. Alc. 24S-T84 Alc. 24S-T81 


(reference 6) 











ocr, = 64,000 Ibs. 


per sq.in. min. 


Ger, = 57,000 Ibs. 
per sq.in. min. 








2 NACA. Find 
A.R.R. -3" _— 
R.R. L4F01 rs 0.040 to 0.187 in.~ 35 to 300 0.032 to 0.187 in. Alc. 24S-T84 Alc, 24S-T81 
(reference 6) | ~ 
a Cer. = 64,000 Ibs. oer, = 57,000 Ibs. 
. Pal per sq.in. min. per sq.in. min. 
3 Republic Aviation | “ | 0.081 in. 45to70 0.081 to 0.125 in. Alc. 75S-T Alc, 75S-T 
tests 2 Cer. = 70,000 Ibs. oer, = 67,500 Ibs. 
‘re per sq.in. per sq.in. 
net ATS { 4 0.040 to 0.064 in. 75to 100 0.032 to 0.040 in. Alc. 24S-T oer, = 40,000 Ibs. 
4 Republic Aviation f Alc. 75S-T per sq.in. 
tests [ R301-T Ser, = 68,000 Ibs. 
os Fined per sq.in. 
& 5 Cer, = 58,000 Ibs. 
1 > per sq.in. 
5 Republic Aviation ‘e | 0.102 in. 125 0.040 in. 75S-T extrusion Alc. 75S-T 
tests ee Ss Ger, = 75,000 Ibs. er, = 66,000 Ibs. 
a _| per sq.in. per sq.in. 
ed 
Series tested at Republic Aviation Corporation verified cr, = 3.101 E(t/b)'-* (1b) 


Fischel’s observation. 


ANALYSIS OF TEST DATA 


Construction of the test panels was typical of current 
design. The important panel parameters are presented 
in Table 1. In all experimental work, panel buckling 
stresses were determined by use of SR-4 resistance-type 
electric strain gages mounted back-to-back at the panel 
center. 

The experimental critical stress data were separated 
according to the geometry of the stiffening elements. 
All J and Z stiffened panels were grouped together; all 
hat-section data were treated as a second group. The 
division of data is justified by the consideration that hat 
sections form closed cross sections when joined -to the 
panel and, consequently, possess considerably greater 
torsional rigidity than either the J or Z elements that 
are of open cross section. The trend of the test data 
definitely confirmed this grouping. . 

Evaluation of the experimental data by the method 
of least squares yielded the following empirical equations 
for flat long plates (a/b > 5); for Z* stiffened panels: 


Cer, = 2.124E(t/b)»-™ (1a) 


and for hat stiffened panels: 


* “Z” is used in a general sense—it refers to Z, J, and lipped 
Zand J stringers. 


In the usual form of the plate instability equation, 
the geometric term is squared rather than raised to the 
odd powers of Eqs. (1) 


Ger, = KE(t/b)? (2) 
and, hence, for the Z stiffened panels: 
K = 2.124(b/t)9-4 (3a) 
and for the hat stiffened panels: 
K = 3.101(b/t)®-148 (3b) 


Eqs. (1) are plotted in Fig. 1 and are shown for com- 
parison with the curves for simply supported and fixed- 
edge conditions. A definite transition from fixed-edge 
support at low panel critical compressive stresses to 
simple support at high critical stresses is indicated. 

This transition is empirically stated by Eqs. (3). 
In Fig. 2 this transition as a function of panel geometry, 
b/t, is clearly shown. 


Discussion of Analysis 


All critical stress measurements were performed on 
panels constructed of Alclad sheet. Consequently, be- 
low a stress of approximately 13,000 Ibs. per sq.in., the 
primary modulus of elasticity for aluminum alloys of 
10,500,000 Ibs. per sq.in. holds, while above this 
stress value the secondary modulus of 9,600,000 Ibs. 
per sq.in. is applicable. 








40 


HAT STIFFENED 
PANELS 


20 


° 
“ 
‘3, 
° 


SIPS 2 


Z STIFFENED 


PANELS 


on 
r) 





2 
20 30 . 40 60 80 100 150 §=6200 


% 


Fic. 1. Critical compressive stress of flat Alclad sheet panels. 
E = 9.7 X 106 lbs. per sq.in. a/b > 5. 


In determining the empirical buckling equations, 
Eqs. (1), the panel geometry and critical stress were 
measured. The secondary modulus of elasticity was 
assumed to be applicable throughout the entire stress 
range with the result that Eqs. (3), the K transition 
expressions, include the variation from primary to 
secondary moduli at low stress values. Therefore, in 
conjunction with the empirical critical stress equation, 
the secondary modulus of elasticity is to be used 
throughout. 

If the results presented thus far are employed on 
panels in which the sheet is not of Alclad material, an 
inherent conservatism should be noted. The low yield 
strength cladding is in the region of the most highly 
stressed material and, consequently, reduces the effec- 
tive thickness of the sheet somewhat. As the critical 
stress is a function of the square of the thickness, this 
may cause a significant change in the critical stress. 


Buckling Above the Proportional Limit 


In Fig. 1, critical compression stresses are depicted 
for a large range of b/t values. The largest measured 


buckling stress of the test panels was 50,000 Ibs. per 
sq.in., approximately the proportional limit for the high- 
strength aluminum alloys of which the panels were 
constructed. ; 

To use the results of this investigation with materials 
that possess a lower proportional limit, or to extend 
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these results to buckling stresses above the proportional 
limit of the test materials, the secant modulus method? 
for plate instability above the proportional limit may be 
applied. 

Above the proportional limit, the critical stress de- 
termined by elastic theory must be reduced by a factor 
n; i.e., 


as , 
Sor, = N58 cr, 


where ¢,,.’ is the critical stress calculated by elastic 
theory. The secant modulus method indicates that the 
n term is equal to the ratio of the secant modulus at the 
critical stress to the elastic modulus. 

Critical stress and critical strain are implicitly related 
by the secant modulus according to this theory and 
consequently, 


fer, = “SE =K(t/b)? (4) 
E,; 
where K is determined from Eqs. (3) or Fig. 2. The 
critical stress is the stress value corresponding to the 
critical strain on a stress-strain curve of the particular 
material. 


Buckling of Plates of Short Aspect Ratio 


The results presented thus far were determined from 
measurements on long panels, i.e., plates of aspect ratio 
greater than five. It is probable that these results can 
be applied to panels of shorter aspect ratio. 

Lundquist and Stowell® state that the plate instability 
coefficient K is a function of the restraint coefficient 


v = 4Sb/D (5) 

So is the stiffness of the stringer element per unit length, 
2 

4S. = "2? (oe, — oer) (6) 


where 


X = half-wave length of buckle, in. 
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Fic. 2. Instability coefficients for stiffened panels. 
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Fic. 3. Restraint coefficients for stiffened panels. 
I, = stringer polar moment of inertia about axis 


of rotation, in.* 
str. =. Stringer failing stress, Ibs. per sq.in. 
Cer. = Critical stress of sheet, Jbs. per sq.in. 
Letting B = 7°J,/bD and combining Eqs. (5) and 
(6a) 
(6b) 


(Gste. — Gc) 


eet 2 
(/b)? 
Since \/b = (a/b)/m, where m is some integer of 
half-wave lengths, and a/b is the plate aspect ratio, it 
is possible that the restraint coefficient, y, may be 
changed a significant amount for plates of short aspect 
ratio. 

Since changes in half-wave length cause small varia- 
tion in critical stress of plates with aspect ratios greater 
than 2.5, it is probably permissible for practical con- 
siderations to establish 2.5 as the lower limit of ap- 
plicability of the results presented. 


Torsional Rigidity of Stringers 


A cursory examination of Eq. (6b) indicates that as 
the critical stress of the panel increases, a reduction in 
restraint coefficient is realized. This trend is confirmed 
by data presented in Figs. 1 and 2. 

The variation of restraint coefficient for the test 
panels is depicted in Fig. 3. The restraint coefficient 
is plotted against b/t and was obtained by cross-plotting 
data from Fig. 2 with Fig. 3 of reference 3. Each curve 
has a cutoff at the b/t corresponding to critical stress of 


13,000 Ibs. per sq.in. above which the discontinuity 
caused by the primary modulus of elasticity is elimi- 
nated. 


EFFECTIVE WIDTH 
Action of Sheet 


In the typical sheet-stringer panel three distinct 
phases of sheet action may be observed. The first phase 
occurs in the initial loading stages when buckling has 
not yet begun and consequently the panel can be con- 
sidered isotropic. In the second phase, the stress dis- 
tribution varies across the panel after the sheet has 
buckled. In addition, the stringers damp out buckling 
at the line of attachment of sheet to stringer and thus 
confine instability to the bays between stringers. 

The third phase of sheet activity is concerned with 
the design of the sheet-stringer attachments when the 
stringers do not completely damp sheet buckling. In 
panels of heavy skin and relatively light stringers, fail- 
ure of the attachments in tension caused by excessive 
loads from the sheet buckles may result in premature 
failure of the panel. More frequent is the occurrence 
of interrivet buckling which concludes the ability of the 
sheet to carry additional compression load. The vari- 
ous phases of activity are depicted graphically in Fig. 4 
as a function of strain at the edge of the sheet, «,. 


’ Effective Width Equations 


An efigctive width equation can be derived directly 
from the plate instability expression, Eq. (2). The 
effective width of sheet acts at an edge stress that cor- 
responds to the strain in the stringer. Hence, the edge 
stress can be considered to be uniformly distributed*® 
over a nonbuckled panel of width 2w,: 
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o, = KE(t/2w,)? (7) 
Solving for w, 
. w, = 0.5t WKE/e, 
and letting ed 
c= 05 VK (8) 
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Fic. 4. Various phases of sheet-stringer activity under com- 
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Fic. 5. Coefficients for effective width equation. 





Ww, = ct VE/o, (9a) 


By taking the ratio of Eq. (7) to Eq. (2), another form 
of the effective width equation is obtained. 


Me = 0.5 \" (10) 
b Ce 


Eq. (10) was first derived by von Karman.’ 

Eqs. (3a, b) are analytic expressions for the plate 
instability coefficients of panels stiffened by Z- and hat- 
section elements. These equations can be introduced 
into Eq. (8) to find effective width expressions for panels 
with various stiffening elements. Thus, for Z stiffened 
panels, 


we = 0.73(b/t) °-9% WE/o, (9b) 
and for hat stiffened panels 
Ww. = 0.88(b/t)°"% VE/o, (9c) 


The variation of effective width coefficient, c, is shown 
in Fig.5. Included is the ANC-5® value of the coefficient 
for comparison. 
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Panel Tests Results 


Simultaneously with the measurement of critical 
compression stresses on stringer panel groups 3, 4, and 
5 of Table 1, effective width measurements were made 
with SR-4 strain gages cemented to each stringer, 
These data are presented in Figs. 6 and represent a 
range of b/t values between 45 and 125. For each test 
point, o,,, was measured directly and w,/b was com- 
puted from test data. Included for comparison is the 
theoretical effective width curve calculated from 
Eq. (10). 

The interrivet buckling stress for each panel was 
computed by the method of Howland” and all effective 
width test data at edge stresses greater than the inter- 
rivet buckling stress were omitted from ‘Figs. 6. 

Several effective width theories have been postulated 
at various times and it was thought desirable to check 
both the theory and test data of this paper with this 
previous work. Fig. 7 is taken from a paper by Hoff" 
to which several additions have been made. Also 
plotted is the b/t range of the effective width tests pre- 
sented in Figs. 6 and the theoretical expression, Eq. 
(10). 

Of all the theories shown, Eq. (10) apparently gives 
the best fit to the test data and is still predominantly 
conservative throughout the range of the data. 


Effective Width Above Proportional Limit 


In determining effective sheet width when the edge 
stress exceeds the proportional limit, a close approxi- 
mation to the actual effective width is obtained by using 
a strain analysis. Consider the stress-strain state of a 
buckled sheet. shown in Fig. 8. To the left of the 
centerline the strain is indicated while the right ordinate 
is a stress axis. 

The use of the modulus of elasticity above the pro- 
portional limit results in the unconservative estimate of 
effective width shown as (1) in Fig. 8. The ‘‘exact” 
solution (2) is obtained above the proportional limit if 
the stress corresponding to the strain shown at the left 






























































Fics. 6. Comparison of effective width theory, Eq. (10), and test data. 
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EFFECTIVE WIDTH OF FLAT PLATES 


of the centerline is determined from the stress-strain 
curve of the sheet material. Use of this stress distri- 
bution in determining effective width leads to a graphi- 
cal integration process, however. 

A good approximation to the latter method is ob- 
tained by finding the edge stress corresponding to the 
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Fic. 7. Comparison of experimental data with various theories. 





(1) Cox, theoretical, simple 
support. 

(2) Cox, theoretical, clamped. 

(3) Cox, experimental. 

(4) Marguerre. 

(5) Wagner, edges free. 


(6) Wagner, edges fixed. 
(7) ANC-5. 

(8) Sechler. 

(9) Marguerre, cubic. 
(10) Eq. (10). 
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edge strain from the stress-strain curve and then using 
the usual effective width equations. Method (3) is 
conservative relative to the “‘exact’’ method by the 
shaded area shown on the figure. 


Effective Area Nomograph 


Calculation of effective area by use of Eqs. (9a, b, c) 
has been simplified by construction of an effective area 
nomograph, Fig. 9. The nomograph is entered on the 
left knowing the edge strain. The sheet stress can 
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Fic. 8. Effective width above proportional limit. 











(1) Assuming linear theory. 

(2) Using stress-strain curve. 

(3) Determining o, from stress-strain curve—then using 
linear theory. 
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Fic. 9. Effective area nomograph for aluminum alloys. 
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then be found by use of the typical-stress-strain curves 
drawn, or else the nomograph can be entered at the 
abscissa with this knowledge predetermined. Further 
operations required to determine the effective area for 
a particular panel geometry are indicated in the ex- 
ample included on the nomograph. 


CONCLUSION 


A least squares analysis of the test data results in 
analytic expressions that represent mean values rather 
than lower limits of data. Essentially, these data are 
used to determine variations of the physical constant 
for the well-established plate instability equation, Eq. 
(2). The data are in accord with established theoretical 
concepts and are not employed to postulate new 
theories. Consequently, it is felt that the use of a 
statistical procedure in analyzing the data is justified. 

Furthermore, initial buckling did not result in cata- 
strophic failure of the panels tested. In general, the 
unstable sheet retains the capacity to carry additional 
loads which is insurance against slight errors in estimat- 
ing the critical stress. . 

Above critical stresses of 20,000 Ibs. per sq.in., 
scattering of test data from the least squares curves is 
small. At lower stresses, scattering is partly attributed 
to the experimental errors associated with testing of 
relatively thin sheet structures. In general, relatively 
small scatter throughout the critical stress range and the 
consistent trend of the data indicate that in current 
structures the choice of the stringer to be used with a 
particular sheet b/t ratio is governed by considerations 
that reduce to narrow limits the range within which the 
stringer torsional rigidity may be selected. Therefore, 


the use of the results presented herein for panels that 
are radically different from those presented in Table 1 
is cautioned. 
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Theoretical Principles of Wing-Tip Fins for 
Tailless Airplanes and Their Practical 
Application 


MAURICE A. GARBELL 


Consultant in Aeronautical Engineering and Meteorology 


SUMMARY 


The paper regards the relative merits of wing-tip fins and cen- 
tral vertical fins in their comparative contribution to the static 
directional-stability derivative, the rotary damping in yaw, and 
the side-force slope of tailless aircraft. 

A detailed analysis reveals the manner and degree in which the 
lift and drag of wing-tip fins concur in producing stability and 
damping. Thé measure in which toed-in fins increase the fin 
effectiveness and the useful stable range of angles of sideslip are 
critically examined, with particular consideration afforded the 
drag penalty induced by this arrangement. 

The effect of wing sweepback on the various types of fin con- 
figurations is discussed. The difficulty of obtaining adequate 
fin effectiveness on wings without sweepback, especially if the 
fin drag constitutes a critical design factor, is examined. 

New thoughts are suggested which permit the attainment of 
higher fin effectiveness than heretofore achieved, at a very low 
drag penalty even on wings with no sweepback; all of the reason- 
ings are substantiated by encouraging experimental results. 


SYMBOLS 


The standard N.A.C.A. nomenclature and symbols have been 
adapted for the special purposes of this paper. 


aj = slope of lift curve of fin airfoil (at infinite aspect ratio) 

A = geometric wing aspect ratio (= b?/S) 

Ay; = geometric fin aspect ratio (= ;?/S;) 

b = wing span (ft.) 

by = fin span (ft.) 

Cy, = fin drag coefficient, referred to the fin area Sy 

Cpp = profile drag coefficient, referred to the wing area S 

Coy, = fin profile drag coefficient, referred to the fin area Sy 

Cr = wing lift coefficient 

Ci, = fin lift coefficient, referred to the fin area Sy 

Cry, = fin lift coefficient at zero sideslip (fin preload coeff- 
cient) 

Cm = wing pitching-moment coefficient about the center of 
gravity of the aircraft 

Cm, = wing pitching-moment coefficient about the aero- 
dynamic center of the wing . 

Cy = airplane yawing-moment coefficient dbout the aero- 
dynamic center of the wing 

Cv¢ = yawing-moment coefficient contributed by the fin(s) 

Cy = side-force coefficient of the aircraft, along its Y-axis 

D = drag (Ibs.) 


span efficiency of the wing 

span efficiency of the fin 

fore-and-aft distance between the aerodynamic cen- 
ters of the wing and of the vertical fin (ft.) 

N = yawing moment (ft.lbs.) 

Ny = yawing moment contributed by the fin(s) (ft.Ibs.) 

r = yawing velocity (rad. per sec.) 


e 
ef 
ly 


Presented at the Design Session, Fourteenth Annual Meeting, 
I.A.S., New York, January 29-31, 1946. 


> = wing area (sq-ft.) 
Sy = area of one fin (sq.ft.) 
V = true air speed (ft. per sec.) 
B = angle of sideslip (rad.) 
€8 = angle of sidewash (rad.) 
nm = relative fin efficiency” 
p = standard sea-level air density (0.002378 slugs per 
cu.ft.) 
y = angle of yaw = —8(x/180) (degrees) 
Subscripts: 
f = pertaining to fin(s) or contributed by fin(s) 
l = pertaining to left fin or contributed by left fin 
r = pertaining to right fin or contributed by right fin 
w = pertaining to wing or contributed by wing 
INTRODUCTION 
Purpose of Study 


HE MAJOR PROBLEM that has hindered the general 
T acceptance of tailless aircraft throughout the 
history of aviation is their unsatisfactory lateral oscil- 
latory stability. 

Other stability problems of the tailless airplane— 
namely, the longitudinal phugoid stability, the longi- 
tudinal short-period oscillatory stability, the rolling 
subsidence, the spiral stability, and the general stalling 
characteristics—have been solved satisfactorily; yet, 
the lateral oscillatory behavior during combined rolling 
and yawing oscillations in virtually all practicable de- 
signs of tailless aircraft developed to date has not been 
sufficiently satisfactory to warrant the confidence and 
decisive approval of the final judges—the pilot and the 
customer. 


Longitudinal Phugoid Oscillation 


Briefly, the adequacy of the tailless design with re- 
spect to the longitudinal phugoid oscillation may be 
attributed to the reduction of the longitudinal moment 
of inertia (shortening of the pitching radius of gyration) 
due to the elimination of the fuselage and tail surfaces. 

Static longitudinal stability may be achieved without 
a tail by employing a wing with a positive pitching 
moment coefficient Cy, about the aerodynamic center 
and by locating the airplane center of gravity ahead of 
the wing aerodynamic center; the only penalty in- 
curred through the omission of the horizontal tail is a 
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reduction of the rotary damping in pitch of the aircraft 
by approximately 75 per cent. 

A comparison of the moments of inertia shows that a 
similar reduction of 75 per cent or more is attainable 
through careful design and allocation of weights, so 
that the nondimensional damping-in-pitch factor ap- 
pearing in the biquadratic dynamic-stability equations, 
which is proportional to the aerodynamic damping and 
inversely proportional to the moment of inertia, re- 
mains substantially the same for tailless airplanes as 
in conventional designs, with comparable values of the 
static-stability derivative dCy/dC,. 

Favorable also for the tailless airplane is the fact 
that its longitudinal oscillatory response to a vertical 
gust is less violent than that of the conventional air- 
plane, because the tailless is free from the pitching mo- 
tion that is produced in the conventional airplane by the 
time lag intercurring between the respective impinge- 
ments of the gust on the wing and on the tail. 


Longitudinal Short-Period Oscillation 


Considerations of static longitudinal stability and 
damping in pitch similar to those evolved for the phu- 
goid oscillation indicate the basis for the attainment 
of short-period oscillatory stability in tailless air- 
craft. 

Two peculiar problems in connection with short- 
period longitudinal oscillations are also susceptible to 
relatively simple solutions: 

The problem of the binary short-period oscillation, 
or ‘‘wallowing,’’ involving both the entire craft and the 
hinged wing-trailing-edge elevators, has been satis- 
factorily solved by reducing to zero both the elevator 
mass unbalance and the elevator hinge-moment de- 
rivative with respect to the angle of attack. 

The dangerous characteristic of certain tailless de- 
signs that exhibit negative values of the damping in 
pitch, through trailing-edge airflow separation and wing- 

tip stall during violent nose-up pitching movements 
while the airplane is trimmed for high C;’s, and that at 
near-neutral static longitudinal stabilities (aft centers 
of gravity) lead to an undamped positive pitching rota- 
tion of the airplane known as “tumbling,” may be ef- 
fectively counteracted with comparative ease by in- 
creasing the static longitudinal stability through a for- 
ward displacement of the center of gravity. 


Rolling Subsidence 


The two principal factors that determine the rolling- 
subsidence characteristics—that is, the rotary damping 
in roll and the rate of change of rolling-moment coef- 
ficient with sideslip—are essentially functions of the 
plan form and dihedral of the wing. Therefore, they 
are not affected appreciably by the elimination of the 
tail surfaces; hence, tailless aircraft exhibit virtually 
the same rolling-subsidence characteristics .as do 
conventional airplanes with comparable wing geome- 


try. 





Spiral Stability 


Specifications exacted by pilots and customers usu- 
ally minimize the importance of this criterion, unless 
an airplane is severely spirally unstable. Despite the 
exiguity of their damping in yaw, tailless airplanes as a 
rule exhibit slightly positive or neutral spiral stability 
because of their extremely small static directional sta- 
bility and their comparatively large rolling-moment 
derivative with sideslip (required for good rolling- 
subsidence and handling-near-the-ground characteris- 
tics). 


Stalling Characteristics 


Although stalling characteristics are primarily inher- 
ent in the wing design, the absence of a tail further com- 
plicates the problem of longitudinal and lateral stability 
during the stall, especially in the case of sweptback 
wings. ; 

However, with correct design and appropriate selec- 
tion of airfoils! it is possible to locate the stall-inception 
area in the inboard portion of the wing so that, even in a 
sweptback wing, the pitching moment at the stall be- 
comes decidedly negative (nose-down), thus rendering 
the static longitudinal-stability derivative fully stable 
throughout the entire stall (Fig. 11 of reference 1). 


Lateral Oscillatory Stability 


The unsatisfactory damping of the combined rolling 
and yawing oscillations of tailless airplanes remains 
as the only justifiable objection to the flight character- 
istics of this type of craft. 

Although the oscillations usually are convergent, the 
large oscillatory amplitudes and the great number of 
complete (5- to 10-sec.) cycles required to damp the 
initial amplitude to a tolerable value have heretofore 
limited the usefulness of the tailless airplane for passen- 
ger comfort, precision bombing, or precision gunning. 

An examination of the individual terms appearing in 
the factors of the general biquadratic lateral-stability 


equation reveals that the aerodynamic and mass terms - 


referring to the axis of roll (X-axis) are virtually of the 
same magnitude as the corresponding terms of conven- 
tional airplanes but that the terms referring to the axis 
of yaw (Z-axis) are greatly altered in the adverse sense 
of reduced stability. 

The balance of reduced static stability and rotary 
damping versus an equally reduced moment of inertia, 
which establishes good dynamic Jongitudinal-stability 
characteristics in skillfully designed tailless airplanes, 
has no counterpart in the case of the directional terms. 

The static directional-stability and rotary damping- 
in-yaw derivatives of most tailless airplanes are sub- 
stantially smaller than the corresponding derivatives 
of satisfactory conventional airplanes. Unfortunately, 
this adverse decrease in static stability and damping is 
only partly balanced by a concurrent reduction in the 
yawing moment of inertia, especially on multipropul- 
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WING-TIP FINS FOR TAILLESS AIRPLANES 


sion airplanes where the heavy power units located on 
the wing, together with the heavy wing structure it- 
self, are the principal contributors to the yawing mo- 
ment of inertia. 

Another unfavorable factor is the unbalance between 
the small static directional-stability derivative and the 
comparatively large rolling-moment derivative with 
sideslip (required for good rolling-subsidence, spiral- 
stability, and handling characteristics), which further 
reduces the damping convergence of the lateral oscilla- 
tions and produces the so-called ‘“‘Dutch Roll.” 

Summarizing, the inadequacy of the lateral oscilla- 
tory-stability characteristics in niost tailless aircraft is 
due to the following three individual] causes, listed in 
their order of importance and-effectiveness: 

Static Directional-Stability Derivative—A_ typical 
value of the static directional-stability derivative is 
dCy/d8 = 0.005 for an all-wing tailless airplane, as 
compared with a desirable value of at least 0.030 (com- 
parable to approximately 0.060 on a conventional de- 
sign).* 

Rotary Damping Dertvative in Yaw.—A characteris- 
tic value for a typical all-wing tailless airplane is 
dCy/d(rb/2V) = —0.002 at a high-speed C, and 
—0.005 at a cruising C,;, as compared with desirable 
values of —0.005 and —0.010, respectively (comparable 
with —0.015 and —0.030, respectively, on conven- 
tional designs).* 

Side- Force Derivative with Sideshp.—This derivative 
affects the character of the lateral oscillatory stability 
only to a minor degree but has a more appreciable ef- 
fect on the linear amplitude of the horizontal displace- 
ments perpendicular to the intended course and on the 
handling characteristics of the craft during controlled 
maneuvers, with a value of dCy/d8 = —0.010¢ of the 
all-wing tailless airplane, as compared with a desirable 
value of —0.300 for both the tailless and the conven- 
tional designs. 

It is possible to achieve adequate lateral oscillatory 
stability fully in tailless airplanes through the use of 
automatic directional controls actuated by displace- 
ment, velocity, and acceleration pickups. 

However, the attainment of inherent stability and 
damping, independent of any control instrument or 
apparatus, is worthy of earnest engineering effort. 

The present study analyzes the various former at- 
tempts to improve the foregoing three critical charac- 


* The difference between the values defined as desirable for the 
tailless designs and those quoted for conventional airplanes is 
justified in two ways: (1) the appreciable decrease in yawing 
moments of inertia (though not comparable to the reduction of 


the pitching moment of inertia) of the tailless airplane; (2) an‘ 


assumed appreciable increase of the “‘effective’’ damping in yaw 
which may be attained through appropriate adjustment of the 
cross derivatives—that is, of the rolling-moment derivative due 
to yaw and of the yawing-moment derivative due to roll. 

t This typical value refers to the flaps-up configuration only. 
With trailing-edge flaps deflected even positive (unstable) values 
are encountered. 
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teristics by means of wing-tip fins and other types of 
vertical fins. A new method supported by experimen- 
tal proof is introduced demonstrating how the desired 
result has now been accomplished. 


DISCUSSION 


The static directional-stability derivative, the rotary 
damping in yaw,f{ and the side-force derivative con- 
tributed by various fin configurations on tailless air- 
planes are systematically examined in the following dis- 
cussion. A further quantitative comparison of the 
most significant configurations is shown in the Synop- 
tic Table on pages 528 and 529. 


Parallel Wing-Tip Fins with Symmetrical Airfoils on 

Wings Without Sweepback 

The most elementary application of wing-tip fins is 
that of parallel fins with symmetrical airfoils on wings 
without sweepback (Fig. 1). 

In unyawed flight the fins operate at nearly zero lift, 
assuming the airflow to be strictly two-dimensional, 
even at the wing tips. Thus the fins add only their 
profile drag to the overall drag of tie airplane. How- 
ever, the end-plate effect** of fins with aspect ratios of 
1.5 to 2.0 is usually found to be sufficiently pronounced 
to nullify this additional fin profile drag through a re- 
duction of the induced wing drag (by increasing the ef- 
fective wing aspect ratio) at C,’s of 0.6 and higher. 

In a sideslip, the two fins operate at the same angle 
of attack; hence, their lift and drag forces are equal. 
The resulting moment about the wing aerodynamic 
center and, therefore, the contribution of this type of - 
wing-tip fins to the static directional stability of the craft 
is zero. 

A simple approximation indicates the extent to which 
the wing-tip fins contribute to the rotary damping in 
yaw. If it is assumed that the craft is yawing at a 
constant angle of zero sideslip and at a constant angular 
velocity 7, a damping moment JN; is produced by the 
disparity of the respective drag forces acting on the two 
fins. Its magnitude is expressed by 


N; = (D, — D,)(b/2) (1) 


and, assuming the same profile drag coefficient for either 
fin, 


t The values of the rotary damping in yaw, in this discussion, 
refer to a yawing motion at constant zero angle of yaw (curved 
flight path). The importance of this conservative assumption has 
been emphasized by Multhopp.”: ® 

** References 3, 4, 5, and 6; also author’s unpublished data. 
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DESIGN i VS. DAMPING SIDE-FORCE 
re ANGLE OF SIDESLIP IN YAW SLOPE 
OVERALL DRAG INCREASE 
| WING ALONE 
A= 10 (REF. 12) dc 
C Cp dp 
d = Dp , “Di 
Cpp = 0.0060 SCN wo ~ oa be ee A 
dp | N d5y 
e = 0.85 
dCy dc 
= 0.005 N TOTAL Sa 
dp a az 
SS 2V 
| AT CL =0.1 
-2.0016 
AT C, =06 
DRAG BASIS - -0.0060 -0.0! 
2 PARALLEL WING-TIP FINS WITH EQUATION 
SYMMETRICAL AIRFOILS ON WINGS A = 10 EQU.«5) ss 
WITHOUT SWEEP-BACK dCn 
i —N\ =-o0010 
Cop 0.0060 (is 
ic 2avit 
e = 0.90 | ae - 
- 
== FS= Sp = 0.05S 5 
Ap = 2.0 
! dc 
+— b | Tota, —N 
“dp, = 0.0100 dt 
‘is ey, = 0.70 AT C, =O1 
4Cpp (SUBSTANTIALLY SAME AS WING ALONE; ~0.0026 
AT C,=0.1: 0.0010 Nt = 0.90 POSSIBLE SLIGHT INSTABILITY AT ZERO AT C,=0.6 
AT C,=0.1: 0.0002 SIDESLIP) - 0.0067 Jaa 
3 CENTRAL VERTICAL FIN 
EQU.(9) 
A =10 
Ce =- 0.0028 
. rT 
Cop = 0.0060 2v/¢ 
EQU. (8) 
e = 0.85 dCy N 
» _ FIN STALL 
re ea ee Sp = 0.10S _—9 
! oo j p 
eating, Sang ome: As = 1.5 +15° 
Dp, = 0.0100 
4v =0.07 ‘i ‘ 
“he dc 
ad TOTAL — 
<8 = 0.035 are 
dp 2v 
FOR 2° WING , AT C_ =O! 
. DIHEDRAL] (AT WING ANGLES OF ATTACK NEAR STALL L’"-0.0044 
THE FIN EFFECTIVENESS IS DANGEROUSLY | ,_ 
AT ANY CL: &Cp, = 0.0010 REDUCED THROUGH BLANKETING BY “104.0085 
P SEPARATED AIRFLOW. -0.28 
4 PARALLEL WING-TIP FINS WITH 
SYMMETRICAL AIRFOILS ON EQU.13) 
SWEPT-BACK WINGS. dCy 
A= 10 —N) =-0.0045 
d 
Cop, =. 0.0060 2v/f 
o. <'\998 EQu.cin 
Sp = 0.05S ley =0.024 N 
dp av 
As = 2.0 
le b »| Cp = 00100 USUALLY Pp 
Pr FLAT SPOT 
WV.007 e, = 0.70 
b 0.0 f 
Np = 0.90 
dc 
TOTAL —N 
Y= Sogh a 
: 2v 
A 
+ = 05882 AT C_=01 
- 0.006! 
4Cop 
AT CL =0.1 0.0010 AT CL=0.6 
AT CL =0.6 0.0002 -00102 ~027 
Nore: All calculated data have been fully substantiated by wind-tunnel tests. 
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TABLE 1 (Continued) 












































SYNOPTIC TABLE 
CONTINUED 
FIN ARRANGEMENT 
sibiden YAWING MOMENT COEFFICIENT DAMPING SIDE-FORCE 
vs. IN YAW SLOPE 
PARAMETERS ANGLE OF SIDESLIP 
OVERALL DRAG INCREASE ; 
5  TOED-IN WING-TIP FINS WITH 
SYMMETRICAL AIRFOILS ON WINGS EQU.(25) = 
WITHOUT SWEEP-BACK A= 10 dCy =-enoes 
gre 
Cpp = 00060 avit 
EQU.18) 
= 088 dCn Ic 
sllTOE-IN Cc, =0. oll og- ” ==N =0.006 N 
at Lp, = 048 S508 IN ae AY 2Eno 
\ , | Sp = 0055 SIDESLIP 
= SS Ap = 10 p 
~—b a 
‘a Cop, = 0.0100 ToTaL SCN 
f gzb 
2V 
er = 0.70 
AT C, =0.1 
4Cpp Ne = 0.90 ~ 0.004! 
AT C.=0.1 0.0025 AT C, =06 
AT C_ =06 0.0020 - 0.0084 apne 
6 TOED-IN WING-TIP FINS WITH 
SYMMETRICAL AIRFOILS ON SWEPT- EQU(36) 
BACK WINGS ¥ 
rb 
9 vl¢ 
£QU.(27) AT C, <0. 
dCn ic -0.0063 
——- 7 N 
se||TOE-IN C4 =0.26 — Se\|TOE-IN ap ae zEA0 AT C06 
+ _ \| SIDESLIP -0.0067 
vt p 
4 lv. 9.07 
| | 
~——b 2 > rota. SCN 
gz 
2v 
AT C =01 
ACpp -0.0079 
AT C_=0.1 0.0025 AT C =06 
AT C_=06 0.0017 -00124 -0.27 
WING-TIP FINS WITH INBOARD A = 10 
CAMBER +ZERO EFFECTIVE TOE-IN EQu«s) 
7A WINGS WITHOUT SWEEP-BACK| Cpp = 00060 SCn\ 0010 
gre 
e = 090 avit 
£QU.(37) 
S, = 0.05S 
J <EN- 0020 Cn 
Ay = 2.0 P 
| “Op, = 0.0100 B 
— 
» e, = 0.70 i 
TOTAL —= 
Nr * 0.90 d5y 
J AT CG =01 
ACpp = soseer L "= 0.0026 
AT CL=01 0.0010 r = 05842 AT C, =06 
AT Ci, =06 0.0002 - 0.0067 -0.27 
7B SWEPT-BACK WINGS 
EQUS.11308(38) 
c 
—N|\ =- 0.0047 
dzB), 
£QU.(37) 2v 
SCN. 0.040 N 
dp 
Bp 
dc 
Tota. —N 
arb 
2V 
AT. C, =01 
4Cp, -00063 
AT C_=0.1 0.0010 AT C =06 
AT C/=06 0.0002 “0.0104 -0.27 
























































Note: 


All calculated data have been fully substantiated by wind-tunnel’ tests. 








ll 


Ny = Cop, S,5 (V2 — Ve) 5 


= Con, 5,8 (V — av)? - (V+ AVS 


= —Con, 5,5 (4VAV) : (2) 
Substituting AV = r(b/2) and rearranging, 
p rb 
Ny; = —2CovbS; 5 V? OV (3) 


or, in nondimensional form, 
Cuz = —2Cdp, (S/S) (rb/2V) (4) 


The differentiation with respect to (rb/2V) yields the 
nondimensional damping-in-yaw derivative: 


[dCy/d(rb/2 V)], = —2Cop, (S;/S) (5) 


The contribution of the wing-tip fins to the side-force 
derivative with sideslip is determined to be: 


dCt) a [ ($6) con p+ ($2) sn 0] S 
a) > 2| (4 cos B + dB sin B 5 (6) 


and, within the limits of practical accuracy, 
(=) a -o Sx) = — as eae (7) 
dp /y dp /S 1 + (ap/aA,) S 


Central Vertical Fin 








Even on tailless airplanes with no wing sweepback, a 
vertical fin may be conveniently located at some dis- 
tance aft of the wing aerodynamic center by placing it 
at or near the vertical plane of symmetry (XZ-plane), 
where the central portion of the wing offers the neces- 
sary structural strength for the support of the fin close 
to the wing trailing edge, an advantage that cannot 
be had at the wing tip (Fig. 2). 














Fic. 2. 


The static directional-stability contribution of the 
central vertical fin is expressed by the following equa- 
tion in which the extremely small effect of the fin drag 
has been disregarded: 


dCy deg\ ly Sy, D4 

at) whic Et} —2-——_ (8) 

dp dB/ b S 1+ (ap/7Ay) 
where the term deg/d@ is the derivative with respect to 
the angle of sideslip of the sidewash due to the effective 


dihedral of the wing,’ which, in the case of a vertical fin 
located at the wing trailing edge and in close proximity 
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to the vertical plane of symmetry, assumes some im- 
portance in reducing the stability efficiency of the fin. 

The rotary damping in yaw contributed by a central 
vertical fin is expressed by® 


dCy ly . af S; 
tetas € 

d(rb/2V) Jy b/ 1+ (a,/rAy) S 
It is herein assumed that the yawing rotation occurs 
at constant zero angle of sideslip, thus keeping the side- 
wash and also the sidewash derivative with respect to 

(rb/2V) equal to zero. 

The side-force derivative of the central vertical fin 
finds an expression similar to that of the wing-tip fins 
[Eq. (7)] but with due consideration for the sidewash: 


(=) (1 “) Qf Sy 

—}) = -» —. ae ae 

dp /, r dB] 1+ (ap/nAy) S 

Parallel Wing-Tip Fins with Symmetrical Airfoils on 
Sweptback Wings 


(10) 


Parallel fins with symmetrical airfoils placed at or 
near the tips of sweptback wings combine the respective 
aerodynamic effects of similar wing-tip fins on wings 
without sweepback and of vertical fins having a fore- 
and-aft “‘tail length.” 

In a sideslip the fins operate at the same angle of at- 
tack; hence, their lift and dtag forces are equal. 

A resolution of the resultant aerodynamic forces in 
chordwise and normal fin-force components (Fig. 3) 
reveals that the moment of the equal chordwise force 
components about the wing aerodynamic center is zero 
and that the rate of change of the yawing-moment coef- 
ficient with angle of sideslip due to the fins (assuming 
then fin normal force to be practically equal to the fin 
lift force) is: 


ee 
dB), ~~ b 1+ (@y/eAy) S 


where the term deg/d8 appearing in Eq. (8) has been 
omitted because the rate of change of the sidewash with 
angle of sideslip is nearly zero at the wing tip. This is 
the total contribution of the fins to the static directional 
stability of the airplane; it is wholly due to the fore- 
and-aft tail length of the fins. 


(11) 
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The rotary damping in yaw of the sweptback parallel 
wing-tip fins is conveniently obtained by resolving the 
motion of the fins, while yawing, into X and Y compo- 
nents. The following reasoning leads to an expression 
for the damping in yaw (bearing in mind the definition 
of S; as the area of ove fin). 

From the previous derivations of Eqs. (5) and (9) 
it is readily seen that 


dCy ly Or Sy 
i ~2[ Coy A 2n/(F i + rae 
(12) 
where 
Coz = Cony + (Cr?/rA fey) 
and 
a rly Q4o 
oS gs (ay,/7A ;) 


hence, substituting 


7 rly 2p —) 2s 
ey ( fe TAs + dp) ey 
in Eq. (12): 
dCy ] | ( af y mAy 
on 3 Coy 7; tAy+ ay) e + 


ly azo Sy 
, u(G ) 1+ maa (13) 


It is apparent that the middle term in the bracket 
refers to the induced fin drag. This term is dependent 
on both the yawing velocity 7 and the air speed V, and 
for relatively slow yawing oscillations its magnitude be- 
comes negligible at high air speed. 

The contribution of the fins to the side-force deriva- 
five remains as expressed by Eq. (7). 


Toed-In Wing-Tip Fins with Symmetrical Airfoils on 

Wings Without Sweepback 

Wing-tip fins on tailless airplanes with little or no 
wing sweepback have been toed-in by some designers 
in order to obtain a fin contribution to the static direc- 
tional stability. The following investigation reveals 
the inefficiency of this arrangement. 

At zero sideslip the fin preload coefficient Cry, due to 
the toe-in induces an additional drag increment 

(ACop)+ = 2(Cz?/7A ze) (S;/S) (14) 
which, in all practical cases, equals or exceeds the drag 
of a fuselage and conventional tail and thus partly in- 
validates the prime justification for the tailless 
design. 

Further analyzing the contribution of toed-in fins 
to the static directional stability of the craft, it is found 
from an inspection of Fig. 4 that, at an angle of sideslip 
8, the preload coefficient Crs, produces a destabilizing 


moment, whereas the induced fin drags contribute a 
stabilizing moment. 

In fact, the yawing-moment coefficient Cy; at an 
angle of sideslip 8 may be expressed by the equation 





Cy, = 

V2. 2b Cry? — Cry? 2 
ae S;- (- 2Cx, sin B + % “Boos 2) /* bS 
Z y WA sey 2 


| —2C1y,, sin B + 
(—Cty, — ACt,)? — (Cry, — 





A 2 
CLy) on |= Sy 


WA sey 2S 
' 4CLi_ ACL, ] Sy 
= = 9 —_——____ oma = 
| 2CLy, sin B + mo cos B 35 (15) 
and substituting 
ayo > 
ACL, = 1, 16) 
0 "TE (an/aAy) 
Eq. (15) is transformed into 
. 2CL 4 ns1 Jo \ os 
=( — eeatetes . Deon = (17) 
Cny ( Cry, sin B + id + Gale B cos B 3 (17) 
and differentiating Eq. (17) with respect to 8, 
@),L 
— } =| —Cz,, cos B + 
& t : 
2CLy 11ly , Sy 
—_——— s — (18 
(aA, + &)e; (cos 8 + £ sin B) 5 (18) 


An inspection of the terms contained in the brackets 


of Eq. (18) shows that, in order to achieve a positive 


directional-stability contribution from the fins, it is 
necessary that 


2nay (1 — 6 tan £) 
(Ay + dp)ey 


and solving Ineq. (19) for Ay, 


>1 (19) 





A;< < “82% (1 — B tan ~) — | (20) 





TOE-IN=| 


“yr roe-m 
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>——— 
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At zero sideslip this condition is fulfilled if 
Ay < (ap/7)[2(n,/e7) — 1 (21) 
and, assuming that n; approximately equals ey, 


A; < dy/m = approx. 2 (21a) 


If a stable range of angles of sideslip of + *'/: rad. 
( + 28°) is desired, then 


a 
A;,< “#(1.5 ~~ 1) = approx. 0.5— = approx. 1 
T 


T ey 
(22) 


This development shows that, regardless of the mag- 
nitude of the fin preload (within the limits of aerody- 
namic linearity), the general character of the fin contri- 
bution to the static directional stability is a function of 
the fin aspect ratio only. 

By rearranging Eq. (18), to read 


dCy = | ( 2nsA 7 -: )- 
@ ) ihre! nang YO Vr Sn 


2 S 
B sin B Ns Jo | oT 
f. 


CA;+ayis ™ 


it becomes conversely apparent that, for a given fin 
geometry and a given angle of sideslip, the quantita- 
tive magnitude of the static directional-stability con- 
tribution is directly proportional to the fin preload 
(toe-in). 

Furthermore, by differentiating Eq. (18) with respect 
to the fin aspect ratio and equating to zero, it is found 
that 





8(dCy/dB)y _ 
5A, 
2n sO, ” . T #4 is 
Tg CLy (cos B — B sin 8) (eA, + a,)° +a, = 0 (24) 


For small angles of sideslip cos 8 + 8 sin 8; hence, 
both solutions of Eq. (24) are Ay = ©. An inspection 
of the subsequent derivative, as well as previous con- 
siderations, characterizes this double solution as a 
minimum of dCy/dB. It is therefore concluded that— 
all other factors remaining equal, notably the fin pre- 
load coefficient Cx, (not the angle of toe-in)—the static 
directional stability increases with decreasing fin aspect 
ratio (at the cost of a linearly growing induced fin drag), 
and there is no optimum fin aspect ratio. 

Analogous considerations show that, in a fin configu- 
ration with negative toe-in (so-called foe-out), the fin 
aspect ratio must exceed an approximate critical value 
of 2 in order to achieve a positive static directional- 
stability contribution within a sideslip range of +'/2 
rad., including the angle of zero sideslip. 

Throughout the foregoing reasoning it has been as- 
sumed that the lift-curve slope of fins with extremely 
small aspect ratio remains constant up to '/, rad. 
geometric angle of attack. The two alternate solu- 


tions offered—toe-in with a fin aspect ratio less than 
unity, or toe-out with a fin aspect ratio greater than 2— 
are not fully equivalent, however, because at angles of 
sideslip greater than the angle of fin stall the sharply 
increasing drag of the leading fin in the case of toe-in 
remains stabilizing, whereas the corresponding stall 
drag of the trailing fin in the case of toe-out is destabiliz- 
ing and thus produces dangerous yawing instability at 
large angles of sideslip. In this connection it should 
also be remembered that the stalling angle of attack de- 
creases rapidly as the fin aspect ratio exceeds 2. 


The rotary damping in yaw contributed by the toed-: 


in (or toed-out) fins is obtained by substituting the total 
drag coefficient to the profile-drag coefficient in Eq. 
(5), yielding 


] = ste, 2 => 
ia 


| dCy 
d(rb/2V) S 





(25) 


Nc 
= Coy + Cin \z 


1A sy, S 


The side-force derivative is not affected by toe-in or 
toe-out and remains expressed by Eq. (7). 


Toed-In Wing-Tip Fins with Symmetrical Airfoils on 

Sweptback Wings 

The characteristics of the preceding arrangement 
are here combined with those of the central vertical 
fin. 

The static directional-stability contribution is hence 
obtained by combining the terms contained in Eqs. 
(11) and (18) :* 


2CLyn A 4 


' (dCy 
(52), ite [ - Cu WRT (aA, + Apes 3 


(cos 8 — B sin 8) + dn, on te | (26) 
1 re a4 


nA, 





or, rearranging, 


dCy v7 {- 2n Af, 
= ) pS se kad wAy + ap " 


CLy, Ws i ly ) Sy (4 
| = (cos 8 — Bsin B) + rA,y ae (27) 





The variation of (dCy/d8); with fin aspect ratio is 
investigated by differentiating Eq. (27) with respect to 
this variable: 


5(dCy/dB)y _ { 
6A; 





= 2nd yo | <4 
(wAy + dp)*L ey 
=| 4 omer (28) 


(cos 8 — B sin B) + Ay > me serene Ee 


Rearranging and equating to zero, 
* The change in the stability contribution of the induced fin 


drags due to the change in drag arms with sideslip is negligible 
for practical values of the sweepback. 
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2n sy | ‘ete! es) (cos B — Bsin B) + 7A, (ly/b) _ 
tA, + ay tA; + dp 


| 
l=0 (29 
y}=0 (29) 


Considering only positive values of A; it is evident 
that Ay = @ is a root of Eq. (29), but whether this 
represents a maximum or a minimum of (dCy/d8), de- 
pends on the sign of the term in brackets which, in 
turn, depends on whether 


< Apps(ly/b) 30 
*h > cos B — Bsin B (30) 





If the < sign applies, (dCy/d8),is maximum for A; = 
©, indicating that the static directional stability im- 
proves with increasing fin aspect ratio. 

Conversely, if the > sign applies, (dCy/d8), is mini- 
mum for Ay = ©, indicating that a reduction of the fin 
aspect ratio produces an increase in static directional 
stability. 

For C1, equal to the term opposite it in Ineq. (30), the 
static directional stability is independent of the fin 
aspect ratio and assumes the value 


dv) ] _ 4 b(_, cos, S, 
IG Ln ™ “| * cos B - fan” 2m) 








(31) 
The critical value of the fin preload coefficient 
Ayes (ly /b) 
CL — — “ 30a 
[ Ly] t. cos 3 a 8 sin 3 ( ) 


for which the sense of the fin aspect-ratio effect on static 
directional stability changes sign, is well within the 
practicable range of typical values of airplane geometry 
(see the Synoptic Table). It is therefore not surprising 
that former wind-tunnel investigations have yielded 
contradictory results, or none at all, with reference to 
the effect of fin aspect ratio on static directional sta- 
bility. 
In the case of 


Cipy > [anesly/b)|/[cos 8 — B sin 8B] 


an expression for the maximum fin aspect .ratio for 
which the static stability contribution is still stable— 
comparable with Ineq. (20)—is derived by equating the 
term in braces in Eq. (28) to a value greater than zero 
—that is, 


2n Af k — B tan B 1A i] ie 
tA, + dy er Cry, cos B b 


and solving for A ,, 


2a) 22 — B tan B) — | 
er 


(30b) 





(32) 


Tv 


A; < 





33 
a 2nsA4y ly (38) 
CLy, cos B b 


and by substitution from Ineq. (20), 


2n 49 | 
A,=A i/i-—a4 34 
f tS (no ene | K CLs, cos 8 b (3 ) 


Hence, the maximum value A, for positive static 
directional stability increases from Af (no sweepback) tO in- 
finity as Cty, approaches the value 


2nay (ly/b)cos B (34a) 


Within the range of values of the preioad coefficient 
CL, defined by 


Aper(ly/b) 
cos 8 — sin B 


the sensitivity of (dCy/d8),; to changes of the fin aspect 
ratio decreases as CL, approaches the value of the right 
term in Ineq. (35); meanwhile the absolute values of 
(dCy/d8); approach that expressed in Eq. (31). For 
the range of values of Cz,, defined in Ineq. (35), of course, 
the concept of a maximum critical fin aspect ratio for 
positive directional stability loses its significance, 
because even with Ay = o, (dCy/d8); remains more 
than 0. 

The fin contribution to the rotary damping in yaw is 
readily expressed by an equation similar to Eq. (13), 
wherein the fin lift coefficient (in the induced-drag 
term) is adjusted to include the fin preload coefficient 
CL, as follows: 


Qn gp (ly /b) 
cos 8 





> fie > (35) 








dCy rly a = ’ 
Se Copy + CLy + V ™ 4 + fo 4. 
d= mA 4 
suU/ 7 A sey 
a4 
ly\? ———_- |S 
™ (F) 14/3 0) 
wAy 


The side-force contribution of the fins remains ex- 
pressed by Eq. (7). 


Wing-Tip Fins with Outboard Camber 


Lippisch and others have advocated the use of fins 
with cambered airfoils facing outboard (Fig. 5). 

Such configuration does not offer any improvement 
over that of fins with symmetrical airfoils having the 
same effective toe-in (same fin preload coefficient) inso- 
far as the stability and damping derivatives are con- 
cerned. 





Fie. 5. 
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The contention that appropriately selected cambered 
airfoils present a lower profile drag (at zero sideslip) 
than comparable toed-in symmetrical-airfoil fins, at a 
given fin preload coefficient C1,,, cannot be accepted as 
conclusive in view of the predoirinance of the induced 
fin-drag coefficient Cry = Ciy"/7A se; over the fin pro- 
file-drag coefficient Cop, 


Wing-Tip Fins with Inboard Camber 


Some attempts were also made, by Bader and others, 
to utilize cambered fin airfoils facing inboard (Fig. 6). 








ban 
| p Wing 
y oR 
ly ! 
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Pr . a 
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Most of these proposals remained unsuccessful, be- 
cause this type of configuration in general does not dif- 
fer in any substantial point from the configuratiou em- 
ploying cambered fin airfoils facing outboard; hence, 
no improvement was attained over fins with symmetri- 
cal airfoils having the same effective toe-in. 

On the other hand, a marked improvement in static 
directional stability and damping in yaw can be ob- 
tained with inwardly cambered fins if fin airfoils are 
employed which exhibit a sharp increase in the profile- 
drag coefficient at negative angles of attack just below 
the angle of zero lift, whereas the profile-drag charac- 
teristics at positive lift are normal. Such airfoils can 
actually be designed as illustrated, for example, by the 
well-known NACA 4306 airfoil (Fig. 714). 

The fins are toed-in geometrically so that they oper- 
ate at zero lift when the airplane is flying at zero side- 
slip; the fin preload or the effective aerodynamic toe-in of 
the fins is therefore zero. Hence, the fin profile drag 
is at its minimum and the induced fin drag is zero in 
unyawed flight. 





NACA4306 SECTION CHARACTERISTICS!! 


0.3 


0.2 
CL 





Q° (ANGLE OF ATTACK) 
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When the airplane is yawed, the leading fin (which 
operates at negative lift, see Fig. 6) produces a stable 
yawing moment through its abruptly increasing profile 
drag, whereas the trailing fin (which operates at posi- 
tive lift) suffers no comparable increase in profile drag. 

The induced-drag forces of the two fins are equal and 
do not create any yawing moment (for practical angles 
of wing sweepback). Hence, the principal contribu- 
tion of the fin drags to the static directional stability of 
the craft is due to the unequal fin profile drags. 

The fin lift forces act as in the preceding examples. 

The total contribution of the fins to the static direc- 
tional stability can, therefore, be expressed by 
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It has been asked, at times, whether fins of the ‘‘in- 
board-camber profile-drag type” on tailless airplanes 
with sweptback wings should not be toed-in a small 
additional amount to obtain the optimum rate of sta- 
bility increment versus drag penalty 


. A(dCy/d8) ,/ ACp, 


Quantitative studies based on equations similar to 
Eq. (27) and supported by numerous wind-tunnel tests 
indicate that the opt‘mum stability-versus-drag ratio- 
is attained at an extremely small angle of effective 
aerodynamic toe-in—of the order of magnitude of 
0.1°—and that any further toeing-in of the fins reduces 
the value of the stability-versus-drag criterion. 

It is therefore acknowledged that there is an optimum 
fin setting and that any further increase in stability is 
more economically obtained by increasing the fin area 
as suggested in Fig. 8. 

Wings without sweepback, with cambered fins facing 
inboard and with zero effective aerodynamic toe-in (of 
the zero-lift lines), have damping-in-yaw characteris- 
tics similar to those defined in Eq. (5) for parallel fins 
with symmetrical airfoils. 

On sweptback wings the local fin angle of attack due 
to the angular velocity of yaw 
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creates a difference in fin drag forces which, in Eq. 
(13) replaces the amount —2Cp,,(S,/S) by 


i E's S, 
| so | = (Cory + Com) (38) 


It should be noted here that CD, and Cppjy, vary 
with AS,and, hence, with in a nonlinear fashion. Con- 
sequently, the damping in yaw of this type of fin is not 
constant but increases steeply with augmenting values 
of the yawing velocity r. 

The side-force derivative is not altered by the adoption 
of cambered-airfoil fins and remains as expressed by 


Eq. (7). 


NOTE ON TAILLESS AIRPLANES WITH SWEPT-FORWARD 
WINGS AND TaIL-First AIRPLANES 


Considerations of various nature have influenced 
some airplane designers to resort to tailless layouts with 
swept-forward wings or to tail-first airplanes—that is, 
designs with control surfaces placed ahead of the wing. 

Without entering into a discussion or an evaluation 
of the general problems and aerodynamic characteris- 
tics of these two types of unconventional designs, the 
most expedient arrangements of directional stabilizers 
are here examined. 


Swept-Forward Tailless Airplanes 


The peculiar geometry of this type of tailless design 
(Fig. 9) discards the wing-tip fins entirely and offers an 
excellent opportunity for the installation of a central 
vertical fin from the points of view of both aerodynamic 
efficiency (a tail arm virtually as long as in conventional 
aircraft) and structural expediency. 











Fic. 9. 


Because of the relatively great tail length available 
in tailless airplanes with very low wing aspect ratio 
(Fig. 10), central vertical fins are also most opportunely 
employed in this type of design, in preference over wing- 
tip fins. 


Tail-First Airplanes 


The static directional stability and damping in yaw 
of tail-first airplanes requires, as a rule, the provision of 

















one or more vertical fins on the wing (located in the aft 
portion of the craft). Customarily, the fins are placed 
above or below the wing at a small distance from the 
fuselage (for examples refer to the Focke-Wulf Ente 
and the Stefanutti Anitra S.S.1); the fin effectiveness 
depends entirely on the fin lift which, by reason of the 
comparatively short tail arm and because of the ap- 
preciable magnitude of the sidewash—especially in the 
case of fins placed above the wing—falls short as an ef- 
ficient directional stabilizer. 

Other more recent designs utilize the less efficient 
types of wing-tip fins at a considerable drag penalty. 

The adoption of fins with cambered airfoils facing in- 
wardly that exhibit the aforedescribed high-drag char- 
acteristics at negative fin lift will both improve greatly 
the static directional-stability and rotary damping in 
yaw of tail-first airplanes and minimize the fin drag. 


CONCLUSIONS 


It is revealed that the newly disclosed inwardly fac- 
ing arrangement of fins with certain cambered airfoils 
that afford a sharp profile-drag increase at negative 
angles of attack immediately below the angle of zero 
lift permits the attainment of satisfactory values of 
static directional stability, rotary damping in yaw, and 
side-force slope in tailless airplane designs with straight 
or sweptback high-aspect-ratio wings. In tail-first air- 
planes, it is maintained that the central vertical fin is 
preferable for tailless designs with swept-forward wings 
or wings with very low aspect ratio. 
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Letters to the Editor 


Dear Sir: 


The paper ‘“‘Notes on the Design of Pressurized Cabin Air- 
craft’? by Daniel O. Dommasch, which appeared in the June, 
1946, issue of the JouRNAL, contains a number of inaccuracies 
which call for comment. 

(1) The arbitrary formula, Eq. (1), for reduction of load fac- 
tor with altitude, which is based on the assumption that the load 
factor decreases as the air density, is not consistent with rational 
design criteria. Pressurized transport aircraft will, in general, 
have supercharged engines having high critical altitudes so that 
design speeds will increase with altitude, partially offsetting the 
decrease in density in the load factor relations. Pressurized 
military aircraft may develop the limit maneuvering load factors 
in dive pullouts at altitudes considerably above sea level. The 
argument presented in the paper concerning limitation of load 
factor through loss of controllability with altitude does not 
conform either to theory or to observed experimental data on 
load factors. These factors will considerably alter the conclu- 
sions arrived at in the paper regarding the critical structural alti- 
tude. 

(2) Throughout the paper, arguments are presented against 
the assumption that pressure loads are carried by hoop tension in 
the skin. These arguments are summarized in the abstract by 
the statement: ‘In general, it is not possible to assume that 
diametral pressure loads are carried by hoop tension in the skin 
unless a careful investigation proves that there is no other struc- 
ture capable of carrying these loads (this holds for fuselages of 
any cross section).”” The author seems to favor the assumption 
that the frames carry most of the hoop tension loads: Ka&rman 
and Biot! show that, even for infinitely rigid frames, the frame 
spacing must be of the order 2+/rt/~/3 in order to materially re- 
duce the skin hoop tension stresses. For large fuselages with 
thin skins such as are common in pressurized cabins, even con- 
sidering the effective ¢ to be that corresponding to the increased 
flexural stiffness of the skin due to stringers, this would require 
4- or 5-in. frame spacings, as against the 18- and 20-in. spacings 
now common. Howland and Beed? have shown that the skin 
hoop tension stress is somewhat higher than the average stress 


wae Na 
- (Ap + dt) 


where d is the frame spacing and A; is the frame cross-sectional 
area. Further, Howland and Beed? have indicated a method of 
calculation which leads to results for skin hoop tension in excellent 
agreement with their experimental data. 

(3) Table 1 and Fig. 2 of the paper for the stress coefficients 
of a rectangular membrane appear to be in error. The writer, 
Komm, and Gomza‘ have extended the calculations of Timo- 





shenko’ to a rectangular membrane and have obtained the 
following values: 


b/a o2/W p*Ea*/# oy/W p*Ea?/# 
1 0.252 0.252 
1.5 0.280 0.192 
2 0.374 0.153 
4 0.402 0.117 


It will be seen that these coefficients are considerably lower 
than those given by Dommasch; for instance, at b/a = 1 the 
error is about 75 per cent. These coefficients have been checked 
against calculations made by Féppl using a more accurate method 
and presented later by Sechler and Dunn.’ At b/a = 1, the 
agreement is excellent; as b/a increases agreement is not so good, 
since the representation of the normal deflection along the long 
axis of the plate by a single sine term is not accurate for long 
plates. In all cases, Féppl’s stress coefficients are lower, which 
would further increase Dommasch’s error. 

(4) It is assumed in the paper that the internal pressure does 
not alter the buckling load of the skin in shear and compression, 
and it is stated that ‘the limited amount of test data available to 
date seems to indicate that the buckling load is not reduced by 
internal pressure.”” This is not ai adequate statement of the 
facts. There is available in the literature sufficient evidence, 
both theoretical and experimental, to show that the tensile stress 
produced in the skin by internal pressure exerts a stabilizing ef- 
fect against buckling in either shear or compression. It is only 
necessary to cite several cases discussed by Timoshenko’ covering 
both theoretical and experimental investigations. Actually, the 
problem is greatly complicated in design practice because of 
modification of the tension field postbuckling state of the skin 
by the tension due to pressure. Aside from the effects of increase 
in buckling strength previously mentioned, the tension due to 
pressure tends to stabilize the stiffeners and frames against buck- 
ling under tension field loads. 

(5) In the treatment of the effect of pressurization on skin 
under tension due to external and inertia loads, it is assumed that 
the tensile stress due to bending under external loads, the longi- 
tudinal membrane stress due to pressure, and the longitudinal 
tensile stress due to pressure are all directly additive, as in Eq. 
(7) of the paper. This is far from the true state of affairs. Ac- 
tually, because of the increase in effective skin and the corre- 
sponding shift in neutral axis, the skin tensile stresses due to 
bending under external and inertia loads will be reduced. This is 
clearly shown by the experimental results of Howland and Beed?. 
Further, it can readily be shown‘ that the presence of tension in 
the skin due to the other loads will appreciably reduce the mem- 
brane stresses due to pressure. 


(Continued on page 555) 
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On the Steady Flow of a Gas Through a Tube 
with Heat Exchange or:Chemical Reaction’ 


PAUL CHAMBRE? anp CHIA-CHIAO LINt?t 
California Institute of Technology 


ABSTRACT 


This paper is concerned with a simple phenomenological dis- 
cussion of gaseous combustion through a tube and related phe- 
nomena involving hydrodynamic, thermodynamic, and chemical 
considerations. It is also indicated how the discussion can be 
conveniently developed into a treatment from the point of view 
of microscopic chemical kinetics. Indeed, some interesting 
results are already obtained without such an extension. For 
instance, when there is no chemical reaction, it is found that for 
the range of Mach Numbers of flow 1/>/7 < M < 1 in the 
subsonic region, the temperature of the gas is actually decreasing 
while external heat is being added. It is also shown that, by 
continuous heating, both subsonic flows and supersonic flows 
will eventually approach the sonic state. When the system 
undergoes chemical or physiochemical transformations, the 
situation is much more complicated. In extreme cases, it is 
possible that the temperature of the mixture decreases continu- 
ously for all Mach Numbers while heat is being evolved by the 
transformation. The results apply successfully to the phenome- 
non of detonation. A general discussion of the propagation 
of flame front in a tube is then made, both in the region of slow 
burning and in the region of detonation. 


1. INTRODUCTION 


ij ee PHENOMENON of gas flow through a tube in- 
volving heat exchange has numerous applications. 
Among these may be mentioned the heat exchanger 
and the propagation of combustion in a tube. The 
theory of these phenomena has been studied in some 
detail from various points of view—some stressing the 
thermodynamical aspects,’: ? others the chemical- 
kinetic aspects.4~* 

Certainly a complete understanding of the process 
of combustion and flame propagation can be accom- 
plished only by considering both the chemical process 
and the mechanical process involved. But a thorough 
analysis of this type, involving the detailed molecular 
processes, might obscure the general characteristics of 
the problem. It is, of course, also possible to obtain 
certain results through simple considerations of the 
initial and the final states.. This method has been 
used with some success in dealing with the problem of 
detonation. However, this type of analysis can only 
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furnish incomplete information and does not directly 
pave the way to a more detailed analysis. 

In the present paper a simple phenomenological 
treatment of the problem is discussed which, however, 
leads naturally to the chemical-kinetical treatment as 
the next step. By developing the equations for steady 
flow of a gas in differential forms, it is found possible 
to bring out several essential characteristics of the 
problem not easily obtainable from other treatments. 
On the other hand, the theory as developed here does 
not depend on specific processes and simplifying as- 
sumptions, which would be necessary in a detailed 
treatment including the molecular processes. 

The theory is found to be successful when applied 
to the phenomenon of detonation. The reason 
probably is that detonation does not depend upon the 
rate of chemical reaction so long as that rate is suffi- 
ciently large. Thus, a phenomenological treatment 
would be adequate for such problems if the transfer 
of thermal and mechanical energies is adequately con- 
sidered. In Section 2 the general equations are de- 
veloped in their differential forms, which will be applied 
to the case of a tube of uniform cross section (Sections 
3and 4). When there is no chemical or physiochemical 
transformation, it is shown that by continuous heating, 
both subsonic and supersonic flows will eventually 
approach the sonic state. It is also found that for a 
certain range of Mach Numbers of flow in the subsonic 
region, the temperature of the gas is actually decreasing 
while heat is being added. When the system under- 
goes a themical or physiochemical transformation, 
(Section 4), the situation is much more complicated. 
In Section 5, the propagation of flame front in a tube 
is discussed, both in the region of slow burning and in 
the region of detonation. 


2. GENERAL EQUATIONS FOR STEADY FLOW oF A GAS 
THROUGH A TUBE 


Consider the steady flow of a gas (in general, a 
gaseous mixture) through a tube with slowly varying 
cross section A(x), where x is the distance measured 
along the axis of the tube. Let us assume the gas to 
be an ideal gas, with pressure p, density p, and tem- 
perature T satisfying the relation 


pb = pRT (1) 
Here, R is the gas constant and is given by 


537 








538 OURNAL OF THE AERONAUTICAL SCIENCES—OCTOBER, 1946 


R = nRy (2) 


where m is the number of moles of the gaseous mixture 
per unit mass, and R) is the gas constant per mole 


Ro 


8.3136 X 10’ ergs per °C. per mole 
1.9864 calories per °C. per mole 
1.9864 B.t.u. per °F. per Ib.-mole 

The amount of mass flow per second across any 
section of the tube is a constant for steady flow. Thus, 
if u is the (average) velocity of flow at a given cross 
section, we have 


puA = m, a constant (3) 
The momentum equation is 
pudu = —dp (4) 


Finally, the energy equation for unit mass may be 
written in the differential form 


6Q = cdT — (p/p")dp (5) 


where 6Q is the heat introduced into the system corre- 
sponding to a density change of dp and a temperature 
change dT, and ¢, is the specific heat of the gaseous 
mixture at constant volume. The heat change 6Q 
may be divided into two parts—a part 6Q,, through 
heat conduction across tube wall, and another part, 
5Q2, related to chemical reaction. Thus, 


6Q = 60, + 5Q2, 502 = Qoin (6) 


where 6m is the increase of number of moles per unit 
mass, and Q is the heat /Jiberated by the chemical 
reaction per mole increase of the substance, for constant 
volume reactions at local values of p and 7. For 
example, for the combustion of hydrogen, we have 


He a 1/50. ap H:0 “fb 57.826 kg. cal. 


at 18°C. and 1 atmospheric pressure. There is a 
decrease of half a mole of the reaction mixture in the 
formation of one mole of water for each 57.826 kg. cal. 
of heat liberated. Consequently, the heat of reaction 
at constant volume is 


57.826 — 1/2RoT = 57.527 kg. cal, 


and hence, 


Qo = —115.074 kg. cal. per mole 


In case the chemical reaction does not involve a change 
of the number of moles, we may take the formal pro- 
cedure Qo — © in the following formulas. When there 
is no chemical reaction, Eq. 6 becomes 5Q = 6Q}. 

From the Eqs. (1-6), it can easily be shown that 
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where y, M, H, Ho, and a are, respectively, the ratio 
of specific heats, the Mach Number, the enthalpy per 
unit mass, the enthalpy per mole, and the “velocity of 
sound’’* of the fluid. They are defined by the following 
equations: 


y¥=6,/¢ =1+ R/c, (14) 
M = u/a (15) 
Hy, = H/n = c,T/n (16) 

a = V(b/p) = VyRT (17) 


Note that all these quantities are local values for the 
gaseous mixture. Eqs. (7-17) form the basis of our 
discussion. 


3. TuBES oF UNIFORM Cross SECTION 


These general equations naturally imply the usual 
discussions of the adiabatic flow through a nozzle of 
variable cross section (Laval nozzle). In the present 
paper, however, we are more concerned with the effects 
of external heating and combustion. In order not to 
complicate our problem unnecessarily, we shall restrict 
our discussions to the case of uniform cross section. 
Thus, we may put 


dA =0 (18) 

The Eqs. (7-13) then reduce to 
a 1 5Q én) 19 
= (1 et (19) 


a | (20) 
n 
60 


a (1% 
y(1 — y(1 — M2) 
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* Actually, in a reacting gaseous mixture, even sound waves 
of infinitesimal amplitudes do not propagate with this veloc- 
ity.” 
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dM* _ tae 6Q 2) 
w-7a—-Mm\VAte (24) 
nN Qo H 


To simplify the problem further, let us first consider 
the case in which the number of moles per unit mass 
also remains fixed, i.e., 


dn = 0 (26) 


This corresponds to the case of no chemical reaction 
or to a chemical reaction without change of the total 
number of moles. We distinguish two cases: (a) sub- 
sonic case, M <1; (b) supersonic case, M > 1. 

(a) Subsonic Case M < 1.—In this case it can easily 
be seen that when heat is being added, u and M? 
increase, p and p decrease, while T and a? increase or 
decrease according as 


M <1/V/y or M> 1/7 (27) 


Of particular interest is the fact that the temperature 
actually decreases when heat is being added if 


/V/y <M <1 (28) 


Physically, this means that the kinetic energy of mean 
flow is gained from the heat supplied’ and also from 
the kinetic energy of molecular agitation. To be sure, 
the transfer of heat energy into mechanical energy of 
mean flow is not in contradiction to the second law of 
thermodynamics, because there is a large amount of 
mechanical work involved in such a change (cf. Eq. (5)) 
and the gas is not going through a cycle. 

Thus, as heat is continually being added, the tem- 
perature goes to a maximum when M = 1/,/y and 
then starts to decrease. However, the Mach Number 
M, and the velocity of flow u, increase all the time 
(cf. Eqs. (19 and 24)). For 1/\/y < M <1, the 
velocity of sound decreases, while the velocity of flow 
keeps on increasing. It is therefore expected that 
these velocities will soon become equal. That this 
tendency exists may be seen more clearly from Eq. 
(24). If we consider a supply of heat g per unit length 
of the tube, we have 


1 dMt_1tyMg 
Mt? dx 1-M4H 
For M near unity, and with g ¥ 0, wé have 
2(dM/dx) ~ const./(1 — M) 


which indicates that M approaches unity with an 
infinite rate. The solution is 


(1 — M)? = const. (x) — x) (29) 


where xo represents the point where M = 1. For 
x > xo, the above equation loses its physical meaning. 
Thus, we would expect either that the gas is unable to 
absorb heat (¢ = 0), or rather that the assumed condi- 
tions cease to exist (e.g., combustion stops). 


We note that when g + 0, we have an infinite gradient 
in all the variables u, p, p, T, etc. This corresponds 
to a sharp flame front if the combustion stops beyond 
the point where M = 1. Actually, the effects of 
molecular transfer—molecular diffusion, viscosity, and 
heat conductivity—and of turbulent mixing must 
play a dominant role at such a sharp front; and the 
general equations developed here need modification. 
The study of these transfer effects requires a separate © 
discussion. 

It is important to note that when heat is being con- 
tinually added, a subsonic flow cannot become supersonic 
but will always approach the sonic state. 

(b) Supersonic Case M > 1.—The foregoing dis- 
cussion can also be extended to the supersonic case. 
When heat is being added, u, 7, and a? increase, p 
and p decrease, while the Mach Number M decreases. 
But, since the sign of (1 — yM?*)/(1 — M?) remains 
positive throughout this range, the temperature is 
always increasing when heat is being added. 

The increase of the Mach Number during cooling 
might appear strange at first sight. Evidently this 
is due to a larger reduction of the velocity of sound 
than that of the velocity of mean flow. The velocity 
of mean flow and the flux of kinetic energy both de- 
crease through cooling. These conclusions can easily 
be verified by using the Eqs. (19-23). 

The flow will approach the sonic state in a manner 
similar to the conditions discussed above. Instead of 
Eq. (29), we have 


(M — 1)? = const. (xo — x) (30) © 


Thus, when heat ts continually added, a supersonic flow 
cannot become subsonic but will always approach the 
sonic state. Like the subsonic case, there will be a 
sharp front when the sonic state is reached. 

The fact that both subsonic and supersonic flows, 
when heated, will result in a sonic flow confirms the 
stability of the sonic state. Thus, in order to effect 
a change in the nature of the flow, we must use com- 
bined heating and cooling. Starting with the sub- 
sonic state, the gas is heated to the sonic state. If 
this sonic state is now cooled, the flow may become 
either subsonic or supersonic. Just which state will 
occur depends, probably, on the boundary conditions 
downstream. This is similar to the condition at the 
throat section of a nozzle. The pressure at the exit 
will determine whether the flow becomes subsonic or 
supersonic after the gas passes through the throat 
section. 


4. TuBES OF UNIFORM CROSS SECTION, INVOLVING 
CHANGE OF MOLAR NUMBER OF GASEOUS PHASE DUE 
TO CHEMICAL OR PHYSIOCHEMICAL CHANGES 


We shall now consider the more general case in 
which there is a resultant change of the number of 
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moles per unit mass during the chemical or physio- 
chemical changes considered. 

The chief difference of the present case from the 
simpler case considered in the preceding section lies 
in the presence of the terms involving the dimension- 
less coefficient 


B = Hh/Qo (31) 


This coefficient measures the relative importance of 
the change in the number of moles which occurs in the 
reaction. Since we are now primarily interested in 
the heat evolved from chemical reaction, we shall 
express our equations in terms of 6Q, and 6Q, instead 
of 6Q and dn. Thus, the Eqs. (19-25) become, with 
the help of Eq. (6), 
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In order not to obscure the essential features of the 
effect of molar change, we shall consider the adiabatic 
case in greater detail and discuss the effect of external 
heat exchange briefly at the end of this section. 

As in the preceding section, we shall put more em- 
phasis on the change of temperature and Mach Number. 
We have, with 6Q, = 0, 


dT 1 8Q 


wa ignliines 1-3) 39 
z. - 7 nd 6 Y (39) 


dM? 1+ 7M? 2 *) 
EB IB 40 
M? 1-—-M H dis OY (40) 


It should be noted that 8 and y are not constants, but 
vary during the process of chemical reaction. 

From Eq. (39) it can be seen that the extremum 
of temperature occurs at 


M? = (1/y)[1 — By — 1)/7] (41) 
We see that the extremum of temperature occurs at 
M = 1/Vy 


only when 8 = 0. But it will still be a maximum in 
the subsonic range if 
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—7¥< B< yv/(y — 1) (42) 


holds throughout the reaction process. On the other 
hand, if 8 is always larger than y/(y — 1), there will 
be no maximum in the temperature. The temperature 
is always decreasing while heat is being added. Again, 
when £ lies always in the range —~ < 8 < —y, there 
will be a temperature minimum in the supersonic range. 
If 8 may be equal to y/(y — 1) somewhere during the 
process of reaction, the situation is rather compli- 
‘cated. The detailed discussion will not be included 
here. 

The change in Mach Number with the internal 
addition of heat can also be different from the case 
8 = O discussed in the preceding section. The be- 
havior is similar if 


1+ B/y< 0, or —-y< 8B 


But when 8 < —y, an internal addition of heat will 
decrease the Mach Number in the subsonic range and 
increase it in the supersonic range. Thus, the flow 
will never tend to the sonic state by addition of heat. 
A rapid approach to the sonic state can only be achieved 
if the reaction absorbs heat. However, the authors 
have so far been unable to find any practical application 
of this case. 

It can thus be seen that the exact value of 8B has 
great significance. Therefore, we shall briefly discuss 
its value in connection with the following two reaction 
mechanisms. 


Chemical Reactions 
If an overall reaction, such as the oxidation of butane 
CyHio + 6'/202 — 4CO, + 5H2O 


is considered, the value of 8 will generally be less than 
0.1, and we find that this situation prevails for nearly 
all other overall reaction mechanisms. However, it 
is known that when members of the paraffin series, 
such as butane, pentane, etc., are subjected to tem- 
peratures of 800°C. and higher, they will break into 
smaller fractions. The scission of a C-C bond in 
the above-noted case of butane prior to complete 
oxidation will assertedly cause the molecules to ‘‘crack”’ 
into one molecule of propane and one of methane. In 
such chain reactions 8 can become greater than 1, and 
even exceed the value y/(y — 1). 


Physiochemical Reactions 


These cases occur generally in reactions involving a 
phase transition, such as the condensation of a vapor 
into a liquid. The condensation shock encountered 
in high-speed wind-tunnel operation is of this nature. 
Here 8 can be estimated with the help of Troutons’ 
rule* which states that the molal heat of vaporization 
of any substance at the boiling point under atmos- 
pheric pressure, divided by the absolute boiling tem- 
perature, is approximately constant (16-28). With a 
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value of C, of the order of 5-10 calories per mole 
per °K., one obtains values for 8 between 0.65 and 0.75. 
Thus the discussion in Section 3 holds in such cases. 


Transfer of Heat Across Tube Wall (5Q, ~ 0) 


The effect of transfer of heat across the tube wall 
will alter the conclusions just obtained. The extent 
of this modification will depend on the relative amount 
of heats 6Q,:5Q2. However, the general features 
discussed will probably remain. We shall, therefore, 
not go into the details of this case except to mention 
the following: when all the heat generated by chemical 
reaction is removed across the wall of the tube, 
the temperature may still change as a result of the 
change of number of moles per unit mass. Thus, 
when 





5Q; + 5Q, = 0 
Eq. (35) reduces to 
qT _ _ Biv —1) 8 _ v1) in 
Z, yl — M*) H y(1 — M*) n 


Thus, when there is an increase of the number of moles, 
the temperature increases in the supersonic range but 
decreases in the subsonic range. 


5. COMBUSTION OF GASEOUS MIXTURES THROUGH A 
TUBE 


In this section we shall discuss some aspects of the 
combustion of gaseous mixtures through a tube. It 
will be found that some of the foregoing conclusions 
prove to be significant in such discussions. 

Two types of conditions may be distinguished, 
namely: (1) burning, and (2) detonation. In either 
case, the unburnt gas is at rest (Fig. 1). Consider 
the simplest case where a stationary final state is 
reached, i.e., where the flame front (separating the 
burnt and unburnt gases) is propagating with a uni- 
form velocity D. In the case of burning, the burnt 
gas is flowing in a direction opposite to that of flame 
propagation, and the flame velocity D is low (subsonic 
relative to the unburnt gas). These two states are 
respectively represented by the two branches BKD 
and AJC of the Hugoniot curve (Fig. 2) in a p-V 
diagram, which is constructed for a given initial condi- 
tion O and a given amount of heat liberated per unit 
mass of gas burnt. 

To discuss the problem of combustion by the method 
of the present paper and, in particular, to study the 
nature of the Hugoniot, curve, we consider an observer 
moving with the flame front so that the whole problem 
becomes a steady one, with gas flowing from right to 
left. 

The Hugoniot curve may then be derived from Eqs. 
(1-5). Let us denote the unburnt state by subscript 
( ), and the burnt state by subscript ( )2 For uni- 
form cross section, we have 








unburnt 
(at rest) 


burnt . om 
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() -——> flame front 








. 
Fic. 1. Diagram showing conditions in a tube of burning 


mixture with directions of flow of the burnt gas indicated in the 
cases of (a) slow burning and (b) detonation. 


pi = pa RiTy = mprRoTi (43) 
pe = peR2T2 = meop,RoT>? (44) 

Pilly = polls (45) 
pi + pitts? = pz + prtte? (46) 


To apply Eq. (5), we use the result (Kirchoff’s law) 


2 2b 
—(E; — Bi) = f (8Q — edT) = — f ae (a7) 


for any path. Thus, it is most convenient to take 
Si” (6Q — edT) = Qr — G(T: — Ti) 


where Q, is the heat of reaction liberated at constant 
volume at temperature 7}, and ¢,, is the average specific 
heat of the resultant gaseous mixture (with composition 
corresponding to state 2) between the temperatures 
T; and T>. In calculating “‘the specific heat ¢,, of the 
mixture,’ the composition is assumed to be unchanged, 
and the contributions of the individual components are 
added together. The evaluation of the second integral 
in Eq. (47) may be carried out along the actual path 
satisfying Eqs. (45) and (46). This gives rise to 


2 
f P dp = » (ua? — 7) + a (48) 
&. 9 2 pP2 Ppl 


From these relations it is known that the following 
properties may be ascribed to different portions of the 
Hugoniot curve: 
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Fic. 2. The Hugoniot diagram. 
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Branch AJC AJ uz < a (impossible) 
D= um > aq Y U2 = de 
(Detonation) jc U2 > dz 

Branch AB U2 imaginary 
Branch BKD BK U2 < de 

D= uw < ad K uz = de 

(Slow burning) KD Uz > az (impossible) 





The present results exclude the possibility of the 
branches AJ and KD during a continuous exothermic 
-chemical reaction. Indeed, we would expect the points 
J and K (where M = 1) to be the only conditions 
actually occurring. This is found to be the case with 
‘detonation. The wave front is found to have an 
extremely sharp front in accordance with theoretical 
prediction, and the detonation speed calculated from 
the conditions at J agrees satisfactorily with experi- 
mental results, provided the end products are known.* 

With the case of slow burning, the present theory is 
perhaps oversimplified. Here, the effects of molecular 
transfer and turbulent mixing come in, and the flame 
velocity is usually much lower than that predicted 
from conditions at K, which may, therefore, be re- 
garded as an ideal upper limit. 

The actual conditions of burning in tubes of various 
‘sizes may be described as follows. For an extremely 
small tube, the flame propagates slowly, corresponding 
to a condition along BK close to B. This is the region 
-of laminar flow. For larger tubes, the laminar flow is 
unstable’ and the flow becomes turbulent. The flame 
speed is suddenly increased at the transition. In 
the turbulent regime, the flame speed first increases 
approximately in proportion to the tube diameter,’ 

and then approaches the conditions at K as a limit. 
Landau’ states that if the flame velocity predicted by 
this proportionality law exceeds that given by K, only 
-detonation can occur. It seems conceivable, however, 
that the flame speed may depart from the linear law 
and may even finally take the fixed speed characterized 
by K, analogous to the constant detonation speed at J. 
} The instability theory of Landau neglects the effect 
of viscosity, diffusion, and heat’ conductivity. It is 
interesting to consider whether or not the limit of 
transition (say, the critical Reynolds’ Number) could 
be determined from a theory including those quantities. 


CONCLUSIONS 


This paper is a rudimentary discussion of the com- 
bustion problem from the hydrodynamic and thermo- 


e 
* The conclusion that the gas leaves the detonation front at 


sonic speed has also been obtained from stability arguments 
(see, e.g., page’351 of Becker’s paper’). 





dynamic points of view. It is certainly only the be- 


ginning of a more extensive investigation but already 


it has brought out rather interesting results. For 
example, it must now be kept in mind that the addition 
of heat to a flow system does not necessarily lead to an 
increase of temperature. Such results are somewhat 
surprising if one is not aware of the inertia forces that 
must enter the picture of rapid motion of fluids. There- 
fore, the treatment presented, although rudimentary, 
does indicate the necessity for bringing into considera- 
tion the hydrodynamics of the problem in conjunction 
with the usual thermodynamic and chemical points of 
view. 

The next step of the investigation is perhaps a study 
of the microscopic mechanism. Thus, it might be of 
interest to consider the molecular processes by which 
both the external heat energy supplied and the molecu- 
lar energy of the gas are transformed into the me- 
chanical energy of flow under certain conditions. In 
carrying out the investigation of the chemical-kinetical 
aspects, the equations developed here will all be re- 
tained, and Eq. (6) is the only equation in which chemi- 
cal reaction must be considered. Of course, the actual 
calculation of a given problem still requires the usual 
step-by-step process involved in the determination of 
chemical equilibrium. But as soon as the chemistry 
is known for all the states that may be involved in the 
problem, all the necessary equations will be furnished 
by the Eqs. (1-4), (6), and one of the differential 
equations—Eqs. (5), (7-12). 
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Data to Standard Air Conditions by 
the Temperature Altitude Method 


J. G. BEERER* 


North American Aviation, Inc. 


SUMMARY 


Analysis of flight performance results on modern high-speed 
airplanes has shown that currently used methods for the correc- 
tion of observed performance to standard air conditions are in- 
adequate. Present methods have been derived on the assump- 
tion that the compressibility of air is negligible, and therefore 
unsatisfactory results have been obtained in applying such 
methods to airplanes flying at Mach Numbers of approximately 
0.6 and above. 

A new method for the reduction of observed flight data is sug- 
gested herein, in which allowance is made for Mach Number ef- 
fects on the aerodynamic characteristics of the airplane and pro- 
peller and on the thermodynamic characteristics of the engine. 
This new method is called the “‘temperature altitude’”’ method, 
since the observed performance is corrected for the deviation in 
atmospheric pressure from the standard pressure corresponding 
to the observed air temperature. 

The method is applicable to airplanes driven by Otto cycle 
engines and propellers, by gas-turbine engines and propellers, 
by gas-turbine jet propulsion engines, and by combinations of 
these types of engines. Although intended principally for high- 
speed airplanes, the method is applicable to airplanes of all speed 
classes from low Mach Numbers up to and including supersonic 
flight velocities. 

The analysis and proof of the temperature altitude method are 
somewhat lengthy, but the routine application of the method to 
the reduction of observed flight data is of extreme simplicity. 
Use of the method should result in improved correlation of high- 


’ speed flight performance data. 


INTRODUCTION 


i er DETERMINATION of the performance of an air- 
plane in standard air conditions from the results 
of tests made under nonstandard atmospheric conditions 
has been the subject of many papers. Among the more 
recent works are those of Jones,’ Allen,* Cameron,’ and 
Reed.‘. Every method advanced to date for the cor- 
rection of observed performance to standard has been 
based on the assumption that airplane and propeller 
characteristics are not affected by the compressibility 
of air. As a result, current correction methods give 
unsatisfactory results when used to establish the stand- 
ard air performance of modern high-speed airplanes. 
The advent of the gas-turbine jet-propelled airplane 
has further emphasized the need for a more suitable 
method. 


Presented at the Flight Test Session, Fourteenth Annual Meet- 
ing, I.A.S., New York, January 29-31, 1946. 
* Assistant Chief Technical Engineer. 


In this paper a method of performance correction is 
suggested which allows for the effects of Mach Number 
on the airplane, propeller, and engine characteristics. 
This method is called the ‘‘temperature altitude’”’ 
method because it corrects the airplane performance at 
the observed ambient air temperature for the variation 
of air pressure from the standard value corresponding 
to the observed temperature. Therefore, the basic 
variable governing compressibility effects on the air- 
plane, the speed of sound in air, is held constant be- 
tween observed and standard air conditions. 


NOTATION 

C = airplane true rate of climb, ft. per sec. (ft. per min.) 

Cp = airplane drag coefficient 

Cr = airplane lift coefficient 

Cp = propeller power coefficient 

D = airplane drag, lbs. 

D = propeller diameter, ft. 

F = net thrust of jet, lbs. 

F; = gross thrust of jet, Ibs. 

g = gravitational constant, 32.2 

Hr = temperature altitude, ft. 

Hp = pressure altitude, ft. 

N = engine r.p.s. (r.p.m.) 

n = propeller r.p.s. 

Ne = supercharger r.p.s. 

Pe = atmospheric pressure, absolute, Ibs. per sq.ft. (in. 
Hg) 

Pa = manifold pressure, absolute, lbs. per sq.ft. (in. Hg) 

Pe = supercharger, or compressor, inlet impact pressure, 
absolute, lbs. per sq.ft. (in. Hg) 

R = gas constant for air, 53.3 

Ss = airplane wing area, sq.ft. 

T,(t,) = atmospheric temperature, °R (°F.) 

7; = supercharger, or compressor, inlet impact tempera- 
ture, °R 

Tm(tm) = manifold temperature, °R (°F.) 

= airplane velocity, ft. per sec. (m.p.h.) 

Vay = airplane velocity corrected to standard air condi- 
tions and to reference gross weight, ft. per sec. 
(m.p.h.) 

w = engine charge flow, lbs. per sec. 

We = engine airflow, lbs. per sec. 

wy = fuel flow, lbs. per sec. 

W = airplane weight, Ibs. 

p = air density, slugs per cu.ft. 

n = propeller efficiency 

BHP = brake hp. 

CHP = supercharger hp. 

FHP = rubbing friction hp. 

IHP'- = indicated hp. 
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JHP = exhaust jet hp. 
PHP = pumping hp. 
THP = thrust hp. 


Subscript ‘‘s’”’ refers to standard air conditions. 

Subscript ‘‘t’’ refers to observed air conditions. 

Dimensions given in parentheses, while not dimensionally con- 
sistent, are generally used in quoting airplane performance data. 


Analysis of the 
Temperature Altitude Method 


I. THE PROPELLER-DRIVEN AIRPLANE WITH SUPER- 
CHARGED OTTo CyYcLE ENGINE 


Engine Performance 


Studies on the subject of performance correction 
usually start with an investigation of the effect on en- 
gine performance of changes in atmospheric conditions. 
For the purposes of this study the variation of engine 
characteristics with change of ambient air pressure 
must be known at a fixed value of ambient air tempera- 
ture. Few published data exist on this problem; most 
engine testing has been done to check the effect of the 
standard altitude on engine characteristics and to check 
the effect of variation of ambient air temperature at 
constant ambient pressure. 

The effect of ambient air pressure at constant tem- 
perature on the volumetric efficiency of several Otto 
cycle engines was thoroughly studied by Zeyns.5 The 
tests, made with unsupercharged and supercharged 
engines, showed that the volumetric efficiency of a 
given engine at a fixed r.p.m. varied only with the ratio 
P/P a. This result is in agreement with theory.® 
Therefore, when operating at constant r.p.m. with con- 
stant manifold charge temperature T7;,,, and with con- 
stant supercharger pressure ratio P,,/P,, the weight 
flow of air drawn through an Otto cycle engine will vary 
directly as the absolute ambient air pressure. 

The performance characteristics of centrifugal super- 
chargers were studied theoretically and experimentally 
by Capon and Brooke.’ The dimensional analysis and 
test data presented show that the temperature ratio 
T/T, and the pressure ratio P,,/P,, of a given centrif- 
ugal compressor are functions of the two parameters, 
n.//T, and w./T,/P,. It is apparent, for the case of 
an engine with gear-driven supercharger, that if the 
engine is operating at constant r.p.m., with constant 
ambient air temperature, the compressor temperature 
ratio and pressure ratio will be constant only if the 
weight flow of air varies directly with the ambient air 
pressure. ; 

It may be concluded, therefore, that the weight flow 
of air through an Otto cycle engine equipped with a 
gear-driven supercharger without charge cooling will 
vary directly with the ambient air pressure when 
r.p.m. and ambient air temperature are held constant. 
Further, the manifold temperature will be constant, and 


the manifold pressure will vary directly as the ambient 
air pressure. 

Theoretical considerations show that if the manifold 
temperature, the fuel-air ratio, and the pressure ratio 
P,,/P, are all fixed quantities for a given engine, then 
the instantaneous pressure ratios and temperature 
ratios are also fixed throughout the entire thermody- 
namic cycle. Therefore, the mean effective pressure 
acting on the pistons must vary directly with the am- 
bient air pressure. Or at constant r.p.m. and ambient 
temperature, the indicated horsepower of a given engine 
should be directly proportional to the ambient air pres- 
sure. Test data on this point are scarce; data on obso- 
lete unsupercharged engines shows that IHP does vary 
directly with ambient air pressure within limits of ex- 
perimental accuracy.*: ® 


The power output of an Otto cycle engine is usually 
divided thus: 


IHP = BHP + CHP + PHP + FHP (1) 


It was shown above that IHP varied as the ambient air 
pressure with constant ambient temperature and r.p.m. 
The power absorbed by the gear-driven supercharger 
CHP at constant r.p.m. and ambient temperature varies 
directly as the weight flow of charge air;’? but since the 
charge airflow varies directly with ambient pressure, 
then CHP must vary directly as ambient air pressure. 
The pumping horsepower term PHP is due to the pres- 
sure differentials acting on the pistons during intake 
and exhaust strokes. It was shown above that pres- 
sure ratios are fixed at any point in the engine cycle 
with fixed ambient temperature, therefore the mean 
effective pressure must vary directly with ambient pres- 
sure, and hence at fixed r.p.m. the pumping power term 
must be proportional to ambient pressure. The fric- 
tion power term FHP represents the power absorbed 
in the sliding friction of the pistons, the bearing losses, 
and the accessory drive power. It is often assumed 
that FHP for a given engine varies only with r.p.m. 
but tests have shown that FHP is appreciably affected 
also by the average gas pressure on the pistons. The 
extensive tests reported by Sparrow and Thorne” 
showed that FHP drops with increasing cylinder wall 
temperatures (cooling water). Extrapolation of these 
data to the range of temperatures encountered on mod- 
ern engines shows that the FHP term may be extremely 
small. 


Therefore, at constant r.p.m. and constant ambient 
temperature, IHP, CHP, and PHP vary directly as the 
ambient air pressure, and, since FHP is likely to be 
small, the remaining term of Eq. (1), the brake horse- 
power, should vary nearly as the ambient pressure. 
Few experimental data are available on modern en- 
gines to test this conclusion. Ragazzi reported the 
results of altitude chamber tests conducted on a Fiat 
A80 18-cylinder two-row radial air-cooled engine of 
2,790-cu.in. displacement, equipped with a single-stage 
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Fic. 1. Variation of brake horsepower with ambient air pres- 
sure. Fiat A-80 supercharged radial airtcooled engine tested in 
altitude chamber with fixed ambient air temperature of 15°C., 
2,100 r.p.m., and fixed fuel air ratio. Points shown are taken 
from faired experimental data.™ . 


gear-driven centrifugal compressor." A part of these 
data has been cross-plotted in Fig. 1; it shows the vari- 
ation of brake horsepower at 2,100 r.p.m. with ambient 
air pressure, P, for constant values of P,/P,. Am- 
bient air temperature was held at 15°C. and fuel-air 
ratio was held constant. The normal boost ratio, P,,/ 
P,, for the A80 supercharger at 2,100 r.p.m. is 1.56 with 
15°C. inlet air, therefore, all values of P,,/P, less than 
1.56 represent throttled conditions, while higher values 
correspond to rammed full-throttle operation. Data 
on the throttle positions for the various boost ratios 
were not reported, but from dimensional analysis’ a 
fixed throttle setting would be expected to correspond 
to each of the P,,/P, lines plotted. 


It is apparent from Fig. 1 that this representative 
modern geared supercharged engine shows a BHP 
characteristic which varies practically directly with 
ambient air pressure when r.p.m., ambient air, and 
throttle. are held constant. This indicates that the 
friction horsepower is either extremely small or that it 
varies as the mean effective cylinder pressure (i.e., 
with the ambient pressure). 


The preceding discussion dealt with the case of an 
engine under static atmospheric conditions, that is, with 
zero flight speed, so that the compressor inlet pressure, 
P., could be considered equal to ambient pressure, P,, 
and T, equal to 7,. In high-speed airplanes the air 
compression in the engine air inlet duct is of major im- 
portance. From considerations of the airflow cycles 
involved it can be shown that the duct pressure ratio, 


P./P,, will be constant if the flight and duct flow Mach 
Numbers are fixed. Thus, for fixed ambient air tem- 
perature, P./P, will be constant if the flight speeds are 
equal and if the weight flow of air varies as the ambient 
pressure. The compressor inlet temperature, 7., will 
be constant with constant ambient air temperature, 
T,, only if the true flight velocity is constant, be- 
cause the temperature rise 7-7, can be considered as 
due only to the adiabatic compression in the 
duct. 


The above analyses considered a gear-driven super- 
charged engine without charge cooling. Where charge 
cooling is involved, either during compression or after 
compression, additional considerations are needed. 
Charge cooling systems are of two types—heat transfer 
from charge to cooling air, or charge to liquid to cooling 
air. In both cases the primary heat transfer factors 
are the charge flow and the cooling airflow. With a 
given cooling air duct configuration, the cooling airflow 
can be shown to vary directly as the ambient air pres- 
sure for the case of an airplane flying at constant Mach 
Number and fixed angle of attack, because for these 
flight conditions local pressure ratios about the airplane 
will be constant, and hence also the pressure ratio caus- 
ing flow through the duct. Since, for a given cooling 
system, the rate of heat transfer varies somewhat less 
rapidly than the air and charge flows, the cooling ef- 
fectiveness will change slightly as the flows vary. How- 
ever, for practical purposes it may be assumed that the 
effectiveness will be unchanged, and therefore that the 
charge temperature drop through any heat exchanger 
will be constant for variations ‘in air and charge flow 
with constant cooling air and charge air temperatures 
entering the heat exchanger. 


Exhaust Jet Performance 


An appreciable part of the propulsive effort obtained 
from a supercharged Otto cycle engine in a high-speed 
airplane is derived from the thrust generated by the 
single-cylinder jet exhaust stacks usually found on this 
type installation. It was shown in the section on en- 
gine performance that with fixed ambient air tempera- 
ture, the local pressure ratios and temperature ratios 
in the engine cycle are constant. The theory of the flow 
of gases through nozzles and orifices shows that the 
weight flow of gas through a given nozzle, operating at 
a fixed pressure ratio with fixed inlet temperature, will 
vary directly as the end pressure (i.e., the ambient air 
pressure). This is entirely consistent with the experi- 
mental results of Zeyns on the variation of volumetric 
efficiency.’ The exhaust velocity from the jet stacks, 
operating under the above conditions, will be constant, 
and therefore the thrust reaction wit vary directly as 
the ambient air pressure. With constant flight speed 
the thrust horsepower generated by the jet stacks, JHP, 
will also vary directly with the absolute ambient pres- 
sure. 
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Propeller Performance 


The aerodynamic characteristics of a given propeller 
are fixed by the following operating conditions: ad- 
vance ratio V/nD, the flight Mach Number, and the 
power coefficient c,.'7 The power coefficient is defined 
as 


55 550gRT, (BHP) 
> n3D>—(P,) 





(2) 


where p has been repiaced by P,/gRT,. 

When operating under conditions of fixed r.p.m., 
ambient air temperature, and flight velocity, it is ap- 
parent that the aerodynamic characteristics, including 
the propulsive efficiency, of a given propeller will be 
constant if the engine brake horsepower varies directly 
as the ambient air pressure. 


Airplane Performance 


The aerodynamic characteristics of a given airplane 
.of a fixed geometry are determined by the Mach-Num- 
ber of flight and by the angle of attack of the airplane.!” 
Reynolds Number effects are assumed negligible. 
Therefore, in flight under conditions of constant am- 
bient temperature and fixed angle of attack, the drag 
coefficient, Cp, and lift coefficient, C,, will be constant 
only if the flight speed, V, is constant. 


Flight Operating Conditions 


It is apparent from the preceding analyses that the 
following flight operating conditions should be held 
fixed in correcting obsérved airplane performance to the 
standard air condition: (a) airplane angle of attack, 
(b) engine r.p.m., (c) engine throttle setting, (d) air- 
plane aerodynamic configuration, and (e) atmospheric 
temperature. 

The following flight conditions may be considered 
variable in correcting from test to standard air condi- 
tions: (f) atmospheric pressure, and (g) airplane gross 
weight. 


The “‘Reduced”’ Altitude 


The flight altitude is not specified directly in the fore- 
going flight operating conditions. However, it is re- 
quired that the observed air temperature be held con- 
stant in correcting to standard conditions while the 
air pressure is allowed to change. This means that the 
standard or ‘‘reduced’’ altitude is defined by the ob- 
served air temperature only. Thus, if the observed air 
temperature were 59°F., the standard or “reduced”’ 
altitude would be sea level, regardless of the observed 
pressure altitude, density altitude, or tape-measure 
height. 


The ‘‘Reduced”’ True Air Speed 


The balance of horizontal forces acting on an airplane 
in stabilized level flight can be written 


Thrust = Drag (3) 
or in the usual power form 
THP = (CpSV*p/2)/550 (4) 


Replacing the density, p, by its equivalent P,/gRT,, and 
writing THP as 7BHP + JHP, Eq. (4) becomes 
CpSV*P, 
BHP + JHP = ——— 5 
wee tN 550RT,2¢ (5) 
Eq. (5) is the general airplane power relationship for 
stabilized level flight under any atmospheric conditions, 
Using ‘‘t’” subscript for observed test conditions and 
“s’’ for standard air conditions, Eq. (5) can be written 
for each condition. Then, dividing the standard by 
the test equation and eliminating constants, 


BHP, + JHP,  Co,PaTa,V,* 








n.BHP, + JHP,  Cp,PaTaVe 6) 
A study of Eq. (6) on the basis of the flight operating 
conditions specified shows that it is an equality only for 
the case of 
V. = V, (7) 
because for this condition: 

(a) Co; = Cp,, since standard and test flight Mach 
Numbers are equal and the angles of attack are equal. 

(b) BHP,/P«a, = BHP,/Pa2,, since standard and test 
flight Mach Numbers are equal, engine r.p.m.’s are 
equal, supercharger operating conditions are equivalent, 
and engine throttle setting is unchanged. 

(c) JHP;/Pe, = JHP,/Pa,, since jet stack thrust 
varies as ambient pressure when standard and test air- 
plane velocities are equal. 

(d) ms = my since flight Mach Number, propeller 
advance ratio, and power coefficient for standard air 
are equal to those for test air conditions. 

Therefore, the “reduced” true air speed is equal to 
the observed true air speed. 


The ‘“‘Reduced’’ Gross Weight 


The balance of vertical forces acting on an airplane 
in stabilized level flight can be written 


Weight = Lift (8) 
or in the usual form 

W = C,SV"p/2 (9) 
Eq. (9) can be handled in a manner similar to Eqs. (4) 
and (5), with the following result, 


W, CL; Pas To V,? 
a ee (10) 
W, CL; Pa, Ta; V? 
Analysis of this equation for the specified flight oper- 
ating conditions shows that 
W,/W;, = Pa,/Pa; (11) 


because for this condition: 
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(a) V, = V,, the reduced true speed is equal to ob- 
served. 

(b) Cr, = Cz, since the standard and test flight 
Mach Numbers are equal and the angles of attack are 
equal. 

Therefore, the equivalent or “‘reduced”’ gross weight 
under standard air conditions is given by the test gross 
weight multiplied by the ratio of the standard altitude 
pressure to the observed altitude pressure. 


The ‘‘Reduced’’ Rate of Climb 


The simplified equation for rate of climb universally 
used is 


C= (nBHP — THP) (12) 


where 7BHP is the thrust horsepower available, and 
THP is the thrust horsepower required to balance the 
airplane drag. Eq. (12) may be written for standard 
and test conditions, 

W,.(n,;BHP, — THP,) 


We 
C, W.().BHP, — THP,) (13) 





For the flight operating conditions previously defined, 
it was shown that W,/Pa, = W,/Pa,, n» = 1, BHP,/ 
Psa, = BHP,/Pe,;. And since the equation for thrust 
horsepower required is similar to Eq. (4), it can be 
shown that for the specified flight conditions, THP,/ 
Pa, = THP,/Pa,. Therefore, 


C,= C; (14) 


‘ 


That is, on the temperature altitude basis, the “re- 
duced” rate of climb in standard atmosphere is equal 
to the observed true rate of climb. As in the level 
flight case, the reduced gross weight is equal to the 
product of the test gross weight by the ratio of standard 
to observed altitude pressure. 


The ‘‘Reduced’”’ Manifold Pressure and Temperature 


In the section on engine performance it was shown 
that the manifold pressure would vary as the ambient 
pressure, and the manifold temperature would be con- 
stant, if the following operating conditions were satis- 
fied: (a) constant engine r.p.m., (b) fixed throttle 
setting, (c) constant ambient air temperature, (d) con- 
stant airplane flight Mach Number, and (e) fixed con- 
figuration of any cooling air ducts. 

Since these requirements are all satisfied by the speci- 
fied flight operating conditions, the following manifold 
pressure and temperature relationships will exist, 


Pm,/Pm, = Pa,/Pa, (15) 
tm, = tm, (16) 


Therefore, the “reduced” manifold temperature is 
equal to the test manifold temperature, and the ‘“‘re- 
duced” manifold pressure is equal to the test value 


multiplied by the ratio of the standard to the observed 
altitude pressure. 


Application of Method 


In the foregoing analysis of the temperature altitude 
method the following simple relations were shown to 
exist between standard and observed flight performance 
data: 


V.= V, 

e = W (Pa;/Pa;) 
C, = C, 
Pm; = Pm,(Pa,/Pa,) 


tms = tm, 


Application of the method is perhaps best illustrated 
by correcting an actual airplane flight-test run. The 
following observed flight performance data were ob- 
tained on a North American Aviation experimental 
fighter airplane during a run at 3,000 r.p.m., full throt- 
tle, level flight: 


Observed free air temperature, fa, = —52.5°F. 
Observed pressure altitude, H, = 32,255 ft. 
Observed altitude pressure, Pa, = 8.006 in. Hg 
Observed true air speed, V,; = 464 m.p.h. 
Observed manifold pressure, Pm, = 59.6 in. Hg 
Observed manifold temperature, tm, = 201°F. 
Estimated flight weight, W, = 6,620 Ibs. 


The observed free air temperature of —52.5°F. cor- 
responds to a temperature altitude of 31,300 ft. and a 
standard air pressure of 8.38 in. Hg.'* From the above 
relations, the observed flight conditions can be reduced 
to standard: 


Standard free air temperature, ta, = —52.5°F. 
Standard altitude, Hy; = 31,300 ft. 

Standard altitude pressure, Pa, = 8.38 in. Hg 
Standard true air speed, V, = 464 m.p.h. 
Standard manifold pressure, Pm, = 62.4 in. Hg 
Standard manifold temperature, tm, = 201°F. 
Standard flight weight, W, = 6,940 Ibs. 


It is at once apparent that the method is extremely . 
easy to use. 

The pressure ratio factor, Ps,/P«,, used in correcting 
manifold pressure and gross weight is plotted in Fig. 2. 
against pressure altitude for a range of ambient. air 
temperatures.'* It is obvious that weight corrections 
will be large if air temperatures are far from standard. 
For small temperature variations from standard, a 
rough rule is that the reduced weight. and manifold 
pressure will exceed the test values by one per cent for 
each degree (°F.) the ambient air temperature exceeds 
the standard. 

One of the main objections to the temperature alti- 
tude method lies in the fact that the reduced weights 
may vary widely from the test weights when air temper- 
atures are far from standard. In order that the re- 
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Fic. 2. Correction factor for reduction of test-flight weights 

and pressures to standard air by temperature altitude method. 


duced weight may approximate the reference or design 
weight, the take-off weight of the test airplane should 
be adjusted to suit the anticipated upper air conditions. 
For example, if air temperatures are 20°F. above stand- 
ard at the test pressure altitude, the take-off weight 
should be about 80 per cent of the desired reduced 
weight, plus the fuel needed to climb to the test altitude. 
For many airplanes the large range of take-off weight 
required by this method may be impractical. In cases 
where take-off weight control is not practical, resort 
must be made to a calculated correction for weight 
based on incompressible airflow; this procedure in- 
volves a change in airplane angles of attack which vio- 
lates one of the basic flight operating conditions. In 
cases where the Mach Number effects on airplane drag 
and propeller efficiency are small, the error introduced 
by a calculated weight correction will be a small part 
of the correction itself. But where drag and efficiency 
are varying rapidly with Mach Number, the error may 
approach the value of the correction. 

Another: objection to the temperature altitude 
method arises from the fact that observed air tempera- 
tures will sometimes lie outside the limits of standard 
atmosphere, 59° to —67°F.,'* and thus give fictitious 
reduced altitudes. However, such points may be used 
to advantage in establishing the shapes of performance 


curves. 


Experimental Verification 


Extensive flight trials on one of North American 
Aviation’s experimental fighter-type airplanes have 
provided data with which to check the standard air 
performance obtained by the usual density altitude cor- 


rection method, as compared with the temperature alti- 
tude method. The flight-test airplane had the follow- 
ing physical characteristics: 

Type: Low-wing, single-engine tractor, single- 
place fighter. 

Normal gross weight: 6,900 Ibs. 

Wing: Area, 236 sq.ft.; span, 37.1 ft. 

Airfoil: Low-drag type. 

Engine: Packard-built Rolls-Royce Merlin. 
Water-cooled. Two-stage, two-speed supercharger 
with aftercooler and water injection. 

Three different propellers, each with a diameter of 
11.0 ft., designated A, B, and C, were tried during the 
course of the flight tests. 


The following flight-test instruments were mounted 
in the airplane photorecorder: sensitive air-speed indi- 
cator, sensitive altimeter, tachometer, manifold pres- 
sure, manifold charge temperature, impact air tempera- 
ture, clock; and other special instruments depending 
on the particular test. 

The air-speed indicator and altimeter systems, using 
an underwing pitot tube and an aft fuselage static 
source, were calibrated in flight at sea level in a manner 
identical to that described by Jones.’ A similar system 
in another experimental airplane was flown at 30,000 
ft. in addition to sea level in order to check the Mach 
Number effect on position error; no measurable effect 
was found up to the maximum flight Mach Number of 
0.65. 

The impact air temperature was measured by an 
electrical resistance bulb mounted in the engine air in- 
take duct at the section of maximum area. Calibra- 
tion flight tests made at sea level and at altitude showed 
that this system recorded the true impact temperature 
within 1°F. at all level flight air speeds. 

The ambient air temperature and true flight velocity 
are obtained as follows: The flight Mach Number is 
known from the ratio of impact pressure to static, or 
altitude, pressure. The ratio of impact to static, or 
ambient, air temperature is a function of flight Mach 
Number only, hence the ambient temperature may be 
determined. Then, since the flight Mach Number is 
known, and the speed of sound is determined from the 
ambient air temperature, the true flight velocity is 
easily obtained. 

The deficiencies in the common density altitude 
method for reduction of performance data to standard 
became apparent during early flight trials on this air- 
plane. Table 1 gives full-throttle level-flight data ob- 
tained during Flights 21 and 25; these flights were 
made 3 days apart and with unchanged airplane con- 
figuration. At a given pressure altitude the outside 
air temperatures differed by about 12°F. The data re- 
duced by the density method are plotted in Fig. 3. 

The density altitude method as applied herein is 
similar to the method widely used in Britain’ and is 


‘briefly as follows: The observed true air speed is first 
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TABLE 1 
Experimental Fighter Airplane 
Level-Flight Full-Throttle Performance Data at 3,000 R.P.M. 
Temperature Altitude Reduction Method 





Flight 21 Propeller A. Propeller-gear ratio 0.479. Exhaust stack exit area 2.5 sq.in. Normal gross weight 6,900 Ibs. 
Hy ta, V; Pm W, tm; Hr* Pm, W; Vow t 
Run Ft. "FR; M.p.h. In. Hg Lbs. KS Pa;/Pa; Ft. In. Hg Lbs. M.p.h 
1 34,150 —70.5 454 54.5 6,760 184 0.902 36,300 49.1 6,100 449.5 
2 33,180 —69.0 456 57.2 6,680 185 0.880 , 35,900 50.3 5,870 450.5 
3 32,210 —66.0 459 59.4 6,630 189 0.875 35,050 52.0 5,800 453.5 
4 31,730 —62.5 460 60.7 6,600 192 0.895 34,100 54.4 *5,910 455.5 
Flight 25. Propeller A. Propeller-gear ratio 0.479. Exhaust stack exit area 2.5 sq.in. Normal gross weight 6,900 Ibs. 
Hy ta, Vz Pm W, tm; Hr Pm, W, Vist 
Run Ft. °F. M.p.h. In. Hg Lbs. re, Pa;/Pa, Ft. In. Hg Lbs. M.p.h. 
1 34,290 —60.5 461 54.7 6,760 190 1.036 33,500 56.7 7,010 461.5 
2 33,270 — 57.0 463 56.2 6,670 196 1.032 32,550 . 6,890 463.0 
3 32,260 —52.5 465 59.6 6,620 201 1.045 31,300 62.4 6,920 465.0 
+ 31,730 — 50.5 466 60.9 6,560 203 1.050 30,700 63.9 6,880 466.0 
5 31,310 —49.0 467 61.9 6,510 205 1.050 30,300 65.0 6,830 467.0 





* Runs | and 2, flight 21, were reduced to a ‘‘standard” temperature altitude by assuming that lapse rate continued beyond 35,332 ft. 
t Vs» is the reduced true air speed corrected to the normal weight of 6,900 Ibs. 
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Fic. 3. Experimental fighter airplane. Level flight, full- 
throttle performance at 3,000 r.p.m. Density altitude method. 
Normal gross weight, 6,900 Ibs. 


corrected to a standard or reference weight. The re- 
sultant true air speed is then plotted at the standard 
altitude corresponding to the observed air density. 
The observed manifold pressure is plotted on a “‘half- 
and-half’’ basis, that is, it is plotted at an altitude half- 
way between the observed pressure altitude and ob- 
served density altitude. Variations of the density al- 
titude method in use in this country may differ in detail 
application but the end results are practically identical. 
The reduction method is based on the assumption that 
the airplane and propeller characteristics are unaffected 
by Mach Number. 

Inspection of Fig. 3 shows that the standard speed 
performance as determined by the density method is 
unsatisfactory, indicating therefore that the airplane 
and propeller characteristics are being affected by Mach 
Number. Observed ambient air temperatures are also 
plotted showing that the airplane flies fastest on the 
hottest day. 

Table 1 and Fig. 4 give the performance reduced to 
standard air conditions by the temperature altitude 
method. No attempt was made in this series of tests 


to adjust the take-off weight of the airplane, hence the 
reduced gross weights, shown in Table 1, vary appreci- 
ably between the two flights. The speeds plotted in 
Fig. 4 have been corrected to a normal gross weight of 
6,900 Ibs. To do this it was assumed that the speed in- 
crement due to weight increment could be calculated 
on the basis of incompressible flow. Since this airplane 
is flying at speeds where Mach Number effects are ap- 
preciable, the calculated speed changes will be too 
large. 

The speed correction due to weight increment is made 
as follows: It is assumed that the thrust horsepower 
developed by the propeller and jet exhaust stacks is 
constant between uncorrected and corrected speed, 
and that second order differences can be neglected. 
Then, on the basis of Eq. (4), the velocity increment can 
be written, i'% 


AV./V = —1/;ACp/Cp (17) 


where AC> is the increment in drag coefficient due to 
flying at a changed C,. It is now assumed that the 
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Fic. 4. Experimental fighter airplane. Level flight, full- 
throttle performance at 3,000 r.p.m. Temperature altitude 
method. Normal gross weight, 6,900 Ibs. 








JOURNAL OF THE AERONAUTICAL SCIENCES—OCTOBER, 1946 





TABLE 2 
Experimental Fighter Airplane 
Level-Flight Full-Throttle Performance Data at 3,000 R.P.M. with Water Injection 
Temperature Altitude Reduction Method 





Flight 128 Propeller B. Propeller-gear ratio 0.479. Exhaust stack exit area 2.5sq.in. Normal gross weight 7,200 Ibs. 
H, tay V; Pm, W, Hr Pmg W,; Vag* 
Run Ft. -. M.p.h. In. Hg. Lbs. Pa,/Pa Ft. In. Hg. Lbs. M.p.h. 
1 30,100 —41.5 463 62.1 7,420 1.088 28,200 67.6 8,070 465.5 
2 28,250 —31.0 469 , 67.7 7,350 1.146 25,200 77.5 8,410 471.5 
3 25,810 —18.0 475 75.0 7,250 1.197 21,600 89.8 8,680 477.5 
4 23,830 — 9.0 478 81.0 7,140 1.225 19,100 99.2 8,730 480.0 
5 22,290 — 2.5 480 85.9 7,060 1.232 17,250 106.0 8,700 481.5 
6 21,560 1.0 482 89.1 7,020 1.247 16,250 111.0 8,750 483.5 
Flight 129 Propeller B. Propeller-gear ratio 0.479. Exhaust stack exit area 2.5 sq.in. Normal gross weight 7,200 Ibs. 
H, tay V; Pm W, Hr Pms W,; Vsw* 
Run Ft. °F. M.p.h. In. Hg Lbs. Pa;/P a, Ft. In. Hg Lbs. M.p.h, 
1 30,220 —41.0 463 61.9 7,430 1.103 28,000 68.4 8,200 466.0 
2 28,810 —33.0 467 65.9 7,370 1.143 25,800 75.4 8,430 470.0 
3 24,290 —11.0 478 79.9 7,300 1.220 19,600 96.4 8,910 480.5 
d 23,305 — 7.0 480 82.9 7,240 1.226 18,500 101.6 8,880 482.0 
5 22,850 — 4.0 480 84.5 7,110 1.242 17,700 105.0 8,830 482.0 
6 21,340 2.5 484 89.8 7,030 1.253 15,850 112.6 8,800 485.5 
Flight 130 Propeller B. Propeller-gear ratio 0.422. Exhaust stacks 2.5 sq.in. Normal gross weight 7,200 Ibs. 
Hy tay V; Pm, W, Hr Pm; W, Vsuy* 
Run Ft. °F. M.p.h. In. Hg Lbs. Pa;/Pa, Ft. In. Hg Lbs. M.p.h. 
1 30,240 —46.0 469 62.4 7,380 1.036 29,450 64.7 7,650 470.5 
2 28,750 —39.0 473 66.2 7,330 1.055 27,500 69.9 7,740 474.5 
3 24,300 —16.5 483 79.9 7,290 1.141 21,200 91.2 8,320 484.5 
4 23,300 —12.0 485 83.4 7,240 1.155 19,900 96.4 8,360 486.5 
5 22,740 —10.0 486 85.4 7,150 1.155 19,350 98.5 8,250 487.5 
6 21,310 — 5.5 488 90.0 7,060 1.145 18,100 103.0 8,100 489.0 
Flight 131 Propeller B. Propeller-gear ratio 0.422. Exhaust stacks 4.5 sq.in. Normal gross weight 7,200 Ibs. 
Hy ta; V; Pm W, Hr Pmsg W, Vsup* 
Run Ft. °F. M.p.h. In. Hg Lbs. Pas/Pa, Ft. In. Hg Lbs. M.p.h. 
1 31,210 —48.5 469 59.4 7,380 1.050 30,150 62.4 7,750 471.0 
2 30,220 —45.0 473 62.0 7,330 1.049 29,150 65.0 7,680 474.5 
3 29,240 —41.5 474 64.8 7,300 1.049 28,200 68.0 7,650 475.5 
4 28,740 —40.0 475 66.3 7,280 1.046 27,750 69.3 7,610 476.0 
5 24,280 —19.0 485 79.9 7,200 1.100 21,900 88.5 7,980 486.5 
6 22,800 —10.5 487 84.7 7,130 1.150 19,500 97.5 8,210 488.5 
7 21,310 — 2.0 491 89. 7,060 1.190 17,100 106.8 8,410 492.5 








* Vsy is the reduced true air speed corrected to the normal weight of 7,200 Ibs. 





ACp is the change in induced drag only and that the 


parasite drag is unaffected by Mach Number. Again 
neglecting second order differences: 
ACp = 2C,( AC,)/xAe (18) 


where C, corresponds to the lift coefficient at normal 
gross weight, and AC, is the increment due to weight 
difference. Then, by expressing lift coefficient in terms 
of weight, velocity in terms of Mach Number , and 
density in terms of pressure and temperature, Eqs. 
(17) and (18) can be combined and written, 


av, sW,( AW) 
V  3nAe Cp(yPa;)2M4S? 





(19) 


where the weight increment AW is positive when the 
reduced weight exceeds the normal weight W,. The ef- 
fective aspect ratio is represented by Ae, and y is the 
specific heat ratio for air. 

A study of Fig. 4 shows that the velocity-altitude 
curves are brought into closer agreement by the tem- 
perature altitude method than by the density method 
used in Fig. 3; the difference in standard velocity is 
reduced from 8 m.p.h. to 3 m.p.h. Fig. 4 also shows 
that manifold temperature and manifold pressure data 





are satisfactorily correlated by the temperature altitude 
method. 

It may be noted that two of the observed air tempera- 
tures on Flight 21 are below the minimum existing in 
standard atmosphere. Normally, such points would be 
thrown out, but for purposes of illustration the standard 
lapse rate was carried above 35,332 ft. so that these two 
points would correspond to a ‘‘standard” altitude. 

In order to show the repeatability of test data when 
air temperature conditions do not vary widely between 
flights, the data from Flights 128, 129, 130, and 131 are 
given in Table 2 and plotted in Fig. 5. No manifold 
temperature data are given for Flights 128 and subse- 
quent because these flights were made using water in- 
jection in the engine induction system. The high de- 
gree of saturation of the charge flow gave erratic read- 
ings on the manifold temperature bulb. The tempera- 
ture altitude correction method was used in reducing the 
data shown in Fig. 5. Flights 128 and 129 were made 
3 days apart with the same airplane configuration. As 
shown in the plot, one faired line will represent the two 
sets of speed vs. altitude points within 1 m.p.h. The 
propeller-gear ratio was changed from 0.479 to 0.422 
for Flight 130, and a velocity increase of 7 m.p.h. was 





Qukonwre -y 
f= | 





* 
ia 


obtai: 
For F 
from 
2 m.p 
vs. al 
four 1 
ducec 
the o 
plane 
appre 
fulne: 
Th 
on th 
airpla 
tude 
The « 


32 ¢ 


Hy, 1000 FT. 
~ U) n n $ 
n b ra) @ 
7 T 7 “ 7 


od 
°o 
+ 


@ 
7 


STD. ALTITUDE 





Fic. 
throttl 





MQWNIwoa » ” aoongjcn a * | 


se ee ee ee Sa was ee 


de 








FLIGHT-TEST PERFORMANCE DATA 551 





TABLE 3 
Experimental Fighter Airplane 
Level Flight Full-Throttle Performance Data at 3,000 R.P.M. with Water Injection 
Temperature Altitude Reduction Method 





Flight 134 Propeller C. Propeller-gear ratio 0.422. Exhaust stacks 2.5 sq.in. 


Normal gross weight 7,200 Ibs. 


Hy ‘ tay V; Pm W, Hr Pm, W, Vs * 
Run Ft. °F. M.p.h. In. Hg Lbs. Pas/Pa;, Ft. In. Hg Lbs. M.p.h. 
1 31,280 —46.5 470 59.9 7,350 1.080 29,600 64.7 7,930 472.5 
2 30,270 —48.0 467 62.6 7,330 1.010 30,000 63.2 7,400 467.5 
3 28,320 —43.0 475 68.9 7,280 0.988 28,600 68.0 7,200 475.0 
4 24,330 —43.5 476 83.1 7,240 0.822 28,750 68.3 5,940 473.5 
5 22,870 —37.0 479 88.0 7,180 0.838 26,900 73.8 6,010 477.0 
6 22,520 —35.0 478 89.4 7,140 0.846 26,350 75.6 6,040 476.0 
Flight 144 Propeller C. Propeller-gear ratio 0.422. Exhaust stacks 4.5 sq.in. Normal gross weight 7,200 Ibs. 

Hy ta; V; Pm, W, Hr Pms W; Vsu* 

Run Ft. °F. M.p.h. In. Hg Lbs. Pa,/Pa, Ft. In. Hg Lbs. M.p.h. 
1 31,210 —49.0 467 59.5 7,340 1.042 30,300 62.0 7,650 468.5 
2 30,320 —44.0 470 61.9 7,260 1.068 28,900 66.1 7,750 472.0 
3 29,260 —41.5 472 64.8 7,170 1.050 28,200 68.1 7,540 473.0 
4 25,320 —23.5 481 76.6 7,090 1.100 23,100 84.3 7,800 482.5 
5 23,300 —15.0 482 82.8 6,990 1.113 20,750 92.2 7,790 483.0 
6 21,550 —12.0 485 88.8 6,800 1.072 19,900 95.3 7,300 485.0 





* Vsy is the reduced true air speed corrected to the normal weight of 7,200 Ibs. 


obtained throughout the full-throttle altitude range. 
For Flight 131 the jet exhaust stack area was increased 
from 2.5 to 4.5 sq.in., and a further velocity increase of 
2m.p.h. was recorded. The reduced manifold pressure 
vs. altitude data show an excellent correlation for all 
four flights. From Table 2 it can be seen that the re- 
duced manifold pressures are considerably in excess of 
the observed values at the lower altitudes. The air- 
plane performance at manifold pressures beyond the 
approved limits is of importance in evaluating the use- 
fulness of emergency engine ratings. 

The effects of unusual ‘air temperature lapse rates 
on the velocity vs. altitude performance of a high-speed 
airplane as reduced by the density and temperature alti- 
tude methods are well shown in Flights 134 and 144. 
The data for these flights are given in Table 3. The 
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Fic. 5. Experimental fighter airplane. Level flight, full- 
throttle performance at 3,000 r.p.m., with water injection. 
Temperature altitude method. Normal gross weight, 7,200 Ibs. 
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Fic. 6. Experimental fighter airplane. Level flight, full- 
throttle performance at 3,000 r.p.m., with water injection. 
Density altitude method. Normal gross weight, 7,200 lbs. 
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Fic. 7. Experimental fighter airplane. Level flight, full- 
throttle performance at 3,000 r.p.m., with water injection. Tem- 
perature altitude method. Normal gross weight, 7,200 Ibs. 
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TABLE 4 
Spitfire II Fighter Airplane 
Full-Throttle Climb Performance Data at 3,000 R.P.M.* 
Temperature Altitude Reduction Method 





Temperature Altitude, Hp = 31,200 ft. 


ty = ts = — 52°F. 


Flight Hy C; V; Pm; W, Pm, W, c 
No. Ft. Ft. permin. M.p.h. In. Hg Lbs. Pa;/P a; In. Hg Lbs. Ft. per min. 
3 28,450 1,315 186 28.3 5,620 0.883 24.9 4,960 1,315 
5 26,970 1,110 186 30.1 6,480 0.825 24.8 5,340 1,110 
6 27,460 925 185 29.1 6,475 0.845 24.6 5,470 925 
7 28,940 925 186 27.7 6,460 0.900 24.9 5,810 925 
7 29,430 875 186 27.1 6,455 0.922 25.0 5,950 875 
7 29,920 830 186 26.7 6,450 0.942 % | 6,070 830 
7 30,390 775 186 26.0 6,445 0.963 25.0 6,210 775 
Temperature Altitude, Hyp = 32,300 ft. ¢, = 4, = —56°F. 
Flight Hy C; Vz Pm WwW; Pm, WwW, GS 
No. Ft. Ft. permin. M.p.h. In. Hg Lbs. Pa;/Pa;, In. Hg Lbs. Ft. per min. 
1 28,450 1,515 186 27.7 5,620 0.840 23.2 4,720 1,515 
2 28,450 1,390 186 27.9 5,620 0.840 23.4 4,720 1,390 
3 29,920 1,195 185 26.8 5,615 0.897 24.0 5,040 1,195 
5 27,950 1,030 185 29.1 6,460 0.820 23.9 5,310 1,030 
6 28,450 870 186 28.1 6,465 0.840 23.6 5,440 870 
7 32,820 335 185 23.2 6,410 1.025 23.8 6,560 335 
Temperature Altitude, Hr = 34,400 ft. ¢; = t, = —64°F. 
Flight H, V; Pm, W, Pm, W,; G 
No. Ft. Ft. per min. M.p.h. In. Hg Lbs. Pa;/Pa; In. Hg Lbs. Ft. per min. 
1 30,860 1,125 184 25.3 5,610 0.842 21.3 4,720 1,125 
1 34,860 520 184 21.2 5,570 1.017 21.6 5,670 520 
1 35,360 480 184 20.8 5, 1.042 21.6 5,800 480 
1 35,860 435 184 20.1 5,550 1.067 21.4 5,920 435 
1 36,360 405 183 19.7 5,530 1.092 21.5 6,050 405 
2 30,860 1,150 184 25.2. 5,610 0.842 21.2 4,720 1,150 
2 34,860 490 184 21.2 5,570 1.017 21.6 5,670 4 
2 35,360 425 184 20.8 5,560 1.042 21.7 5,800 425 
2 35,860 365 184 20.1 5,550 1.066 21.4 5,920 365 
2 36,360 285 183 19.7 5,540 1.092 21.5 ,060 285 
3 32,360 895 184 24.2 5,595 0.902 21.8 5,040 895 
4 é 1,485 183 29.6 5,630 0.723 21.4 4,070 1,485 
5 29,430 855 184 27.4 6,455 v.790 21.7 5,100 
6 30,390 695 184 26.2 6,445 0.825 21.6 5,310 695 
Temperature Altitude, Hy = 35,900 ft.t ¢t, = t, = —69°F. 
Flight Hy C; V; Pm, WwW, Pmys W, G 
No. Ft. Ft. permin. M.p.h. In. Hg Lbs. Pa;/Pa, In. Hg Lbs. Ft. per min. 
1 32,820 820 183 23.0 5,590 0.863 19.8 4,820 820 
3 34,860 655 183 21.2 5,570 0.951 20.1 5,300 655 
3 35,860 475 182 20.1 5,550 0.997 20.0 5,530 475 
4 28,450 1,370 183 28.5 5,620 0.707 20.1 3,970 1,370 
4 36,860 280 180 19.1 5,530 1.047 20.0 5,780 280 
5 30,390 740 183 26.2 6,450 0.771 20.2 4,970 740 
6 31,900 550 183 24.6 6,430 0.826 20.3 5,310 550 
Temperature Altitude, Hr = 38,400 ft. #4 = ¢; = —78°F.. 
Flight Hy C; V, Pm W, Pm, W, oe 
No. Ft. Ft.permin. M.p.h. In. Hg Lbs. Pas/Pa, In. Hg Lbs. Ft. per min. 
4 30,390 1,140 181 26.2 5,610 0.683 17.9 3,830 1,140 
+ 30,860 1,080 180 25.6 5,610 0.697 17.9 3,910 1,080 
4 34,340 705 181 22.1 5,575 0.823 18.2 4,590 705 
5 31,900 575 181 24.6 6,340 0.733 18.0 4,700 575 
5 32,360 530 181 24.2 6,425 0.747 18.1 4,800 530 
5 33,820 305 180 22.2 6,410 0.802 17.8 5,150 305 
5 34,340 250 181 21.8 6,390 0.822 17.9 5,250 250 
6 33,820 305 180 22.6 6,410 0.802 18.1 5,150 305 








* Take-off weight: 


5,760 Ibs. for Flights 1, 2, 3, and 4; 6,632 Ibs. for Flights 5, 6, and 7. 


t.Normal lapse rate carried above 35,332 ft. to give the “standard” altitudes shown. 


only difference in airplane configuration was the ex- 
haust stacks, which were 2.5 sq.in. on Flight 134 and 
4.5 sq.in. on Flight 144. From the results of Flights 
130 and 131 the larger stacks would be expected.to show 
approximately 2-m.p.h. increase in speed at all full- 
throttle altitudes. Upon reducing the observed data 
of Flights 134 and 144 by the density method, the rather 
startling results shown in Fig. 6 were obtained. The 
peculiar lapse rate during Flight 134 is apparent in the 
air temperature plot. The same observed data re- 


duced to temperature altitude are plotted in Fig. 7; 
although the anticipated 2-m.p.h. increment due to 
larger exhaust stacks is not shown by the reduced data, 
it is obvious that the velocity vs. altitude curves of 
Flights 134 and 144 are brought into good agreement 
by the temperature altitude method. 

A rather compiete series of flight tests has been run 
by the British on a Spitfire II fighter airplane to study 
engine performance at high altitudes.14 Data obtained 
on these flights provide an interesting check on the re- 
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Fic. 8. Spitfire II fighter airplane. Rate of climb performance 
Full throttle at 3,000 r._p.m. Temperature altitude method. 


duction of rate-of-climb performance by the tempera- 
ture altitude method. The flight trials were made with 
two different gross weights and during the course of the 
tests widely different atmospheric conditions were en- 


countered. Seven full-throttle climbs were made from 


‘ approximately 26,000-ft. pressure altitude to near 


service ceiling; four of the flights were with a take-off 
weight of 5,760 lbs., and three with 6,632 Ibs. Temper- 
ature inversions existed during several of the flights. 
In the analysis of the test data, several temperature 
altitudes were selected so as to provide the greatest 
number of test points. The flight data for five temper- 
ature altitudes are tabulated in Table 4. As in the 
case of the level flight runs made with the experimental 
fighter, in this case also air temperatures below the 
standard altitude minimum were encountered. And 
again, since these data are useful, the standard lapse 
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rate was carried beyond 35,332 ft. in order to define a 
“standard” altitude corresponding to the low air tem- 
peratures. 

In Fig. 8 the reduced rates of climb are plotted against 
reduced gross weights for the several temperature alti- 
tudes. Because of the wide variations in air pressures 
at the temperature altitudes, the reduced flight weights 
cover a wide range, and it should be noted that the 
flight data at the two test weights actually overlap and 
coincide when plotted as reduced weight vs. rate of 
climb. Therefore, it is evident that the effects of 
flight weight and altitude pressure are correctly ac- 
counted for in applying the temperature altitude 
method to the reduction of climb performance trials. 


Ii. THe Gas-TURBINE JET-PROPELLED AIRPLANE 


Engine Performance 

The performance parameters for gas-turbine jet pro- 
pulsion engines have been determined by the method 
of dimensional analysis, and the generalized perform- 
ance data on gas-turbine jets are sometimes presented 
in curves similar to those of Fig. 9.4% Shown in this 
manner, these data are independent of inlet duct pres- 
sure loss characteristics. 

As for the case of the Otto engine inlet duct, the pres- 
sure loss in the inlet ducting to the gas turbine is a 
function only of the flight and duct Mach Numbers, if 
Reynolds Numbers’ effects are taken to be small. The 
loss is due to the failure to realize full air-stream impact 
pressure at the duct inlet and to the wall friction and 
separation losses within the duct. Fig. 10a shows a 
nondimensional plot of the duct pressure loss charac- 
teristics; here the abscissa is proportional to the duct 
Mach Number and the parameter is proportional to 
the airplane flight Mach Number. By combining 
Figs. 9b and 10a, Fig. 10b is obtained, and from this 
type of plot the operating points of a given gas-turbine- 
duct unit may be determined. 

A study of Figs. 9 and 10 will show that for a given 
installation, if the three factors V, V, and f, are fixed, 
then the jet thrust, net thrust, airflow and fuel flow are 
all directly proportional to the ambient air pres- 
sure. 

Airplane Performance 

The airplane performance factors are the same for the 
gas-turbine jet-propelled airplane as for the Otto engine 
airplane with propeller. 


Flight Operating Conditions 


These will be identical to those specified in Part I. 


The ‘‘Reduced”’ Altitude 


As in Part I, the standard or ‘‘reduced”’ altitude is 
that altitude in standard air which is defined by the 
observed ambient air temperature. 
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Fics. 10a and 10b. Nondimensional representation of en- 
gine air inlet duct performance characteristics, and combined per- 
formance characteristics of inlet duct plus gas turbine. 


The ‘‘Reduced’”’ True Air Speed 


The balance between the airplane drag and enigne 
thrust forces in stabilized horizontal flight can be 
written in the usual form, 


F = CpS(p/2) V? (20) 


Eq. 20 can be written for standard air conditions (sub- 
script s) and divided by the corresponding equation for 
observed air conditions (subscript 4). Then, replac- 
ing p by P,/gRT, and eliminating constants, the follow- 
ing equation is obtained, 
F, ae i (21) 
F, Co,P,Ta,V? 
Analysis of Eq. (21) for the specified flight operating 
conditions shows that it is valid only when 


V; = V; 22 


because for this condition: 

(a) Cp, = Cp, since the standard and test flight 
Mach Numbers are equal and the angles of attack are 
equal. 

(b) F,/Pa, = F,/Pa,, since the flight Mach Number 
and the compressor tip Mach Number in standard air 
are equal to those observed during test conditions. 

Therefore, in correcting the observed flight speed of 
a gas-turbine jet-propelled airplane to standard air 
conditions by the temperature altitude method, the 
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“reduced” true air speed is equal to the observed true 
air speed. 


The ‘‘Reduced’’ Gross Weight 


As in Part I of the Analysis, the ‘‘reduced”’ gross 
weight under standard air conditions varies from the 
observed weight as the standard altitude pressure varies 
from the observed altitude pressure. That is, 


W,/W, = Pa;/Pa;, (11) 


The ‘‘Reduced’”’ Rate of Climb 
The simplified equation for rate of climb in terms of 


engine thrust and airplane drag is 
C = (F — D)V/W (23) 


Eq. (23) may be written for standard and test condi- 
tions, 
CG, (F,-—D,) V,.W 
bf eee 2 ad: :) = pail (24) 
C. (F, = D,) V, W; 
For the specified flight operating conditions previously 
defined, it was shown that F,/Pa, = F,/Pa,, D./Pa, = 
D,/Pa, V; = Vi, W;,/Pe, = W,/Pa:. 
Therefore, 


t = C. (25) 


That is, the “‘reduced”’ rate-of climb in standard air is 
equal to the observed true rate of climb under test con- 
ditions. 


Application of Method and Experimental Verification 


The detail application of the method is identical to 
that for the propeller-driven airplane with Otto cycle 
engine. No actual flight-test data on jet airplanes are 
available to the author and, therefore, no experimental 


verification is possible. 


III. THe PROPELLER-DRIVEN AIRPLANE WITH GAS- 
TURBINE ENGINE 


Engine Performance 

The generalized performance data on gas-turbine 
engines, determined by dimensional analysis methods, 
may be plotted in a manner similar to that for gas-tur- 
bine jet engines. Gas-turbine engines for airplane in- 
stallations will drive variable pitch propellers, and 
therefore a wide range of power may be absorbed at a 
givenr.p.m. Hence, an additional parameter is needed 
to completely define the engine performance charac- 
teristics. This factor is the maximum temperature 
ratio; that is, the ratio between the gas temperature 
after combustion, 7’;, to the air temperature at the com- 
pressor inlet, 7.. A set of curves similar to those of 
Fig. 9 is derived for each value of the maximum temper- 
ature ratio. Since the major portion of the useful work 
of the gas-turbine engine appears in the form of shaft 
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work, curves similar to Fig. 9a are plotted for each 
value of the maximum temperature ratio with the or- 
dinate BHP/P, in place of F,/P,. The inlet duct per- 
formance and the combined performance of the gas- 
turbine engine and duct can be as shown in Fig. 10. 

An analysis of the several performance parameters 
shows that if the four variables V, NV, T,, and T,; are 
fixed, then the shaft brake horsepower output, the ex- 
haust jet thrust, the airflow, and fuel flow are all di- 
rectly proportional to the ambient air pressure. 


Propeller and Airplane Performance 


These factors are identical to those given in Part I. 


Flight Operating Conditions 

These will be identical to the flight operating condi- 
tions specified in Part I with the addition of the maxi- 
mum gas temperature, 7’, to the fixed operating condi- 
tions. 


The ‘‘Reduced’’ Altitude, True Air Speed, Gross Weight, 
and Rate of Climb 
The respective analyses of Part I apply directly in 
these instances. 


Application of Method and Experimental Verification 


The application of the temperature altitude perform- 
ance reduction method to the case of the propeller- 
driven airplane with gas-turbine engine is identical to 
the case of the propeller-driven airplane with Otto cycle 
engine as given in Part I. 

Experimental verification of the method in this case 
is not yet possible because of lack of sufficient flight- 
test data. 


or 


REFERENCES 


1 Jones, E. T., Flight Testing Methods, Journal of The Royal 
Aeronautical Society, Vol. XLVIII, No. 401, May, 1944. 

2 Allen, E. T., Flight Testing for Performance and Stability, 
Journal of the Aeronautical Sciences, Vol. 10, No. 1, pp. 1-30, 
January, 1943. 

3 Cameron, D., British Performance Reduction Methods for 
Modern Aircraft, Report No. AAEE/Res/170, November 22, 
1942. Aeroplane and Armament Experimental Establishment, 
Boscombe Down, England. 

4Reed, A. C., Airplane Performance Testing at Altitude, 
Journal of the Aeronautical Sciences, Vol. 8, No. 4, pp. 135-150, 
February, 1941. 

5 Zeyns, J., The Air Consumption of Aircraft Engines at Alti- 
tude, D.V.L. Motortechnische Zeitschrift, November, 1939. 

® Taylor, C. F., and Taylor, E. S., The Internal Combustion 
Engine; International Textbook Company, Scranton, Pa., 
1938. 

7 Capon, R. S., and Brooke, G. V., A pplication of Dimensional 
Relationships to Air Compressors, British Aeronautical Research 
Committee R. & M. No. 1336, June, 1930. 

8 Sparrow, S. W., Relation of Fuel-Air Ratio to Engine Per- 
formance, N.A.C.A. Report No. 189, 1924. 

® Brooks, D. B., and Garlock, E. A., The Effect of Humidity 
on Engine Power at Altitude, N.A.C.A. Report No. 426, 1932. 

10 Sparrow, S. W., and Thorne, M. A., Friction of Aviation En- 
gines, N.A.C.A. Report No. 262, 1927. 

Ragazzi, Paolo, The Power of Aircraft Engines at Altitude, 
N.A.C.A. T.M. No. 895, May, 1939. 

12 Durand, W. F., Aerodynamic Theory, Division L, Chapter 1; 
Division I, Chapter 2; Julius Springer, Berlin, 1934 

13 Diehl, Walter S., Standard Atmosphere—Tables and Data. 
N.A.C.A. Report No. 218, 1925. 

14 Abbott, C. V., and Hislop, G. S., The Variation of Power 
with Height of the Merlin XII Engine as Determined by Flight 
Tests on a Spitfire IT, Report No. AAEE/Res/166, May 26, 1942. 
Aeroplane and Armament Experimental Establishment, Bos- 
combe Down, England. 

15 Warner, D. F., and Auyer, E. L., Contemporary Jet-Propul- 
ston Gas Turbines for Aircraft, Mechanical Engineering, Vol. 67, 
No. 11, November, 1945. 


Letters to the Editor 


(Continued from page 536) 


(6) In the discussion of the effects of pressurization on the 
stringers and frames, no account is taken of the fact that many 
stringers and frames are critically loaded by inward forces be- 
cause of tension field formation. The stabilizing effect of pres- 
surization would be extremely important in determining the de- 
sign conditions for the stringers and frames. 
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A. H. FLtax 
Piasecki Hellicopter Corporation 


Dear Sir: 

Since my paper ‘“‘Analysis of Relay Servomechanisms”’ (Jour- 
NAL OF THE AERONAUTICAL SCIENCES, Vol. 13, No. 7, p. 364, 
July, 1946) was submitted, two excellent references not men- 
tioned in my paper have been brought to my attention. Servo- 
mechanisms, by*L. A. MacColl (D. Van Nostrand, 1945), indi- 
cates the application of graphical analysis to a simple relay servo- 
mechanism. Introduction to Non-Linear Mechanics, Part I, 
“Topological Method of Non-Linear Mechanics,”’ by N. Minor- 
sky (David Taylor Model Basin Report 534, December, 1944), 
has been declassified to permit circulation without restrictions 
as to security. The latter paper is an excellent presentation of 
‘phase plane”’ studies of nonlinear systems, describing, in addi- 
tion, much Russian work in the field which is not available to the 
average engineer. It is believed that these references should be 
brought to the attention of those interested in the analysis of 
nonlinear servomechanisms. 


HERBERT K. WEIsS 
Ballistic Research Laboratory 
Aberdeen Proving Ground 
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AERONAUTICAL SCIENCES and the AERONAUTICAL ENGI- 
NEERING REvIEW. The Institute, following the prac- 
tice of other societies, does not pay for contributions. 


The following directions for the preparation of papers, 
if followed by authors, will save correspondence, avoid 
the return of papers for changes, minimize the work of 
preparation for the printer, and save the expense due 
to the charges made for ‘‘author’s corrections.” 
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papers submitted be in final form. Typographical errors may 
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rolled, but manuscripts should be sent flat. Send by first class 
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torial Office, Institute of the Aeronautical Sciences, 2 East 64th 
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quent duplication of initials and surnames which sometimes 
makes it difficult to establish the identity of the author. This is 
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services. All titles and degrees or honors are omitted. The name 
of the organization with which the author is associated should be 
placed after his name on a separate line: The date on which the 
paper is received will be inserted by the Editor. 


SUMMARIES OR AssrRacts: An abstract to be printed at the 
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quate as an index and as a summary. It should contain a state- 
ment of major conclusions reached, since summaries in many 
cases constitute the only source of information used in compiling 
scientific reference indexes. Abstracts printed in other journals, 
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possible the major conclusions to a nonspecialist in the subject. 
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on the length of the paper. 


SUBHBADINGS: Subheadings should be inserted by the au- 
thor at frequent intervals. The work of editorial preparation 
will be simplified by the author providing many subheadings. 
Because of the breaking of columns and the insertion of illustra- 
tions, some of the subheadings may have to be omitted. 


SHORTENING OF PAPERS: Some papers, at the end, fill in only 
a portion of a page. This leaves much wasted blank space as 
succeeding articles are started at the top of a.page. Authors 
should indicate by notation on the left-hand side of the page 
what matter may be omitted when “runovers” occur. This 
request is important as the Journal cannot afford in the future, 
as it has in its earlier issues, to have blank half pages or more at 
the end of papers. 


MATTER USUALLY DELETED: Acknowledgments of assistance 
in preparation of paper, except by collaborators. Photographs 
or illustrations of little technical interest and not showing ad- 
vances in general practice. Too detailed tabular matter (gen- 
eral results of such tables may be included in the text). Lengthy 
descriptions of materials or processes or of preliminary experi- 
ments or theories that preceded final results; salient features 
only are of interest. 


REFERENCES AND FOOTNOTES: References should appear as 
footnotes only, numbered consecutively, grouped together at the 
end of the manuscript. The arrangement should be as follows: 
(for books)—! Durand, W. F., Aerodynamic Theory, 1st Ed., 
Vol. 1, p. 23; Julius Springer, Berlin, 1934. (For magazines)— 


1 Englund, C. R., Crawford, A. B., and Mumford, W. W., Some 
Results of a Study of Ultra-Short-Wave Transmission Phenomenon, 
Proc. I. R. E., Vol. 20, No. 12, pp. 481 and 482, March, 1933. 
Please give author, title, edition, volume, page, publisher, and 
date of publication as indicated. Omission of one required fact 
causes much extra editorial work and possible inaccuracies. All 
references should be grouped at the end of the article. 


ILLUSTRATIONS: Illustrations should accompany manuscripts 
and each should always be referred to in the text, preferably by 
number. Drawings or graphs should not be larger than 12 X 16 
inches, and must be made with jet black India Ink on white 
paper or tracing cloth, the latter being preferred. Do not use 
typewriter for lettering. The smallest lettering on 8 X 10 inch 
figures should be no less than !/, inch high. Cross-section paper 
(white with black lines) may be used, but should not have more 
than 4 lines per inch. If finer ruled paper is used, the major 
division lines should be drawn in with black ink, omitting the 
finer divisions. In the case of finely ruled paper, only blue- 
lined paper can be accepted. Tracing paper and blueprints are 
not acceptable. Lettering and all markings must be large enough 
to be readable after reduction. Mail rolled or flat, never fold. 
Drawings that cannot be reproduced (including pencil draw- 
ings) will be returned to the author for redrawing, thus delaying 
publication of the paper. Photographs should be distinct and 
show clear black and white contrasts. They must be on glossy 
white paper. Avoid round and oval photographs. 


CAPTIONS AND LEGENDS: Legends or captions must accom- 
pany each drawing or photograph submitted. If written on the 
drawing or photograph, they should be placed below and well out- 
side the part to be reproduced. It is better to place them on 
separate sheets of paper pasted to the back of the drawings or 
photographs. Each table should have a caption such as Table 1, 
Table 2, Table 3, etc. Captions should be complete in them- 
selves so as to make the data intelligible to the reader without 
reference to the text. A duplicate list of captions for figures 
should be included as the last page of the manuscript. Use “‘Fig. 
1” (not Figure 1), Figs. 3 and 4, etc., in both the text and the 
numbering of illustrations. In the text, ‘Eq. (1),’’ or ‘“‘Egs. (1) 
and (2)” are preferable to ‘Equation (1).’’ In captions and 
legends, except for ‘‘Fig.”’ and ‘‘Eq.,”’ and in table headings, write 
all words in full; do not abbreviate. Avoid placing explanatory 
written matter in the drawings; it should be in the text. 


MATHEMATICAL WorK: Only the simplest formulas should 
be typewritten; all others should be carefully written in pen 
and ink, the writing to be large enough so that ample room is 
provided to mark mathematical matter for the printer. A con- 
siderable space for marking should be allowed above and below 
all equations. All complicated equations should be repeated on 
separate sheets with plenty of space left for marking. The solidus 
should be used for simple fractions appearing within the text. 
Make all expressions clear to the typesetter. Greek letters used 
in formulas should be clearly designated by name on the margin 
of the manuscript. All symbols should be clearly written and 
carefully checked. The difference between capital and lower- 
case letters should be clearly distinguished and care taken to avoid 
confusion between zero (0) and the letter (0), between the numeral 
(one) and the letter (ell) and the prime (’), between alpha and a, 
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NOMENCLATURE AND ABBREVIATIONS: The National Advisory 
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PRATT & WHITNEY AIRCRAFT DIVISION 
PRATT & WHITNEY AIRCRAFT CORPORATION OF 
MISSOURI 
SIKORSKY AIRCRAFT DIVISION 
UNITED AIR LINES, INC. 
UNITED STATES AVIATION UNDERWRITERS, INC 
UNITED STATES RUBBER COMPANY 
DOMINION RUBBER COMPANY, LTD. 
THE VARIETY AIRCRAFT CORPORATION 
VICKERS, INC. 
VIDAL CORPORATION 
THE WACO AIRCRAFT COMPANY 
WARNER AIRCRAFT CORPORATION 
THE WEATHERHEAD COMPANY 
WESTERN AIR LINES, INC. 
WESTINGHOUSE ELECTRIC CORPORATION 
WESTON ELECTRICAL INSTRUMENT CORPORATION 
WYMAN-GORDON COMPANY 
YOUNG RADIATOR COMPANY 





























